Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,856,594 A→G H212H (CAT→CAC)  unknown ← dTDP‑4‑dehydrorhamnose reductase (EC 1.1.1.133)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,856,5940AG97.5% 114.9 / NA 40H212H (CAT→CAC) unknowndTDP‑4‑dehydrorhamnose reductase (EC 1.1.1.133)
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/14);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTACTCCGCATACGCTTTAATATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAATGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGCTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACAACGGCTTGTTCGC  >  1/1856348‑1856835
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
tCTTTACTCCGCATACGCTTTAATATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑t                                                                                                                                                                                                                                               >  1:763813/1‑251 (MQ=255)
tCTTTACTCCGCATACGCTTTAATATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑t                                                                                                                                                                                                                                               >  1:420650/1‑251 (MQ=255)
          gCATACGCTTTAATATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTACGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCata                                                                                                                                                                                                                                     >  1:444848/1‑251 (MQ=255)
                 cTTTAATATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAAc                                                                                                                                                                                                                               <  1:321695/250‑1 (MQ=255)
                     aaTATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGAtt                                                                                                                                                                                                                          <  1:466962/251‑1 (MQ=255)
                             tCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGa                                                                                                                                                                                                                                                 <  1:186189/221‑1 (MQ=255)
                                       ttGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGAtttttt                                                                                                                                                                                                        >  1:314757/1‑251 (MQ=255)
                                                     gttggATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAg                                                                                                                                                                                                                             <  1:731948/216‑1 (MQ=255)
                                                     gttggATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAAtt                                                                                                                                                                                          >  1:667366/1‑251 (MQ=255)
                                                     gttggATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAAtt                                                                                                                                                                                          >  1:669294/1‑251 (MQ=255)
                                                            aCCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGc                                                                                                                                                                                     >  1:776793/1‑249 (MQ=255)
                                                                    cTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAata                                                                                                                                                                           <  1:390843/251‑1 (MQ=255)
                                                                       gttgAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAatat                                                                                                                                                                          >  1:122634/1‑249 (MQ=255)
                                                                        ttgAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATAtct                                                                                                                                                                        <  1:620107/250‑1 (MQ=255)
                                                                                                 gAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTccgcaaggtat                                                                                                                                              >  1:415138/1‑244 (MQ=255)
                                                                                                                       gtgtTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAgtg                                                                                                                        >  1:753913/1‑251 (MQ=255)
                                                                                                                       gtgtTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAgtg                                                                                                                        >  1:856669/1‑251 (MQ=255)
                                                                                                                            ggATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACg                                                                                                                   >  1:698452/1‑251 (MQ=255)
                                                                                                                             gATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACg                                                                                                                   <  1:598267/250‑1 (MQ=255)
                                                                                                                                   ccgccgTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTg                                                                                                            >  1:646878/1‑251 (MQ=255)
                                                                                                                                            tAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACgc                                                                                                   <  1:666264/251‑1 (MQ=255)
                                                                                                                                                      gACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAg                                                                                                                                                                                                                             >  1:100394/1‑119 (MQ=255)
                                                                                                                                                                     ttCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCAt                                                                                                                                                                 >  1:484823/1‑164 (MQ=255)
                                                                                                                                                                     ttCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAggtc                                                                                                                                                            >  1:467656/1‑169 (MQ=255)
                                                                                                                                                                     ttCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATc                                                                                                                                           >  1:300302/1‑186 (MQ=255)
                                                                                                                                                                         aaCATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTg                                                                      >  1:77229/1‑251 (MQ=255)
                                                                                                                                                                           cATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATAtct                                                                                                                                                                        >  1:816517/1‑151 (MQ=255)
                                                                                                                                                                                 ttGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGa                                                              <  1:357289/251‑1 (MQ=255)
                                                                                                                                                                                    cGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAAt                                                            <  1:502593/250‑1 (MQ=255)
                                                                                                                                                                                       aaccgctttattcaaaaTTGGTCCGGCAAAATCCCCCCAGGTTGCATTAGGTTCCCCCGTAAAATGATTAAAT‑CCAAATTCAACAGGTTCGGCTGGAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATAccc                                                        <  1:459009/238‑1 (MQ=255)
                                                                                                                                                                                           gCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAg                                                     >  1:26599/1‑250 (MQ=255)
                                                                                                                                                                                                atcaAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGCTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAg                                                                                                                                                                                                                             >  1:426639/1‑77 (MQ=255)
                                                                                                                                                                                                          gCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGAt                                     >  1:535900/1‑251 (MQ=255)
                                                                                                                                                                                                                gCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGATAAAtt                                <  1:757602/250‑1 (MQ=255)
                                                                                                                                                                                                                              gCTGACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGc                                             <  1:496384/223‑1 (MQ=255)
                                                                                                                                                                                                                                 gACATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACaa              >  1:752537/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aCATAAGGTTCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACaa              <  1:523689/250‑1 (MQ=255)
                                                                                                                                                                                                                                           tCCCCCGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACAACGGCTTGTTc    <  1:766034/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cccGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACAACGGCTTGTTCGc  >  1:406747/1‑250 (MQ=255)
                                                                                                                                                                                                                                              cccGTAAAGTGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGTTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACAACGGCTTGTTCGc  >  1:539030/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
TCTTTACTCCGCATACGCTTTAATATTTTTCAAAGCTCTTTGCCAATCACTTGGTTGGATACCGAATGCTTGTTGAATTTTGGTCAAATCTAAGCAAGAATTCGCAGGGCGTTTGGCTGGTGTTGGATAATCCGCCGTGATAATTGCATTGACAAGCGGTGCTTTTTCTAACATATTTTGCGAAACAGCTTCATCAAAAATTGCTCTGGCAAAATCACACCAGCTGACATAAGGTTCCCCCGTAAAATGA‑TAAATACCATATTCAACAGATTCGCCTGCAATGATTTTTTCTGCAATTTGAATTAACGCAGCGGCAATATCTCCTGCATAGGTCGGTCCGCCGATCTGATCAGCAACTACGCCTAAAGTGTCACGAGTTTTGGCTAAACGCAGCATTGTTTTAACGAAATTGCTGCCGTGCTCGCAGAATACCCAAGCGGTACGCAATACGATAAATTTATCGTTTGCTTCTACAACGGCTTGTTCGC  >  1/1856348‑1856835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑