Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,228 G→C 8.4% 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,2280GC 8.4% 746.4 / 19.3 224intergenic (‑189/+179)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base G (74/128);  new base C (19/3);  total (93/131)
Fisher's exact test for biased strand distribution p-value = 1.07e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

TCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGC  >  NC_000913/1351080‑1351377
                                                                                                                                                    |                                                                                                                                                      
   ggCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                                                                                                                                    >  2:3818641/1‑150 (MQ=255)
     cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCCTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGc                                                                                                                                                  <  2:3619401/150‑1 (MQ=255)
      gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGGTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCCCCCGAACCCTCACCCCAGAGGGGCg                                                                                                                                                 >  2:1393092/1‑150 (MQ=255)
        ttCGTCGTTTTGGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                               >  2:4131622/1‑150 (MQ=255)
        ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATCGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                               <  1:127520/150‑1 (MQ=255)
        ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                               <  1:4114077/150‑1 (MQ=255)
        ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                               <  1:1687317/150‑1 (MQ=255)
         tcgtcgTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                                                                                                                              <  2:3029315/150‑1 (MQ=255)
          cgtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCGGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAggg                                                                                                                                             <  1:2426249/150‑1 (MQ=255)
            tcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGa                                                                                                                                           <  2:3706673/150‑1 (MQ=255)
               ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttagg                                                                                                                                        <  1:2791785/150‑6 (MQ=255)
               ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCGAACCCACACCCCACAGGGGCGAGGGGACCg                                                                                                                                        >  1:3659530/1‑150 (MQ=255)
               tttaAGTAAACAGTACGAACAGATAAACGGATAGTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCCGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCg                                                                                                                                        >  2:3129494/1‑150 (MQ=255)
                    gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCTGGTCCGACCGGAGCAGGCTGCGGCTGGGGGGGCGCTTTTGCGCTGACCCTAACCCTCTGCCCAGAGGGGCGTGGGGACCGTATtg                                                                                                                                   >  2:959457/1‑150 (MQ=255)
                      aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAACCCACTCCCCAGAGGGGCGAGGGGACCGTATtg‑tg                                                                                                                                >  1:662759/1‑150 (MQ=255)
                        aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCa                                                                                                                              <  2:1210966/150‑1 (MQ=255)
                         cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCaa                                                                                                                             <  1:1469047/150‑1 (MQ=255)
                          aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCaaa                                                                                                                            <  1:3511855/150‑1 (MQ=255)
                          aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCAACCAAAAACACAAACCACACGGGCGAGGGGACCGAATTG‑TGCaaa                                                                                                                            >  1:102122/1‑150 (MQ=255)
                            tACGAACAGATAAACGGTTATGATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAata                                                                                                                          <  1:156891/150‑1 (MQ=255)
                              cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATAtt                                                                                                                        <  2:2774989/150‑1 (MQ=255)
                              cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATAtt                                                                                                                        <  1:1587266/150‑1 (MQ=255)
                               gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCCCCCTAACCCTCAACCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTg                                                                                                                       >  2:3624441/1‑150 (MQ=255)
                               gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAACCCACACCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAAAATTg                                                                                                                       >  2:3290562/1‑150 (MQ=255)
                                 aCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCAAAACCACACCCCACCCGGGCGAGGGGACCGTATTG‑TGCAAATATTGtt                                                                                                                     >  1:1971368/1‑150 (MQ=255)
                                   aGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCTGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTCAc                                                                                                                   <  1:1344407/150‑1 (MQ=255)
                                    gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTAcc                                                                                                                  <  1:239188/150‑1 (MQ=255)
                                    gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCGAACCCTATCCCCCGAGGGGCGAGGGGTCCGTATTG‑TGCAAATATTTTTAcc                                                                                                                  >  2:2298340/1‑150 (MQ=255)
                                    gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTAcc                                                                                                                  >  2:1009125/1‑150 (MQ=255)
                                       aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCa                                                                                                               <  1:3943168/150‑1 (MQ=255)
                                       aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCTAAACCACACCCCAGAGGGGCGAGGGGCCCGTATTG‑TGCAAAAAATGTTACCCCa                                                                                                               >  2:1248210/1‑150 (MQ=255)
                                       aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCACCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATAATGTTAACCCa                                                                                                               >  2:2054658/1‑150 (MQ=255)
                                       aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCACACCCTAACCCACTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCCAATATTTTTACCCCa                                                                                                               >  2:1206615/1‑150 (MQ=39)
                                        aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAg                                                                                                              <  2:1462663/150‑1 (MQ=255)
                                        aaCGGTTATTAAAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAg                                                                                                              <  2:857588/150‑1 (MQ=255)
                                              tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGcaacaa                                                                                                        <  1:2937430/150‑1 (MQ=255)
                                               attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCaacaaa                                                                                                       <  1:4157104/150‑1 (MQ=255)
                                               attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCaacaaa                                                                                                       <  1:3536344/150‑1 (MQ=255)
                                                ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTGCCCACACCCAAACCCACTACCCACAGGGGCGGGGGGGCCGGATTGTTGCAAA‑AATGTTACCCCACCaacaaaa                                                                                                      >  1:1967196/1‑149 (MQ=255)
                                                 tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCaacaaaca                                                                                                     <  2:3681252/150‑1 (MQ=255)
                                                  ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAg                                                                                                    <  1:442504/150‑1 (MQ=255)
                                                   taATCAACCTGGCTGTGAGTAGCTATCGTTGCCAGGTCCGCCCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAgg                                                                                                   <  2:4092313/150‑1 (MQ=255)
                                                   taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAgg                                                                                                   <  2:2978066/150‑1 (MQ=255)
                                                     aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCACCAAACACAAAACCACAACCCAGAGGGGCGAGGGGACCCGATTG‑TGCAAAAAAAGTTAAACCAGCaaccaaaaagca                                                                                                 >  2:542359/1‑144 (MQ=255)
                                                      tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTGGCGGGGCCCTAACCCTATGCCCAGAGGGGGGTGGGGTGCGTTTTT‑TTGAAATATTTTTTAACCAACCACAAACAGGCTc                                                                                                >  1:461081/1‑150 (MQ=255)
                                                      tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCCCCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACACACAGGcgc                                                                                                >  2:3489802/1‑148 (MQ=255)
                                                      tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTc                                                                                                <  2:1734370/150‑1 (MQ=255)
                                                                tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGcccc                                                                                      <  2:3496247/150‑1 (MQ=255)
                                                                    aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCACCTAAc                                                                                  <  2:417132/150‑1 (MQ=255)
                                                                     gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAcc                                                                                 <  1:905954/150‑1 (MQ=255)
                                                                     gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAcc                                                                                 <  1:972968/150‑1 (MQ=255)
                                                                       aGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCt                                                                               <  2:2794222/150‑1 (MQ=255)
                                                                         cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCAAACCCttt                                                                             >  2:675575/1‑150 (MQ=255)
                                                                            taGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCAt                                                                          <  2:2186007/150‑1 (MQ=255)
                                                                              gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAAACCACTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTTTTAACCCAGCAACAAACAGGCTCATACAGCCCCAAACCCATTAAttg                                                                        >  2:1150671/1‑148 (MQ=38)
                                                                               ttGCCCGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTACCCTCAACCGAACCCTCTCCCCTGAGGGGCGGGGGGTCCGTATTG‑TGCAAAAATTGTTGCCCCCTCCACCAAAAGGCTCATTAAGACCCAAACCCAGTCATGGc                                                                       >  2:777723/1‑150 (MQ=255)
                                                                               ttGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc                                                                       <  1:4135911/150‑1 (MQ=255)
                                                                               ttGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc                                                                       <  2:1902037/150‑1 (MQ=255)
                                                                                 gCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGa                                                                     <  1:1645958/150‑1 (MQ=255)
                                                                                  ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTATCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCAAACCCTTTCATGGCGAt                                                                    >  2:2073956/1‑150 (MQ=255)
                                                                                  ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGAt                                                                    <  1:3298658/150‑1 (MQ=255)
                                                                                  ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGAt                                                                    <  2:4052596/150‑1 (MQ=255)
                                                                                  ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGAt                                                                    <  1:4077966/150‑1 (MQ=255)
                                                                                  ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCACTCCCCAGAGGGGCGAGGGGACCGTAGTG‑TGCCAATATTGTTACCCCCGCAACAAACACGCGCATACAGCCCCTAACCCGTTCATGGcccg                                                                    >  2:825956/1‑147 (MQ=38)
                                                                                            ccGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACg                                                          <  2:2015223/150‑1 (MQ=255)
                                                                                                     gCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCAAGGCGATGGCTGGGACGGTTCAGAcc                                                 >  2:595420/1‑150 (MQ=255)
                                                                                                     gCTGCGGCAGGGGGGGCGCTTTTCCCCCCACCCAAACCCGCTCCCCAGAGGGGCGAGGGGACCCTATTG‑TGCAAATATTGTTACCCCCGCAACAAACAGGCTCCTACAGCCCCGAACCCGTTCATGGCGCGGGCTGGGACGGTTCAGAcc                                                 >  2:3983548/1‑150 (MQ=38)
                                                                                                       tGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                               <  2:1196591/150‑1 (MQ=255)
                                                                                                       tGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                               <  2:2871699/150‑1 (MQ=255)
                                                                                                         cGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc                                             <  1:1865599/150‑1 (MQ=255)
                                                                                                           gCAGGGGGGGCGCTTTTCCCCTCACCCAAACCCTAACCCCAGAGGGGCGTGGGGTCCGTATTT‑TGCAAAAATTTTTAACACACCAACAAACAGGCTCAAAAAGCCCCAAAACCATTAATTGCGATGGCTGGGACGGTTAAGAACTTGCCg                                           >  2:425581/1‑150 (MQ=255)
                                                                                                              gggggggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAat                                        <  2:3939667/150‑1 (MQ=255)
                                                                                                                  gggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtc                                    <  1:231890/150‑1 (MQ=255)
                                                                                                                    gcgcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtctc                                  <  2:746160/150‑1 (MQ=255)
                                                                                                                      gcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCa                                <  1:2681678/150‑1 (MQ=255)
                                                                                                                       cTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAg                               <  1:3037679/150‑1 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                        >  1:3876064/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                        <  1:3627454/150‑1 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtct                       >  2:2026988/1‑150 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtct                       >  1:3966627/1‑150 (MQ=255)
                                                                                                                                tCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtctc                      <  1:2992134/150‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGttt               <  2:1338810/150‑1 (MQ=255)
                                                                                                                                       aaCCCTCGCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTGACCCCAGCAACAAACAGGCTCATACAGCCCCTACCCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGttt               <  2:1338644/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTc             >  2:1685234/1‑150 (MQ=255)
                                                                                                                                         cccTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTc             <  1:2676818/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATACTCCAGCACCGTCTCCATGTTTTc             >  2:2264059/1‑150 (MQ=255)
                                                                                                                                            tctcCCCAGAGGGGCGAGGGGACCGGATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcacc          <  1:734869/150‑1 (MQ=255)
                                                                                                                                                cccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa      >  1:631052/1‑150 (MQ=255)
                                                                                                                                                 ccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCactac     >  1:3202701/1‑150 (MQ=255)
                                                                                                                                                 ccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac     >  2:24634/1‑150 (MQ=255)
                                                                                                                                                 ccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac     >  2:4039354/1‑150 (MQ=255)
                                                                                                                                                 ccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac     >  2:1574730/1‑150 (MQ=255)
                                                                                                                                                 ccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac     >  2:2353475/1‑150 (MQ=255)
                                                                                                                                                 ccAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac     >  2:1779568/1‑150 (MQ=255)
                                                                                                                                                 acAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGGCTCCATCTTTTCACCacaac     >  1:3734633/2‑150 (MQ=255)
                                                                                                                                                  cAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca    <  2:461081/150‑1 (MQ=255)
                                                                                                                                                  cAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca    <  2:3255948/150‑1 (MQ=255)
                                                                                                                                                  cAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca    >  2:2392552/1‑150 (MQ=255)
                                                                                                                                                    gagGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGc  >  2:1062645/1‑150 (MQ=255)
                                                                                                                                                    gagGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGc  >  2:3100136/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTG‑TGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGC  >  NC_000913/1351080‑1351377

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: ‑