Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,228 G→C 14.0% intergenic (‑189/+179) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,2280GC14.0% 469.8 / 22.0 156intergenic (‑189/+179)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base G (45/87);  new base C (19/3);  total (66/90)
Fisher's exact test for biased strand distribution p-value = 6.67e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

AGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  NC_000913/1351082‑1351376
                                                                                                                                                  |                                                                                                                                                    
aGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg                                                                                                                                                   <  2:792031/150‑1 (MQ=255)
aGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg                                                                                                                                                   <  2:2715762/150‑1 (MQ=255)
 ggCGATTCGTCGTTTTATTAAAAAGTACGAACAGATAAAAGGTTATTATAAACAACCTGGCTGTGAGTAGCTATAGTTTCCAGGTACGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTACCCAAACCCAAAACCAAACCCCACAggg                                                                                                                                                  >  1:448298/1‑150 (MQ=255)
 ggCGATTCGTCGTTTTAGTAAACATTAAGAACAGATAAAACGTTATTATAAAAAACCAGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCCGAGCAGGCTGCGGCAGGGGGGGCGCGTTTTCCCACACCCAAAACCAATCACCACAggg                                                                                                                                                  >  1:5425/1‑150 (MQ=255)
  gCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                                                                                                                                 <  1:2631/150‑1 (MQ=255)
   cGATTCGTCGTTTTAGTAAACAGTAAGAACAGATAAACGGTTATTATAATCAAACAGGCTGTGAGTAGCTATAGTTTCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTGCCCACACCCAAAACCACACCCCCCAGGGGc                                                                                                                                                >  1:877459/1‑150 (MQ=255)
    gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCGAACCCGAGCCCCACAGGGGCg                                                                                                                                               >  2:2364231/1‑150 (MQ=255)
    gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTCCCCCCACCCAAACCCCCACCCCCCAGGGGcc                                                                                                                                               >  2:1234490/1‑149 (MQ=255)
    gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTACCCCCCCCCCAACCCAAAACCCAGGGGGGCg                                                                                                                                               >  2:700853/1‑150 (MQ=255)
    gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCGATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTGCCCCCCCCCAAAACCACACCCCCGCgggggg                                                                                                                                               >  2:700855/1‑148 (MQ=255)
      ttCGTCGTTTTAGTAAACAGTACGAACAGATACACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  2:1554779/150‑1 (MQ=255)
      ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATACTCAACCTGGCTGTGAGTAGCTATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  2:1088323/150‑1 (MQ=255)
      ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  2:1732492/150‑1 (MQ=255)
      ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  2:1364377/150‑1 (MQ=255)
      ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  1:1781641/150‑1 (MQ=255)
      gtcgtcgTTTTAGTACACAGTACGAACAGATAACCGGTTATGATCATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                             <  1:166757/149‑1 (MQ=255)
       tcgtcgTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                                                                                                                            <  1:918861/150‑1 (MQ=255)
          tcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGa                                                                                                                                         <  1:2576831/150‑1 (MQ=255)
             ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttagg                                                                                                                                      <  1:1076030/150‑6 (MQ=255)
             ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttagg                                                                                                                                      <  2:2358347/150‑6 (MQ=255)
                    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCACACCCCAGAGGGGCGAGGGGCCCGTATtggg                                                                                                                               >  2:1074796/1‑148 (MQ=255)
                    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCTAACCCTATCCCCAGAGGGGCGAGGGGACCGTATtggg                                                                                                                               >  2:629701/1‑148 (MQ=255)
                    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCGAACCCTCACCCCAGAGGGGCGAGGGGACCGTATtgtg                                                                                                                               >  2:1982933/1‑150 (MQ=255)
                    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCCAAACCACACCCCCGACGGGCGGGGGGGCCCTAGtgtg                                                                                                                               >  2:1856258/1‑150 (MQ=255)
                    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAACCCGCTCCCCAGAGGGGCGAGGGGACCGGATtgtg                                                                                                                               >  2:306351/1‑150 (MQ=255)
                    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCACACCCAAACCCGCACCCCAGAGGGGCGAGGGGACCCTATtgtg                                                                                                                               >  2:40729/1‑150 (MQ=255)
                           aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                                                                       <  1:2477562/150‑1 (MQ=255)
                            cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                                                                       <  1:2348966/150‑1 (MQ=255)
                            cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTGCCCAACCCCAAAACCACAACCCCCACGGGCGAGGGGCCCCGATTTTTGAAATAtt                                                                                                                       >  1:13850/1‑150 (MQ=255)
                            cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTCCCCACCCCCAAACCCACGCCCCCCCCGGGCGCGGGGCCCGGAGTGTGCcaaaaat                                                                                                                       >  1:1893226/1‑146 (MQ=255)
                            cGAACAGATAAACGGTTATTATAATCAACCTGGCGGTGAGTAGCTATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                                                                       <  1:241752/150‑1 (MQ=255)
                             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAAACCTCAACCCAGAGGGGCGAGGGGACCGTATTGTGGAAATAttt                                                                                                                      >  1:352460/1‑149 (MQ=255)
                              aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGt                                                                                                                     <  1:1276936/150‑1 (MQ=255)
                              aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGt                                                                                                                     <  1:2122112/150‑1 (MQ=255)
                              aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGt                                                                                                                     <  1:2307390/150‑1 (MQ=255)
                               aCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt                                                                                                                    <  2:923886/150‑1 (MQ=255)
                                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTGCCCACACCCAAAACCACAACCCACCGGGGCGAGGGGGCCGGATTTTGGAAATATtttttac                                                                                                                 >  1:786940/1‑147 (MQ=255)
                                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCTAACCCACACCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcc                                                                                                                 >  1:2299111/1‑150 (MQ=255)
                                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAAACCAAAACCCACAGGGGCGAGGGGGCCCGATTGTGCAAATATTGTTAcc                                                                                                                 >  1:2164508/1‑150 (MQ=255)
                                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTGTTCCCCTCACCCAAACCCACTCCCCAGAGGGGCGAGGGGCCCGGATTGTGGAAATATTGTTaaa                                                                                                                 >  1:2159678/1‑148 (MQ=255)
                                     aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa                                                                                                              <  1:1658320/150‑1 (MQ=255)
                                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaaca                                                                                                    <  2:2486571/150‑1 (MQ=255)
                                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaaca                                                                                                    <  1:2343409/150‑1 (MQ=255)
                                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaaca                                                                                                    <  1:177607/150‑1 (MQ=255)
                                                ataATCAACCTGGCTGTGAGTAGCTATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTACCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAg                                                                                                   <  2:1544319/150‑1 (MQ=255)
                                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAAACCACTCCCCAGAGGGGCGAGGGGACCCTATTGTGCAAATATTTTTACCCCAGCAACAAAAAGGCt                                                                                                >  2:2683720/1‑150 (MQ=39)
                                                   aTCAACCTGGCTGTGAGTAGCTATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCt                                                                                                <  2:2320939/150‑1 (MQ=255)
                                                      aaCCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCACACCCCACAGGGGCGAGGGGCCCCGATTGTGCAAATAATGTTAACCCCCCAACAAACACGcacaa                                                                                             >  2:2126642/1‑149 (MQ=255)
                                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACa                                                                                          <  2:341520/150‑1 (MQ=255)
                                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACa                                                                                          <  1:329528/150‑1 (MQ=255)
                                                          tGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAg                                                                                         <  2:535552/150‑1 (MQ=255)
                                                              tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGcccc                                                                                     <  1:1493094/150‑1 (MQ=255)
                                                                 gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTaa                                                                                  <  2:808642/150‑1 (MQ=255)
                                                                    tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAccc                                                                               <  1:335268/150‑1 (MQ=255)
                                                                    tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAccc                                                                               <  2:2196619/150‑1 (MQ=255)
                                                                    tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATAGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAccc                                                                               <  2:2558698/150‑1 (MQ=255)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                                                            <  1:1201677/150‑1 (MQ=255)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                                                            <  1:172983/150‑1 (MQ=255)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCTCTCCCCCGAGGGGCGAGGGGACCGTATTGTGCCAATATTGTTACCCCAGCAACAAACCGGCGCATACAGCCCCGAAcccggt                                                                            >  2:196104/1‑147 (MQ=255)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCGAACCCACTCCCCAGAGGGGCGAGGGGACCCGATTGTGCAAAAATTGTTACCCCAGCAACAAACACGCTCAGACAGCCCCGAAaccaat                                                                            >  2:2351972/1‑147 (MQ=38)
                                                                       cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTCCCCTCACCCAAACCCACTCCCCAGAGGGGCGAGGGGACCCGAGGGTTCAAAAATTGGTACCCCAGCCACACACAGGCGCAGACAGCCCCAAACCCctt                                                                            >  2:1331088/1‑150 (MQ=25)
                                                                             ttGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCAAACCCTTTCATGGc                                                                      >  2:2296322/1‑150 (MQ=255)
                                                                             ttGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTACCCGCACCCCAACCCACGACCCCCAGGGGCGGGGGGGCCCGATTGTGCCCCTATTGTTAACCCAGCCACAAACAGGCACAAACACCCCCAAACCCGTTTATGGc                                                                      >  1:575532/1‑150 (MQ=255)
                                                                                      ccGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCGCAACCGAAACCACTCCCCAGAGGGGCGCGGGGTCCCTATTGTTCAAATATTTTTAACCCAGCAACAAAAAGGCTCATAAAGCCCCAAACCCCTTGAAGGCGATGGCggg                                                             >  2:2573279/1‑150 (MQ=255)
                                                                                      ccGACCGGAGCAGGCTGCGGCAGGCGGGGCGCTTAAAACCTCACCCAAAAAAACTACCCAAAGGGGCGACGGGACCGTAATTAGAAAATATAGATAAACAACCACAAAACAGGATAATACAGACCCTAACCCTATAATTGCGATAGCTgg                                                             >  2:1907547/1‑150 (MQ=255)
                                                                                                cAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAg                                                   <  2:1393465/150‑1 (MQ=255)
                                                                                                  ggCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGAc                                                 <  2:1404097/150‑1 (MQ=255)
                                                                                                    cTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCt                                               <  1:2539570/150‑1 (MQ=255)
                                                                                                     tGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                              <  2:1602360/150‑1 (MQ=255)
                                                                                                     tGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                              <  2:2119056/150‑1 (MQ=255)
                                                                                                       cGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc                                            <  1:100218/150‑1 (MQ=255)
                                                                                                       cGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc                                            <  1:2339155/150‑1 (MQ=255)
                                                                                                       cGGCAGGGGGGGCGCTTTTCCCCTCACCCGAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc                                            <  1:1951106/150‑1 (MQ=255)
                                                                                                        ggCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGcc                                           <  1:618632/150‑1 (MQ=255)
                                                                                                             ggggggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAata                                      <  1:538497/150‑1 (MQ=255)
                                                                                                             ggggggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAata                                      <  1:153345/150‑1 (MQ=255)
                                                                                                              gggggCGCTTTTACCCTCACCCTAAACCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTTTTACCCCAGCAACAAACAGGCTAAAAAAGCCCCAAACCCATTCATGGCCATTGCTGGGACCGTTCAGACCATGCCGAatat                                     >  2:1480078/1‑150 (MQ=37)
                                                                                                               ggggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATAtt                                    <  1:1648984/150‑1 (MQ=255)
                                                                                                                 ggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAGTGTGCAAATATTGTTACCCCAGCAACAAACAGGCGCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtct                                  >  2:2660029/1‑150 (MQ=255)
                                                                                                                 ggCGCTTTTCCCCTCACCCTAACCCTCACCCCAGAGGGGCGAGGGGACCCTATTGTGCAAATATTTTTAACCCAGCAACAAAAAGGCTCATACAGCCCCTCACCCTTTCATGGCCATGGCTGGGACGGTTCAGACCTTGCCGAATATtct                                  >  2:2326751/1‑150 (MQ=38)
                                                                                                                      ttttCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGc                             <  1:606144/150‑1 (MQ=255)
                                                                                                                      ttttCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGc                             <  2:1322402/150‑1 (MQ=255)
                                                                                                                      ttttCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGc                             >  2:1564702/1‑150 (MQ=255)
                                                                                                                         tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCAcc                          <  2:2299111/150‑1 (MQ=255)
                                                                                                                         tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCAcc                          <  2:472233/150‑1 (MQ=255)
                                                                                                                         tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCAcc                          <  2:2164508/150‑1 (MQ=255)
                                                                                                                         tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCAcc                          <  2:2299152/150‑1 (MQ=255)
                                                                                                                                aCCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACACACAGGCTCATACAGCCCCTAACCCTGTCATGGCGATGGCTGGGCCGGTTCAGACCGTGCCGACTAGTCTCCAGCACCGTCGccc                   >  1:2508496/1‑149 (MQ=38)
                                                                                                                                aCCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCGCCa                   >  1:1699267/1‑150 (MQ=255)
                                                                                                                                   cTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGt                <  2:1490730/150‑1 (MQ=255)
                                                                                                                                        ccTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           <  1:1839854/150‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGTCCCTAACCCTTTCATgg                                                                       <  2:396124/89‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGTCCCTAACCCTTTCATgg                                                                       >  1:396124/1‑89 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTACCCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          <  2:5425/150‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          <  2:448298/150‑1 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          >  1:723818/1‑150 (MQ=255)
                                                                                                                                         ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          >  1:1109191/1‑150 (MQ=255)
                                                                                                                                             ccccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACcaca      >  2:2469182/1‑150 (MQ=255)
                                                                                                                                               ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac    >  1:2719414/1‑150 (MQ=255)
                                                                                                                                                 agagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAg  <  1:2126642/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  NC_000913/1351082‑1351376

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