Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,332 A→T 5.6% 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.6% 443.3 / 13.4 125A212A (GCT→GCAhdfRflhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (54/64);  new base T (4/3);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

Only 100 of 133 total aligned reads displayed.

ATCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAAT  >  NC_000913/3947185‑3947480
                                                                                                                                                   |                                                                                                                                                    
aTCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcag                                                                                                                                                     <  2:566318/149‑1 (MQ=255)
aTCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcag                                                                                                                                                     <  4:602289/149‑1 (MQ=255)
 tCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcagc                                                                                                                                                    <  3:416731/149‑1 (MQ=255)
     tATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcagca                                                                                                                                                   <  3:83380/146‑1 (MQ=255)
     tATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTcagca                                                                                                                                                   >  4:83380/1‑146 (MQ=255)
         ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                            <  2:397480/149‑1 (MQ=255)
         ttGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTgg                                                                                                                                            <  3:373549/149‑1 (MQ=255)
           gCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGtt                                                                                                                                          >  2:168401/1‑149 (MQ=255)
           gCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGtt                                                                                                                                          <  3:545543/149‑1 (MQ=255)
             aGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGt                                                                                                                                        >  1:1143018/1‑149 (MQ=255)
              gCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTg                                                                                                                                             <  3:303907/143‑1 (MQ=255)
                     ggCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATgg                                                                                                                                <  1:151316/149‑1 (MQ=255)
                         tAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCAc                                                                                                                            >  3:594364/1‑149 (MQ=255)
                         tAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCAc                                                                                                                            <  3:453008/149‑1 (MQ=255)
                          aaaGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACt                                                                                                                           <  3:247132/149‑1 (MQ=255)
                          aaaGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACt                                                                                                                           >  4:851204/1‑149 (MQ=255)
                            aGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTc                                                                                                                         <  1:151883/149‑1 (MQ=255)
                            aGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTc                                                                                                                         >  1:754536/1‑149 (MQ=255)
                              cGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGt                                                                                                                       <  4:35533/149‑1 (MQ=255)
                              cGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGt                                                                                                                       <  4:577365/149‑1 (MQ=255)
                               ggCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTc                                                                                                                      <  3:152650/149‑1 (MQ=255)
                                gCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCa                                                                                                                     <  2:754536/149‑1 (MQ=255)
                                 ccGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAg                                                                                                                    <  2:945756/149‑1 (MQ=255)
                                   gTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCa                                                                                                                  <  1:515755/149‑1 (MQ=255)
                                       aaGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGa                                                                                                              >  3:764837/1‑149 (MQ=255)
                                            ttgtGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCaaa                                                                                                         <  4:614511/149‑1 (MQ=255)
                                                 cTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGc                                                                                                    <  3:487082/149‑1 (MQ=255)
                                                  tATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGcc                                                                                                   >  2:266374/1‑149 (MQ=255)
                                                   aTCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCt                                                                                                  <  4:1095569/149‑1 (MQ=255)
                                                      gACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCAt                                                                                               <  3:508415/149‑1 (MQ=255)
                                                       aCAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTgg                                                                                                                                            <  1:327271/103‑1 (MQ=255)
                                                         aaCGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCt                                                                                            >  2:543429/1‑149 (MQ=255)
                                                            ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGaa                                                                                                                                   >  1:123690/1‑107 (MQ=255)
                                                            ggTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGtt                                                                                         >  4:1061891/1‑149 (MQ=255)
                                                              tATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGa                                                                                       <  3:329577/149‑1 (MQ=255)
                                                              tATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCAATAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGa                                                                                       >  2:212824/1‑149 (MQ=255)
                                                                 gCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGaaaa                                                                                    <  1:280614/149‑1 (MQ=255)
                                                                  cAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAAt                                                                                   >  3:195845/1‑149 (MQ=255)
                                                                   aGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATc                                                                                  >  1:619480/1‑149 (MQ=255)
                                                                     gccgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGAt                                                                                                             <  2:841631/120‑1 (MQ=255)
                                                                     gccgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGAt                                                                                                             >  1:841631/1‑120 (MQ=255)
                                                                       cgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggcac                                                                                                                                         <  2:443731/90‑4 (MQ=255)
                                                                       cgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCTGCTGAACTggcac                                                                                                                                         >  1:443731/1‑87 (MQ=255)
                                                                       cgccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGc                                                                              >  1:118029/1‑149 (MQ=255)
                                                                         ccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccc                                                                            >  4:387061/1‑149 (MQ=255)
                                                                         ccTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGccc                                                                            >  3:856854/1‑149 (MQ=255)
                                                                               ttACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTc                                                                      <  4:9551/149‑1 (MQ=255)
                                                                                 aCGCGCCCAGCTGACGGGTAGCCAGGTGCAACCACTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCaa                                                                    <  3:938001/149‑1 (MQ=255)
                                                                                    cgcCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTc                                                                 <  3:851204/149‑1 (MQ=255)
                                                                                    cgcCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTc                                                                 <  1:486104/149‑1 (MQ=255)
                                                                                      cccAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGc                                                               <  3:892381/149‑1 (MQ=255)
                                                                                            tGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGCTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATaa                                                         >  3:878468/1‑149 (MQ=255)
                                                                                            tGACGGGTAGCCAGGTGCAACCAATAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATaa                                                         <  1:212824/149‑1 (MQ=255)
                                                                                             gACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAAt                                                        >  4:27982/1‑149 (MQ=255)
                                                                                              aCGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAAtt                                                       <  2:118029/149‑1 (MQ=255)
                                                                                              aCGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAAtt                                                       >  1:306431/1‑149 (MQ=255)
                                                                                                gggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                     >  2:364234/1‑149 (MQ=255)
                                                                                                gggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                     <  4:881267/149‑1 (MQ=255)
                                                                                                gggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                     >  2:1098683/1‑149 (MQ=255)
                                                                                                gggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                     <  2:1029533/149‑1 (MQ=255)
                                                                                                gggTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                     >  3:70551/1‑149 (MQ=255)
                                                                                                   tAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGAt                                                  <  2:1143018/149‑1 (MQ=255)
                                                                                                    aGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGAtc                                                 >  1:940302/1‑149 (MQ=255)
                                                                                                      ccAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGAtctc                                               >  4:37524/1‑149 (MQ=255)
                                                                                                       cAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTcc                                              <  1:252166/149‑1 (MQ=255)
                                                                                                          gTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCtt                                           <  4:478425/149‑1 (MQ=255)
                                                                                                            gCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTa                                         <  4:168385/149‑1 (MQ=255)
                                                                                                              aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATa                                                          >  1:660501/1‑130 (MQ=255)
                                                                                                              aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATa                                                          <  2:660501/130‑1 (MQ=255)
                                                                                                              aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGt                                       >  1:1088410/1‑149 (MQ=255)
                                                                                                              aaCCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGt                                       >  3:762209/1‑149 (MQ=255)
                                                                                                                ccATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGttt                                     <  2:619480/149‑1 (MQ=255)
                                                                                                                 cATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTcc                                              <  3:797875/139‑1 (MQ=255)
                                                                                                                 cATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGtttt                                    <  1:762780/149‑1 (MQ=255)
                                                                                                                 cATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGACAATTAAGATCTcc                                              >  4:797875/1‑139 (MQ=255)
                                                                                                                  aTTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTg                                   >  3:427789/1‑149 (MQ=255)
                                                                                                                     aaGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAg                                <  3:675585/149‑1 (MQ=255)
                                                                                                                       acATCGCAATCTGGTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCAGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAggg                              >  4:991607/2‑149 (MQ=255)
                                                                                                                          tCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGc                                                               <  2:734311/113‑1 (MQ=255)
                                                                                                                          tCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGc                                                               >  1:734311/1‑113 (MQ=255)
                                                                                                                          tCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCa                           >  3:795897/1‑149 (MQ=255)
                                                                                                                           cGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCAc                          <  4:70551/149‑1 (MQ=255)
                                                                                                                              aaTCTGTTGCTGTGCCAGTTCAGCAGCACAGGTTGTGAGAATGGGCACTTCGTCAGCACTGATAAAACCCGCCTGATGCTCTTGAAAATCTGGCCCCCACTCAAGTCGCAGATGATTATCATCTCCCTTTAGTTTTGAAGGGGCACTgg                       <  4:710726/149‑1 (MQ=255)
                                                                                                                                 cTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGtat                    <  1:543429/149‑1 (MQ=255)
                                                                                                                                   gTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATaa                  >  1:255828/1‑149 (MQ=255)
                                                                                                                                    ttGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATaaa                 <  4:764837/149‑1 (MQ=255)
                                                                                                                                     tGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAg                >  3:494157/1‑149 (MQ=255)
                                                                                                                                      gCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAga               >  4:636065/1‑148 (MQ=255)
                                                                                                                                        acaaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAgcgc             >  2:305119/5‑149 (MQ=255)
                                                                                                                                          aaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAggg                              <  4:706827/128‑1 (MQ=255)
                                                                                                                                          aaccAGTTCTGCTGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGACCTCCCTTTAGTTTTGAAggg                              >  3:706827/3‑130 (MQ=255)
                                                                                                                                            ccAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTaa                                                     >  1:1062913/1‑105 (MQ=255)
                                                                                                                                            ccAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTcc                                              <  4:849104/112‑1 (MQ=255)
                                                                                                                                            ccAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTaaa         >  3:799241/1‑149 (MQ=255)
                                                                                                                                            ccAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTaaa         <  4:799241/149‑1 (MQ=255)
                                                                                                                                               gTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCTCTAAAGTg      >  2:935193/1‑149 (MQ=255)
                                                                                                                                                 tCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGaa    <  2:12884/149‑1 (MQ=255)
                                                                                                                                                   agcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGcc                                                                             <  4:43140/74‑1 (MQ=255)
                                                                                                                                                   agcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGcc                                                                             >  3:43140/1‑74 (MQ=255)
                                                                                                                                                   agcagAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAat  >  2:27341/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCGCTATTTTGCAGCCATATGGCATAAAGCGGCCGTGAAAGTGTTGTGCTATCGACAACGGTATGCAGGCCGCCTTTTTTACGCGCCCAGCTGACGGGTAGCCAGGTGCAACCATTAAGCATCGCAATCTGTTGCTGTGCCAGTTCAGCAGAACTGGTTGTCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAAT  >  NC_000913/3947185‑3947480

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: