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

Only 100 of 236 total aligned reads displayed.

ATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGC  >  NC_000913/1351079‑1351377
                                                                                                                                                      |                                                                                                                                                     
aTCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGGTATTCTAATCAACCTGGCTGTGAGTAGCGATAGTTGCCAGGGCCGACCGGAGCAGGCTGCGGCCGGGGGGGCGCGGTTCCCCTCGCCCTAACCCGCGCCCCag                                                                                                                                                       >  2:4293220/1‑150 (MQ=255)
 tCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                                                                                                                                      <  1:1682249/150‑1 (MQ=255)
   aGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAACCCACACCCCAGCgg                                                                                                                                                    >  2:45018/1‑150 (MQ=255)
          tcgtcgTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                                                                                                                             <  1:243176/150‑1 (MQ=255)
          tcgtcgTTTTAGTAAACAGTACGAACAGAGAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  2:2908667/150‑1 (MQ=255)
           cgtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  1:2591882/150‑1 (MQ=255)
           cgtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  1:1617637/150‑1 (MQ=255)
            gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCACACCCGAACCCACTCCCCAGAGGGGCGAgggg                                                                                                                                           >  2:537657/1‑150 (MQ=255)
               gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAcc                                                                                                                                        <  2:1461168/150‑1 (MQ=255)
                 tttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGt                                                                                                                                      >  2:2829819/1‑150 (MQ=255)
                 tttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGt                                                                                                                                      >  2:1606054/1‑150 (MQ=255)
                  ttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                                                                                     <  1:487518/150‑1 (MQ=255)
                   tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAAACCTCGCCCCAGAGGGGCGAGGGGACCGGAt                                                                                                                                    >  1:2952873/1‑150 (MQ=255)
                    aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                                                                   <  2:4013051/150‑1 (MQ=255)
                    aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                                                                   <  2:1660593/150‑1 (MQ=255)
                      tAAACAGTACGAACAGATACACGGTTATGATAATCCACCTGGCTGTGCGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgt                                                                                                                                 <  2:1058118/150‑1 (MQ=255)
                      tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCGAACCCACACCCCAGAGGGGCGAGGGGACCCTATtgt                                                                                                                                 >  1:2805653/1‑150 (MQ=255)
                       aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgtg                                                                                                                                <  2:3379028/150‑1 (MQ=255)
                       aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgtg                                                                                                                                <  2:1383805/150‑1 (MQ=255)
                          cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCaa                                                                                                                             <  1:2081087/150‑1 (MQ=255)
                           aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCaaa                                                                                                                            <  1:1359114/150‑1 (MQ=255)
                             tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCCGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCGTTTCCCCGCACCCCAACCCGCTCCCCCCAGGGGCGCGGGGCCCGGAGTGGGCCAata                                                                                                                          >  2:4737186/1‑150 (MQ=255)
                              aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCGAACCCACACCCCCGAGGGGCGAGGGGACCGTATTGTGCAAatat                                                                                                                         >  1:4203720/1‑150 (MQ=255)
                               cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                                                                        <  1:1028577/150‑1 (MQ=255)
                                gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTg                                                                                                                       <  2:3904005/150‑1 (MQ=255)
                                gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTg                                                                                                                       >  2:3190009/1‑150 (MQ=255)
                                gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGCCTTGTGCAAATATTg                                                                                                                       >  2:2656194/1‑150 (MQ=255)
                                gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAAACCACACCCCACAGGGGCGAGGGGGCCGGATTGTGCAAATATTg                                                                                                                       >  1:308879/1‑150 (MQ=255)
                                gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCTCACCCAAACCCACTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTg                                                                                                                       >  1:3066923/1‑150 (MQ=255)
                                  aCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCAACCAAAACCACACCCCACAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt                                                                                                                     >  1:3109098/1‑150 (MQ=255)
                                    aGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAACCCACACCCCAGAGGGGCGAGGGGGCCCGATTGTGCCAAAATTTTTAc                                                                                                                   >  1:3350573/1‑150 (MQ=255)
                                    aGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTGTTCCCCCCACCCAAAACCACAACCCACAGGGGCGGGGGGGCCCGATTGTGCAAATATTGTTAc                                                                                                                   >  1:1424992/1‑150 (MQ=255)
                                     gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCAAACCCACTCCCCAGAGGGGCGAGGGGACCGGATTGTGCAAAAATTGTTAc                                                                                                                   >  2:2304050/1‑150 (MQ=39)
                                     gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCCCACCCAAACCCACACCCCACAGGGGCCACGGGCCCCGAATGTGCAAAAATTGTTaac                                                                                                                  >  2:4143591/1‑148 (MQ=255)
                                       tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcccc                                                                                                                <  1:3579591/150‑1 (MQ=255)
                                       tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcccc                                                                                                                <  1:1176625/150‑1 (MQ=255)
                                        aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATGGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCACACCCAAAACCACTCCCCAGAGGGGCGCGGGGACCGTATTGTGCCAATATTTTTACCCCa                                                                                                               >  1:3375483/1‑150 (MQ=255)
                                        aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa                                                                                                               >  2:2766065/1‑150 (MQ=255)
                                        aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa                                                                                                               >  2:724466/1‑150 (MQ=255)
                                        aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCGCACCCGAACCCACACCCCCGAGGGGCGAGGGGGCCGGATTGGGCAAATATTGTTGCCCCa                                                                                                               >  1:138153/1‑150 (MQ=255)
                                         aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAg                                                                                                              >  2:3567862/1‑150 (MQ=255)
                                         aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAg                                                                                                              <  1:2353625/150‑1 (MQ=255)
                                         aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCACCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAg                                                                                                              >  2:4227494/1‑150 (MQ=255)
                                               aattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGGCCGTATTGTGCAAATATTGTTACCCCAGcaacaa                                                                                                        >  1:3810112/2‑150 (MQ=255)
                                                attatACTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                                                                                       <  2:1729020/150‑1 (MQ=255)
                                                attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCAAACCCAAACCCACACCCCAGAGGGGCGAGGGGTCCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                                                                                       >  1:2313703/1‑150 (MQ=255)
                                                 ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaaca                                                                                                     <  1:182968/150‑1 (MQ=255)
                                                    taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                                                                                   <  2:2747576/150‑1 (MQ=255)
                                                    taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCAGACCGGTTCAGGCTTCGGCAGGGTGGGAGCTTTTCGCCTAACCCTAACCCTCTCCCAAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTTCCCCAGCaacaaactgt                                                                                                   <  2:2063165/150‑4 (MQ=255)
                                                       tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTc                                                                                                <  2:1369065/150‑1 (MQ=255)
                                                       tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGCGCTTTTCCCCTCACCCTAACCCGCTCCCCAGGGGGGCGAGGGGACCCGATTGTGCCAATATTGTTTCCCCCGCAAGAAACCgggtga                                                                                               >  2:463719/1‑146 (MQ=37)
                                                            cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTGCCCACCCCCGAACCCACACCCCACAGGGGGGCGGGGCCCCGATTTTGCAAAAATTTTTACCCCACCAACAAAAAGCCTCaaaaa                                                                                           >  1:78407/1‑146 (MQ=255)
                                                                 tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTGCCCTCACCCTAAACCTCTACCCCGAGGGGGGGGGGGACCGTTTTTTGGTAATTTTTTTTACCCAGCAACAAAAAGGGTCATACAGcccc                                                                                      >  1:3383961/1‑150 (MQ=255)
                                                                 tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCAAACCCAAACCCACAACCCACACGGGCCACGGGACCCAATTGTGCAAAAAATTTTAACCCACCAAAAAAAAAGCACAAACAAcccc                                                                                      >  1:1181347/1‑150 (MQ=255)
                                                                 tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCGAACCCAAAACCAAAACCCACAGGGGCCACGGGACCGTAATGTTCAAAAAATTTTAACCCAGCAACAAAAAGGCTCATACAAcccc                                                                                      >  2:161711/1‑150 (MQ=255)
                                                                  gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCt                                                                                     <  2:2059403/150‑1 (MQ=255)
                                                                    gAGTAGCTGTAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTACCCGCACCCGAACCCGCACCCCACGGGGGCGGGGGGGCCGGATTTTGCAAAAAATGTTTCCCCAGCAACAAAAAGGCGAAAAAAGCCCCaaa                                                                                   >  2:900120/1‑150 (MQ=255)
                                                                    gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTaa                                                                                   <  1:2463004/150‑1 (MQ=255)
                                                                    gAGTAGCTATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTaa                                                                                   <  1:4737186/150‑1 (MQ=255)
                                                                     aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTGCCCACACCCAAAACCACACCCCACACGGGCGACGGGCCCCTATTGTGCAAAAAATTTTTCCCCAGCAACAAAAACGCGCATAAAGCCCCAAAc                                                                                  >  1:3202222/1‑150 (MQ=255)
                                                                      gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAcc                                                                                 <  1:3025782/150‑1 (MQ=255)
                                                                      gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAAcc                                                                                 <  1:3585852/150‑1 (MQ=255)
                                                                              aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                         <  2:1853541/150‑1 (MQ=255)
                                                                               ggggCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTGTCCCCTCACCCTAACCCTCGCCCCCGAGGGGCGAGGGGACCGAATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATgg                                                                        <  2:3319699/147‑1 (MQ=255)
                                                                                  gCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCGCACCCTAAACCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTTTTACCCCAGCAACAAACAGGCTCATACAGCCCCGAACCCGTTCAGGGCgg                                                                     >  2:78551/1‑149 (MQ=38)
                                                                                     aGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATgg                                                                  <  2:4314430/150‑1 (MQ=255)
                                                                                       gTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTATCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCt                                                                <  1:2673317/150‑1 (MQ=255)
                                                                                         ccGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTgg                                                              <  2:2229746/150‑1 (MQ=255)
                                                                                            aCCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGCCCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGAc                                                           >  2:4345115/1‑150 (MQ=255)
                                                                                            aCCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGGTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCGTTCATGGCCATGGCTGGGag                                                           >  2:3457286/1‑149 (MQ=255)
                                                                                                gAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGtt                                                       <  1:2303582/150‑1 (MQ=255)
                                                                                                gAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGtt                                                       >  2:2418598/1‑150 (MQ=255)
                                                                                                  gCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCa                                                     <  2:867758/150‑1 (MQ=255)
                                                                                                     ggCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGAc                                                  <  2:2633703/150‑1 (MQ=255)
                                                                                                         gCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTg                                              <  1:2194510/150‑1 (MQ=255)
                                                                                                                        gcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCa                                <  2:3091684/150‑1 (MQ=255)
                                                                                                                          ttttCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCCAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGc                              >  2:251723/1‑150 (MQ=255)
                                                                                                                             tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCAcc                           <  2:3085999/150‑1 (MQ=255)
                                                                                                                                 cTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtct                       >  2:2799388/1‑150 (MQ=255)
                                                                                                                                 cTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtct                       >  1:39777/1‑150 (MQ=255)
                                                                                                                                     cccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCAt                   >  1:1798469/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCCTGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATg                  >  2:358730/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATg                  >  1:1665777/1‑150 (MQ=255)
                                                                                                                                       cTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGt                 >  2:2363335/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtt                <  2:1571002/150‑1 (MQ=255)
                                                                                                                                          aCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtttt              <  1:2766077/150‑1 (MQ=255)
                                                                                                                                          aCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtttt              <  1:724465/150‑1 (MQ=255)
                                                                                                                                          aCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtttt              <  1:724954/150‑1 (MQ=255)
                                                                                                                                            ccTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAACCAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCGGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca            >  2:1881381/1‑150 (MQ=255)
                                                                                                                                             ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac           <  1:2485860/150‑1 (MQ=255)
                                                                                                                                               ctcCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcacca         >  1:3524884/1‑150 (MQ=255)
                                                                                                                                                tcCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcaccac        >  1:4537873/1‑150 (MQ=255)
                                                                                                                                                tcCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCAAGTTTTcaccac        >  1:4383909/1‑150 (MQ=255)
                                                                                                                                                  cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa      >  2:3829704/1‑150 (MQ=255)
                                                                                                                                                  cccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa      >  1:2586099/1‑150 (MQ=255)
                                                                                                                                                   ccAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac     >  1:426364/1‑150 (MQ=255)
                                                                                                                                                     agagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAg   <  2:1194147/150‑1 (MQ=255)
                                                                                                                                                      gagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGc  >  2:2990838/1‑150 (MQ=255)
                                                                                                                                                      gagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGc  <  1:1864412/150‑1 (MQ=255)
                                                                                                                                                      gagGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGc  >  2:2593388/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAGC  >  NC_000913/1351079‑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: