Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,896,621 A→G D147D (GAT→GAC bglB ← 6‑phospho‑beta‑glucosidase BglB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,896,6210AG98.2% 370.2 / ‑5.1 111D147D (GAT→GACbglB6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (75/34);  total (76/35)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 111 total aligned reads displayed.

CATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896537‑3896707
                                                                                    |                                                                                      
cATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGAt                                                                                   >  1:1912017/1‑90 (MQ=255)
 aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  <  1:677534/90‑1 (MQ=255)
 aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  >  2:1175837/1‑90 (MQ=255)
 aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  >  1:3025762/1‑90 (MQ=255)
 aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  <  2:1357699/90‑1 (MQ=255)
 aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  <  2:2229228/90‑1 (MQ=255)
   gTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAc                                                                                >  1:352150/1‑90 (MQ=255)
      gATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGTGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  1:2382030/90‑1 (MQ=255)
      gATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  1:2346502/90‑1 (MQ=255)
       aTTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCt                                                                            >  1:76436/1‑90 (MQ=255)
           cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                                        <  2:2165373/90‑1 (MQ=255)
            aTTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAtt                                                                       <  1:1162650/90‑1 (MQ=255)
                    tCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAc                                                               >  2:948967/1‑90 (MQ=255)
                     cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:1662031/1‑90 (MQ=255)
                           tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             >  2:1401010/1‑69 (MQ=255)
                           tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  1:1401010/69‑1 (MQ=255)
                           tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTa                                                        >  1:1593377/1‑90 (MQ=255)
                            aaCGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAg                                                       >  2:1629887/1‑90 (MQ=255)
                             aCGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGt                                                      >  2:2679245/1‑90 (MQ=255)
                                 cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                                                  >  2:1128903/1‑90 (MQ=255)
                                 cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                                                  >  2:1678859/1‑90 (MQ=255)
                                 cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                                                  >  2:2569435/1‑90 (MQ=255)
                                  aCTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTc                                                 >  1:2313621/1‑90 (MQ=255)
                                     ttttgtgGGTAGCGAATAAAGCCGGGGGGGGCGAATTGCTGGAAGGGGCCGATGCCCGCTCGTTTACCCCACCCGCCGTAGTTTTTCAcc                                              <  1:1956337/87‑1 (MQ=255)
                                     ttATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCAcc                                              >  2:2516421/1‑90 (MQ=255)
                                      tttGTTGGTAGCGAGTAAACCCCGGGCGGGCGTAATGCCCGAGGTGGTCGATGCCCGCCGGATTACCCCAACCGCCGTAGTTTTTCACCa                                             <  2:1527292/88‑1 (MQ=255)
                                      tATGTTGGTAGGGAGTAAAGCCCGGGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGTTTACCCCAACCGCCGTAGTTTTTCACCa                                             <  2:432177/90‑1 (MQ=255)
                                      tATGTAGGTAGCGAGAAAAGCCCGTGGGGGGGAAATCCTCAAAGGGGGCGATCCCCGCGGGTTTACCCCAACCGCCGTAGTTTTTCACCa                                             <  2:1499251/90‑1 (MQ=255)
                                       tttttGGTAGGGAGAAAAGCCCGGCCGGGCGTAATGCCAGAAGTGGTCGTTGCCCGCTCGATTCCCCCAACCGCCGTAGTTTTTCACCAg                                            <  2:2089479/87‑1 (MQ=255)
                                       aTGTTGGTAGCGAGAAAAGCCCGTGCGGGCGTAATGCCCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAg                                            <  2:21510/90‑1 (MQ=255)
                                        tGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGc                                           >  1:97927/1‑90 (MQ=255)
                                         gTTGGTAGCGAAAAAAGACCGGGGGGGTGTAAGCCCCGAAGGGGTCGATGCCCGCTCTATTAGCCCAACCGCCGTAGTTTTTCACCAGcc                                          <  2:2026418/90‑1 (MQ=255)
                                            ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:310972/1‑90 (MQ=255)
                                            ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:1993595/1‑90 (MQ=255)
                                            ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:1088066/1‑90 (MQ=255)
                                             gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata                                      >  2:972527/1‑90 (MQ=255)
                                             gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata                                      >  2:463380/1‑90 (MQ=255)
                                              tAGCGAGAAAAGCCCGTGCGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCatat                                     <  2:315650/90‑1 (MQ=255)
                                               aGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATg                                    >  2:1264886/1‑90 (MQ=255)
                                                 cGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGc                                  >  1:67221/1‑90 (MQ=255)
                                                 cGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGc                                  >  1:1813051/1‑90 (MQ=255)
                                                 cGAGAAAACCCCGGGCGGGCGAAATGCTGGAGGGGGCGGATCCCCGCTGGATAACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGc                                  <  1:2461790/90‑1 (MQ=255)
                                                  gaaaaaaGCCCGGGGGGGGGTAATGCCCAAGGTGGTCGATGCCCCCTCAATTACCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCa                                 <  2:2886321/86‑1 (MQ=255)
                                                   ggTAAAGCCTGGGCGGGCGTAATCCCCGAAGTGGTCGAGGCCCGCTGGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCAt                                <  2:795805/89‑1 (MQ=255)
                                                   agaaaaGCCCGGGCGGGCGTAATGCCCGAGGGGGCCGATGCCCCCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCAt                                <  2:1839598/87‑1 (MQ=255)
                                                   agaaaaGCCCGGGCGGGCGAAATGCTCGAAGGGGTCGATGCCCGCCCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCAt                                <  2:954621/87‑1 (MQ=255)
                                                    gTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCAtt                               >  1:1578596/1‑90 (MQ=255)
                                                          cccGGGCGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTaa                         <  2:873864/89‑1 (MQ=255)
                                                          aCCGGGCGGGCGTAATGCTCGAAGTGGTGGATGACCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTaa                         <  1:459245/90‑1 (MQ=255)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  1:2957/36‑1 (MQ=38)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             >  2:2957/1‑36 (MQ=38)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:569740/1‑90 (MQ=255)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:2754072/1‑90 (MQ=255)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:2215807/1‑90 (MQ=255)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:1861564/1‑90 (MQ=255)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:180006/1‑90 (MQ=255)
                                                            cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:2778087/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             >  1:2952196/1‑35 (MQ=38)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  2:2952196/35‑1 (MQ=38)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGTCCATATGGCATTTCGTAATgg                      >  2:544709/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTTGTAAggg                      >  1:2378068/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:585631/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:1130837/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:1119360/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:2671188/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:1336775/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:678801/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:1405025/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:722234/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:2376070/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:807941/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:141456/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:530931/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:2052659/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:393467/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:2160030/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:1325624/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:2227057/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:1372511/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:1464926/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:1380757/1‑90 (MQ=255)
                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAACCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:2101997/1‑90 (MQ=255)
                                                              tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:395073/1‑90 (MQ=255)
                                                              tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:1105344/1‑90 (MQ=255)
                                                              tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:1151369/1‑90 (MQ=255)
                                                              tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  1:2275110/1‑90 (MQ=255)
                                                              tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggc                     >  1:795805/1‑89 (MQ=255)
                                                               gCGGGGGTAATGCTGGAAGGGGCCGATCCCCGCTCGATTACCCAACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa                    <  1:2467512/90‑1 (MQ=255)
                                                               gCGGGCGTAATGCTCGAGGGGGTCGATGCCCCCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa                    <  2:28560/90‑1 (MQ=255)
                                                                 gggCTTAATGCTGGAAGTGATCGATCCCCGCTCGTTTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa                  <  1:2071320/90‑1 (MQ=255)
                                                                   gCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAc                >  1:979432/1‑90 (MQ=255)
                                                                   gCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAc                >  1:1548189/1‑90 (MQ=255)
                                                                            tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc       >  2:1474809/1‑90 (MQ=255)
                                                                            tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc       >  1:1161689/1‑90 (MQ=255)
                                                                             ggaggTGGTCGATGCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCg      <  1:313323/86‑1 (MQ=255)
                                                                             cGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCg      >  1:1056546/1‑90 (MQ=255)
                                                                               aaGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTa                                                        >  2:89345/1‑38 (MQ=38)
                                                                               aaGTGGTCGATGACCGCTCGATTACCCCAACCGCCGTa                                                        <  1:89345/38‑1 (MQ=38)
                                                                                 gTGGTCGATGCCCGCTCTATTCCCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt  <  1:948967/90‑1 (MQ=255)
                                                                                 gTGGTCGATCCCCCCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt  <  2:1130837/90‑1 (MQ=255)
                                                                                    |                                                                                      
CATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896537‑3896707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: