Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,332 A→T 5.9% A212A (GCT→GCA hdfR ← flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,3320AT 5.9% 412.4 / 12.5 118A212A (GCT→GCAhdfRflhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (57/54);  new base T (4/3);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 126 total aligned reads displayed.

TCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAA  >  NC_000913/3947186‑3947479
                                                                                                                                                  |                                                                                                                                                   
tCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcagc                                                                                                                                                   >  1:542938/1‑149 (MQ=255)
 cGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcagc                                                                                                                                                   >  3:550040/1‑148 (MQ=255)
   cTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGTTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGa                                                                                                                                                >  2:131943/1‑149 (MQ=255)
    tATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                                                                               <  1:119487/149‑1 (MQ=255)
     aTTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                                                                               >  4:199504/1‑148 (MQ=255)
     aTTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGaa                                                                                                                                               >  3:894927/1‑148 (MQ=255)
        ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                           <  2:447161/149‑1 (MQ=255)
        ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                           <  2:12837/149‑1 (MQ=255)
        ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTg                                                                                                                                            <  4:550040/148‑1 (MQ=255)
         tGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                           <  3:199504/148‑1 (MQ=255)
         tGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                           <  4:379470/148‑1 (MQ=255)
         tGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                          <  2:992391/149‑1 (MQ=255)
         tGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAATGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                          <  3:662281/148‑1 (MQ=255)
          gCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggca                                                                                                                                         <  2:1041998/149‑3 (MQ=255)
          gCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGt                                                                                                                                          <  3:898182/148‑1 (MQ=255)
                       aTAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                                                                     <  2:174648/124‑1 (MQ=255)
                       aTAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTca                                                                                                                                                     >  1:174648/1‑124 (MQ=255)
                       aTAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCa                                                                                                                            >  2:371346/1‑149 (MQ=255)
                        tAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGtt                                                                                                                                         <  2:720708/135‑1 (MQ=255)
                        tAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGtt                                                                                                                                         >  1:720708/1‑135 (MQ=255)
                             cGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGt                                                                                                                      >  2:509364/1‑149 (MQ=255)
                              ggCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGt                                                                                                                      <  3:815346/148‑1 (MQ=255)
                                     aaaGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCg                                                                                                              >  1:701235/1‑149 (MQ=255)
                                      aaGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGa                                                                                                             >  2:513992/1‑149 (MQ=255)
                                         tgttgtGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCa                                                                                                          >  2:855514/1‑149 (MQ=255)
                                           ttgtGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCaa                                                                                                         <  4:499122/148‑1 (MQ=255)
                                              tgCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAAcc                                                                                                      <  4:597723/148‑1 (MQ=255)
                                                cTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTg                                                                                                    <  3:763797/148‑1 (MQ=255)
                                                  aTCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGcc                                                                                                  >  4:508638/1‑148 (MQ=255)
                                                  aTCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGca                                                                                                  >  3:25724/1‑147 (MQ=255)
                                                   tCGACAACGGTATGCAGGCCGCCTTTTTTACGCTCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCt                                                                                                 <  4:698857/148‑1 (MQ=255)
                                                   tCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCt                                                                                                 >  3:321543/1‑148 (MQ=255)
                                                     gACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCa                                                                                                                                     <  2:954586/110‑1 (MQ=255)
                                                     gACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCAt                                                                                              >  2:373042/1‑149 (MQ=255)
                                                     gACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCACAt                                                                                              <  1:131943/149‑1 (MQ=255)
                                                        aaCGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCt                                                                                           >  1:633698/1‑149 (MQ=255)
                                                          cGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGt                                                                                         <  1:704456/149‑1 (MQ=255)
                                                           ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGt                                                                                         <  4:368832/148‑1 (MQ=255)
                                                             tATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGa                                                                                      >  1:211037/1‑149 (MQ=255)
                                                                gCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGaaa                                                                                    <  3:773072/148‑1 (MQ=255)
                                                                 cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGaaaa                                                                                   <  4:894927/148‑1 (MQ=255)
                                                                  aGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                                                                 <  2:888538/149‑1 (MQ=255)
                                                                  aGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                                                                 >  1:888538/1‑149 (MQ=255)
                                                                   ggCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                                                                 >  4:42706/1‑148 (MQ=255)
                                                                   ggCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCt                                                                                <  2:83454/149‑1 (MQ=255)
                                                                     ccgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTg                                                                               >  3:254996/1‑148 (MQ=255)
                                                                       gccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGTcc                                                                            >  1:223575/1‑149 (MQ=255)
                                                                        ctttttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccc                                                                           >  2:869436/3‑149 (MQ=255)
                                                                         cTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGcccc                                                                          >  2:1046023/1‑149 (MQ=255)
                                                                          ttttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                         >  1:322738/1‑149 (MQ=255)
                                                                           tttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                         >  3:325897/1‑148 (MQ=255)
                                                                           tttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                         >  4:256041/1‑148 (MQ=255)
                                                                           tttttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccccc                                                                         >  4:218914/1‑148 (MQ=255)
                                                                                aCGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCaa                                                                   <  2:633698/149‑1 (MQ=255)
                                                                                  gcgcCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGt                                                                 <  1:855514/149‑1 (MQ=255)
                                                                                     cccAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCg                                                               <  3:218914/148‑1 (MQ=255)
                                                                                     cccAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCg                                                               <  4:814553/148‑1 (MQ=255)
                                                                                     cccAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCg                                                               >  4:159801/1‑148 (MQ=255)
                                                                                      ccAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCa                                                             <  1:509364/149‑1 (MQ=255)
                                                                                       cAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCa                                                             >  3:542355/1‑148 (MQ=255)
                                                                                       cAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCa                                                             >  3:859288/1‑148 (MQ=255)
                                                                                           tGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATaa                                                        <  1:315021/149‑1 (MQ=255)
                                                                                               gggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                    <  1:391137/149‑1 (MQ=255)
                                                                                                 gTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAg                                                   <  3:42706/148‑1 (MQ=255)
                                                                                                  tAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGAt                                                 <  1:513992/149‑1 (MQ=255)
                                                                                                      cAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGa                                                                                                             >  4:276766/1‑85 (MQ=255)
                                                                                                      cAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGa                                                                                                             <  3:276766/85‑1 (MQ=255)
                                                                                                      cAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGTCCCCCACTCAAGTCGCAGATAATTAAGATCTcc                                             <  2:223575/149‑1 (MQ=255)
                                                                                                            cAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTa                                        >  4:87379/1‑148 (MQ=255)
                                                                                                             aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGtt                                                                                        <  1:889772/99‑1 (MQ=255)
                                                                                                              aCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGt                                      >  4:130156/1‑148 (MQ=255)
                                                                                                                cATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGttt                                    <  4:414102/148‑1 (MQ=255)
                                                                                                                   ttaGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTATCAGAATGGGCACTTGGTCAGCCCGTTTCACAGCAGCCTCATGCTGTTGGAAATCTTGGCCCCCCTCAATTTGCAAATAATTAAGATCTTCCATAAGTTTTGa                                 >  4:625462/3‑148 (MQ=255)
                                                                                                                   tAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGa                                 >  4:134905/1‑148 (MQ=255)
                                                                                                                   tAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCCATCAAACGTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGaa                                <  2:1023033/149‑1 (MQ=255)
                                                                                                                      gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAggg                             <  1:41716/149‑1 (MQ=255)
                                                                                                                      gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAggg                             <  1:1059946/149‑1 (MQ=255)
                                                                                                                      gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAgg                              <  4:246147/148‑1 (MQ=255)
                                                                                                                      gCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAgg                              <  3:130156/148‑1 (MQ=255)
                                                                                                                         tCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGc                           >  3:922000/1‑148 (MQ=255)
                                                                                                                          cGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCa                          >  4:351981/1‑148 (MQ=255)
                                                                                                                            cAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTg                       <  1:1072689/149‑1 (MQ=255)
                                                                                                                            cAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTg                       >  1:374543/1‑149 (MQ=255)
                                                                                                                             aaTCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTgg                      <  2:161488/149‑1 (MQ=255)
                                                                                                                              aTCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTgg                      <  3:256041/148‑1 (MQ=255)
                                                                                                                               tCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGta                    >  1:28165/1‑149 (MQ=255)
                                                                                                                               tCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGta                    >  1:785104/1‑149 (MQ=255)
                                                                                                                                    tGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATaaa                >  4:1028309/1‑148 (MQ=255)
                                                                                                                                      cTGTGCCAGTTCAGCAGAGCTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAgc              >  3:622205/1‑148 (MQ=255)
                                                                                                                                       tgtgCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATaaa                <  2:745619/145‑1 (MQ=255)
                                                                                                                                         aaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCGGCCACATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTa          >  2:752125/3‑149 (MQ=255)
                                                                                                                                          accAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTa          >  3:437003/2‑148 (MQ=255)
                                                                                                                                           ccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTaa         <  4:437003/148‑1 (MQ=255)
                                                                                                                                            cAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTaaa        >  3:882794/1‑148 (MQ=255)
                                                                                                                                               ttCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGa                                                  >  3:1006869/1‑103 (MQ=255)
                                                                                                                                               ttCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGa                                                  <  4:1006869/103‑1 (MQ=255)
                                                                                                                                               ttCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTg     >  3:156209/1‑148 (MQ=255)
                                                                                                                                                 cagcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGaaa  <  2:322738/149‑1 (MQ=255)
                                                                                                                                                 cagcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGaaa  <  2:559651/149‑1 (MQ=255)
                                                                                                                                                 cagcagAACTGGTTGTCAGAAGGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGaa   <  3:508638/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAA  >  NC_000913/3947186‑3947479

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