Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,335 A→T 6.7% S211S (TCT→TCA hdfR ← flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,3350AT 6.7% 370.1 / 13.5 105S211S (TCT→TCAhdfRflhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (48/50);  new base T (4/3);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 110 total aligned reads displayed.

CGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATA  >  NC_000913/3947187‑3947481
                                                                                                                                                    |                                                                                                                                                  
 gCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcagcag                                                                                                                                                   >  4:508905/1‑149 (MQ=255)
   tATTTTGCAGCCATATGGCATAAAGCGGTCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                                                                                 <  4:580449/149‑1 (MQ=255)
   tATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGa                                                                                                                                                  >  1:583043/1‑148 (MQ=255)
    aTTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                                                                                 >  2:432288/1‑148 (MQ=255)
    aTTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                                                                                 >  1:336430/1‑148 (MQ=255)
       ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTgg                                                                                                                                             <  3:845891/149‑1 (MQ=255)
       ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTg                                                                                                                                              <  2:785331/148‑1 (MQ=255)
       ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                             <  3:508905/149‑1 (MQ=255)
       ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTg                                                                                                                                              <  1:633616/148‑1 (MQ=255)
        tGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                            <  3:867306/149‑1 (MQ=255)
          cAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGtt                                                                                                                                           <  1:119607/148‑1 (MQ=255)
             ccgtatGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGCAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTc                                                                                                                                        <  2:392699/145‑1 (MQ=255)
                  tGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAAt                                                                                                                                   <  2:540623/148‑1 (MQ=255)
                  tGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATg                                                                                                                                  <  3:855432/149‑1 (MQ=255)
                    gCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATggg                                                                                                                                <  4:388845/149‑1 (MQ=255)
                    gCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATggg                                                                                                                                <  4:200596/149‑1 (MQ=255)
                    gCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATgg                                                                                                                                 <  1:409121/148‑1 (MQ=255)
                       tAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCa                                                                                                                              >  1:904360/1‑148 (MQ=255)
                       tAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCAc                                                                                                                             >  3:423184/1‑149 (MQ=255)
                         aaGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACt                                                                                                                            <  2:904360/148‑1 (MQ=255)
                            cGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCg                                                                                                                         <  2:443501/148‑1 (MQ=255)
                                 gTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCa                                                                                                                   <  3:34562/149‑1 (MQ=255)
                                  tGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCa                                                                                                                   >  1:470801/1‑148 (MQ=255)
                                   gAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCAcc                                                                                                                 >  3:229864/1‑149 (MQ=255)
                                    aaaGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                            <  1:194086/121‑1 (MQ=255)
                                    aaaGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                            >  2:194086/1‑121 (MQ=255)
                                     aaGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCg                                                                                                                >  2:685657/1‑148 (MQ=255)
                                      aGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGa                                                                                                               >  1:631742/1‑148 (MQ=255)
                                      aGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGAt                                                                                                              <  4:334287/149‑1 (MQ=255)
                                        tgttgtGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCa                                                                                                            >  3:657418/1‑149 (MQ=255)
                                          ttgtGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCaa                                                                                                           <  1:674787/148‑1 (MQ=255)
                                          ttgtGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCaa                                                                                                           >  2:237228/1‑148 (MQ=255)
                                           tgtgCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCaaa                                                                                                          <  2:336430/148‑1 (MQ=255)
                                               cTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTg                                                                                                      <  1:544991/148‑1 (MQ=255)
                                                 aTCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGt                                                                                                                        <  2:262033/128‑1 (MQ=255)
                                                 aTCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGt                                                                                                                        >  1:262033/1‑128 (MQ=255)
                                                   cGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTc                                                                                                  >  1:79520/1‑148 (MQ=255)
                                                     aCAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATg                                                                                               >  4:679235/1‑149 (MQ=255)
                                                     aCAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATg                                                                                               >  3:234990/1‑149 (MQ=255)
                                                      cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                            >  2:506217/1‑103 (MQ=255)
                                                      cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                            <  1:506217/103‑1 (MQ=255)
                                                      cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATg                                                                                               <  2:600697/148‑1 (MQ=255)
                                                      cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATg                                                                                               >  2:819893/1‑148 (MQ=255)
                                                      cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATg                                                                                               >  2:513735/1‑148 (MQ=255)
                                                      cAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTCTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCTTg                                                                                               >  2:46768/1‑148 (MQ=255)
                                                         cGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTg                                                                                            <  1:685657/148‑1 (MQ=255)
                                                         cGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGt                                                                                           <  3:99349/149‑1 (MQ=255)
                                                          ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGtt                                                                                          <  4:234990/149‑1 (MQ=255)
                                                          ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGt                                                                                           <  1:232617/148‑1 (MQ=255)
                                                          ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGt                                                                                           <  1:513735/148‑1 (MQ=255)
                                                                cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGaaaa                                                                                     >  3:282939/1‑148 (MQ=255)
                                                                cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGaaaa                                                                                     <  4:282939/148‑1 (MQ=255)
                                                                cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAAt                                                                                    <  4:423184/149‑1 (MQ=255)
                                                                   gccgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCt                                                                                  <  1:237228/148‑1 (MQ=255)
                                                                    ccgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTgg                                                                                <  4:657418/149‑1 (MQ=255)
                                                                     cgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGc                                                                               <  4:818418/149‑1 (MQ=255)
                                                                      gccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGc                                                                               >  1:253633/1‑148 (MQ=255)
                                                                        cTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTgg                                                                                <  4:547687/145‑1 (MQ=255)
                                                                          tttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                           >  2:340989/1‑148 (MQ=255)
                                                                          tttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                           >  2:15120/1‑148 (MQ=255)
                                                                          tttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCAATTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                           <  2:583043/148‑1 (MQ=255)
                                                                              tACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTc                                                                       >  2:564982/1‑148 (MQ=255)
                                                                               aCGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCa                                                                      <  1:340989/148‑1 (MQ=255)
                                                                                          tGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATa                                                           <  2:470801/148‑1 (MQ=255)
                                                                                                       ggTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTccc                                              >  1:327548/1‑148 (MQ=255)
                                                                                                         tGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCtt                                            >  1:235494/1‑148 (MQ=255)
                                                                                                           cAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTa                                          >  2:800367/1‑148 (MQ=255)
                                                                                                            aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGt                                        >  3:362529/1‑149 (MQ=255)
                                                                                                             aCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGtt                                       <  4:229864/149‑1 (MQ=255)
                                                                                                               cATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGtttt                                     <  3:503501/149‑1 (MQ=255)
                                                                                                               cATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGttt                                      <  2:235494/148‑1 (MQ=255)
                                                                                                                aTTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTg                                    >  3:398080/1‑149 (MQ=255)
                                                                                                                aTTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTg                                    >  4:563100/1‑149 (MQ=255)
                                                                                                                 ttAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTg                                    >  1:138733/1‑148 (MQ=255)
                                                                                                                 ttAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTg                                    >  1:166831/1‑148 (MQ=255)
                                                                                                                  tAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGa                                   >  1:766882/1‑148 (MQ=255)
                                                                                                                    aGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTc                                                                  <  3:824350/115‑1 (MQ=255)
                                                                                                                    aGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTc                                                                  >  4:824350/1‑115 (MQ=255)
                                                                                                                    aGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAgg                                >  3:808491/1‑149 (MQ=255)
                                                                                                                    aGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAg                                 <  2:79520/148‑1 (MQ=255)
                                                                                                                     gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                                                                   <  1:187807/97‑1 (MQ=255)
                                                                                                                     gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                                                                   >  2:187807/1‑97 (MQ=255)
                                                                                                                          gCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCAc                           >  2:386292/1‑148 (MQ=255)
                                                                                                                           cAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTg                         >  4:629940/1‑149 (MQ=255)
                                                                                                                                  ttGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                      >  4:735002/1‑113 (MQ=255)
                                                                                                                                  ttGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                      <  3:735002/113‑1 (MQ=255)
                                                                                                                                  ttGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGttt                                      <  4:574718/129‑1 (MQ=255)
                                                                                                                                  ttGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGttt                                      >  3:574718/1‑129 (MQ=255)
                                                                                                                                  ttGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATaa                   >  2:345713/1‑148 (MQ=255)
                                                                                                                                     cTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAgc                >  1:873138/1‑148 (MQ=255)
                                                                                                                                      acaaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAgcgc              >  3:889462/5‑149 (MQ=255)
                                                                                                                                        aaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAgc                <  4:461434/143‑1 (MQ=255)
                                                                                                                                        aaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAgc                >  3:461434/3‑145 (MQ=255)
                                                                                                                                         accAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTa            >  2:299182/2‑148 (MQ=255)
                                                                                                                                          ccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCCCCGATCACACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCCGATAATTACGATCTCCCTTTAGTTTTGACGGGGCACTGGTCTAACGCGCtca           >  2:40366/1‑146 (MQ=255)
                                                                                                                                           cAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTaaa          >  2:330090/1‑148 (MQ=255)
                                                                                                                                           cAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTaaa          >  1:329378/1‑148 (MQ=255)
                                                                                                                                                  gcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAata  >  3:108667/1‑149 (MQ=255)
                                                                                                                                                  gcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAata  >  4:397149/1‑149 (MQ=255)
                                                                                                                                                  gcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTACCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAat   <  2:296470/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATA  >  NC_000913/3947187‑3947481

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