Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,388 T→C 5.8% intergenic (+397/+10) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,3880TC 5.8% 334.1 / 7.8 104intergenic (+397/+10)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/43);  new base C (2/4);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 108 total aligned reads displayed.

TTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGC  >  NC_000913/1552240‑1552536
                                                                                                                                                    |                                                                                                                                                    
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                                                                                                                                                      <  1:414422/149‑1 (MQ=255)
  ttttCATGCTGGATCAGGCGTTCACGCCGCATCCTGCATGAAAACTCTTCCAATTCAATAGAATGTAGTGAACGCGGAGGCCTGATCAGCGTAGCGCATCAGGCAATGTTGGGTTGGTCATCAGTTCTAAATAGCGCTTTATGAAGTGc                                                                                                                                                    >  4:296965/3‑149 (MQ=9)
       aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTTATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                                                                                                                                               >  1:662412/1‑149 (MQ=255)
       aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                                                                                                                                               <  2:221339/149‑1 (MQ=255)
         gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                                                                                                                                             >  1:374047/1‑149 (MQ=255)
           cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttagaa                                                                                                                                           <  4:799506/149‑5 (MQ=255)
            ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTACAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                                                                                                                                          >  4:728047/1‑149 (MQ=255)
              aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                        <  4:564569/149‑1 (MQ=255)
              aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                        <  1:779799/149‑1 (MQ=255)
               tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTg                                                                                                                                       >  2:281639/1‑149 (MQ=255)
                   gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa                                                                                                                                   >  1:327967/1‑149 (MQ=255)
                     gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                                                                 <  2:642496/149‑1 (MQ=255)
                         aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGa                                                                                                                             >  4:656984/1‑149 (MQ=255)
                         aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGa                                                                                                                             >  4:715187/1‑149 (MQ=255)
                         aCGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGa                                                                                                                             <  1:438370/149‑1 (MQ=255)
                            ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACgg                                                                                                                          >  2:169412/1‑149 (MQ=255)
                            ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACgg                                                                                                                          <  1:720041/149‑1 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                                                                 >  2:155566/1‑141 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                                         >  2:320308/1‑149 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                                         <  3:827628/149‑1 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttag                                                                                                                                             <  2:316081/129‑3 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttag                                                                                                                                             >  1:316081/1‑127 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                                                                 <  1:155566/141‑1 (MQ=255)
                              ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg                                                                                                                        >  3:578517/1‑149 (MQ=255)
                              ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttagaa                                                                                                                                           <  1:92981/130‑5 (MQ=255)
                              ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttagaa                                                                                                                                           >  2:92981/1‑126 (MQ=255)
                                 tCCGGCATGAAAACCCTTCAAATCCAATAGATTACAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTaa                                                                                                                     >  3:557963/1‑149 (MQ=255)
                                  ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                                                                         >  4:204583/1‑128 (MQ=255)
                                  ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                                                                         <  3:204583/128‑1 (MQ=255)
                                   cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATc                                                                                                                   <  1:848214/149‑1 (MQ=255)
                                   cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATc                                                                                                                   <  1:590196/149‑1 (MQ=255)
                                       aTGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                               <  2:788131/149‑1 (MQ=255)
                                           aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGAcaca                                                                                                           >  3:601233/1‑149 (MQ=255)
                                            aaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACAt                                                                                                          <  3:715187/149‑1 (MQ=255)
                                               ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCa                                                                                                       >  3:733275/1‑149 (MQ=255)
                                               ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCa                                                                                                       >  3:177025/1‑149 (MQ=255)
                                                 ttCAAATCCAATAGATTTTAGTGAAGGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTGTATAAAGTGCCATTTTTTTATTGTGTAACCAGACGGCATAATCGCGACACATCCACt                                                                                                     >  2:624084/1‑149 (MQ=21)
                                                 ttCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCCCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACt                                                                                                     >  2:271012/1‑149 (MQ=255)
                                                 ttCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATAAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACt                                                                                                     >  4:161652/1‑149 (MQ=255)
                                                    aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                                                                                                                                              >  1:496402/1‑105 (MQ=255)
                                                    aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                                                                                                                                              <  2:496402/105‑1 (MQ=255)
                                                    aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGtt                                                                                                  <  4:365936/149‑1 (MQ=255)
                                                         cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACg                                                                                                                           >  2:808584/1‑119 (MQ=255)
                                                         cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACg                                                                                                                           <  1:808584/119‑1 (MQ=255)
                                                         cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                             <  3:75593/149‑1 (MQ=255)
                                                           aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctg                                                                                           >  3:402504/1‑149 (MQ=255)
                                                           aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctg                                                                                           >  1:343385/1‑149 (MQ=255)
                                                            tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgc                                                                                          <  1:160826/149‑1 (MQ=255)
                                                                ttGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctct                                                                                      <  1:169412/149‑1 (MQ=255)
                                                                   cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctctgcc                                                                                   >  4:353787/1‑149 (MQ=255)
                                                                      tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGa                                                                                >  3:460628/1‑149 (MQ=255)
                                                                           gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCa                                                                           >  1:704078/1‑149 (MQ=255)
                                                                            tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAg                                                                          <  1:636419/149‑1 (MQ=255)
                                                                                   cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                            >  2:557117/1‑124 (MQ=255)
                                                                                   cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                            <  1:557117/124‑1 (MQ=255)
                                                                                   cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                   >  4:371205/1‑149 (MQ=255)
                                                                                   cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                   >  1:616262/1‑149 (MQ=255)
                                                                                    tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTg                                                                  >  4:575440/1‑149 (MQ=255)
                                                                                     gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGc                                                                 <  4:578517/149‑1 (MQ=255)
                                                                                      aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa                                                                <  1:355318/149‑1 (MQ=255)
                                                                                      aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa                                                                >  2:353528/1‑149 (MQ=255)
                                                                                      aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCa                                                                <  3:656984/149‑1 (MQ=255)
                                                                                       tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAg                                                               >  1:131695/1‑149 (MQ=255)
                                                                                       tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAg                                                               <  2:818168/149‑1 (MQ=255)
                                                                                        aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGa                                                              <  1:320308/149‑1 (MQ=255)
                                                                                              aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGc                                                        <  3:525830/149‑1 (MQ=255)
                                                                                              aGCGCACCAGGCAATGTTGCGTTTGTGATCAGTTCTAAATAGCGGTTTATAAAGTGCCATTTTTTTATTGCGTAACCGGACGCCGTAATCGCGACACCTCCACTGTTTTTTCTGCCTCTGCCAGTCTCCAGGGGTTTTGCGGATTCAGc                                                        <  2:98987/149‑1 (MQ=255)
                                                                                               gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc                                                       <  3:30736/149‑1 (MQ=255)
                                                                                                 gcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCcac                                                     <  2:689254/149‑1 (MQ=255)
                                                                                                   aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt                                                                                                                                         >  3:823893/1‑63 (MQ=255)
                                                                                                   aTCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTACCCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctctgc                                                                                    <  2:410937/116‑1 (MQ=255)
                                                                                                    tCAGGCAATGTTGCGATTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTACCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCcacact                                                  >  4:456193/1‑147 (MQ=255)
                                                                                                       ggCAATGTTGCGTTTGTTATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgc                                               <  2:662412/149‑1 (MQ=255)
                                                                                                       ggCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgc                                               >  1:560224/1‑149 (MQ=255)
                                                                                                        gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcg                                              >  2:898143/1‑149 (MQ=255)
                                                                                                              ttGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAAc                                        >  4:728886/1‑149 (MQ=255)
                                                                                                                 cGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGc                                     <  3:371205/149‑1 (MQ=255)
                                                                                                                 cGTTTGTCATCAGTTCTAAATGGCCCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGc                                     <  1:271012/149‑1 (MQ=255)
                                                                                                                       tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATca                               >  1:830167/1‑149 (MQ=255)
                                                                                                                         atcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcacc                             >  3:495554/1‑149 (MQ=255)
                                                                                                                         atcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcacc                             >  2:403718/1‑149 (MQ=255)
                                                                                                                          tcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcacca                            >  2:358641/1‑149 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcgg                                             >  3:753422/1‑128 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcgg                                             <  4:753422/128‑1 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGa                        >  4:238301/1‑149 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGa                        >  3:698366/1‑149 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGACAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGa                        >  4:632853/1‑149 (MQ=255)
                                                                                                                                   aaTGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGtt                   <  3:864103/149‑1 (MQ=255)
                                                                                                                                    aTGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGttt                  >  1:16374/1‑149 (MQ=255)
                                                                                                                                      ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTaa                <  1:403718/149‑1 (MQ=255)
                                                                                                                                       gcgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAAt               >  3:62752/1‑149 (MQ=255)
                                                                                                                                        cgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATc              >  3:107534/1‑149 (MQ=255)
                                                                                                                                        cgcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATc              <  2:616262/149‑1 (MQ=255)
                                                                                                                                           tttATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGcc           >  4:326598/1‑149 (MQ=255)
                                                                                                                                           tttATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGcc           >  1:888338/1‑149 (MQ=255)
                                                                                                                                            ttATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCa          <  3:575440/149‑1 (MQ=255)
                                                                                                                                                   gTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTgg   <  2:895552/149‑1 (MQ=255)
                                                                                                                                                    tGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGc  >  2:235357/1‑149 (MQ=255)
                                                                                                                                                    tGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGc  <  4:62752/149‑1 (MQ=255)
                                                                                                                                                    tGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGc  >  2:331845/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGC  >  NC_000913/1552240‑1552536

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