Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,388 T→C 5.7% 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.7% 457.1 / 12.7 139intergenic (+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 (72/59);  new base C (5/3);  total (77/62)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 143 total aligned reads displayed.

TTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGG  >  NC_000913/1552240‑1552535
                                                                                                                                                    |                                                                                                                                                   
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                                                                                                                                                     >  1:1409325/1‑149 (MQ=255)
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                                                                                                                                                     <  3:1411114/149‑1 (MQ=255)
tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                                                                                                                                                     >  1:1118719/1‑149 (MQ=255)
      cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                                                                                                                                               <  3:808018/149‑1 (MQ=255)
      cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtt                                                                                                                                               <  4:1179321/149‑1 (MQ=255)
        tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                                                                                                                                             >  3:546611/1‑149 (MQ=255)
        tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                                                                                                                                             <  2:1204700/149‑1 (MQ=255)
           cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCa                                                                                                                                                 >  2:496853/1‑142 (MQ=255)
           cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt                                                                                                                                          <  2:289922/149‑1 (MQ=255)
            ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                                                                                                                                         >  1:1071532/1‑149 (MQ=255)
            ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                                                                                                                                         <  2:615670/149‑1 (MQ=255)
            ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                                                                                                                                         <  1:768358/149‑1 (MQ=255)
              ataaggcgTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                       >  1:1300497/5‑149 (MQ=17)
              aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                       <  4:140386/149‑1 (MQ=255)
              aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                       <  3:1235116/149‑1 (MQ=255)
                     gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc                                                                                                                                >  1:39029/1‑149 (MQ=255)
                      tGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAcc                                                                                                                               >  4:388033/1‑149 (MQ=255)
                           gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACg                                                                                                                          >  1:298303/1‑149 (MQ=255)
                             cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                                        <  3:1456969/149‑1 (MQ=255)
                              ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg                                                                                                                       >  4:415513/1‑149 (MQ=255)
                               tATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGt                                                                                                                      <  1:405657/149‑1 (MQ=255)
                                 tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTaa                                                                                                                    <  2:732400/149‑1 (MQ=255)
                                   cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATc                                                                                                                  >  2:1474542/1‑149 (MQ=255)
                                     gCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCa                                                                                                                              >  2:410763/1‑135 (MQ=255)
                                     gCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcgc                                                                                                                <  2:1071532/149‑1 (MQ=255)
                                     gCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATcgc                                                                                                                >  2:1327868/1‑149 (MQ=255)
                                          aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGacac                                                                                                           >  3:1128470/1‑149 (MQ=255)
                                            aaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACAt                                                                                                         >  2:503448/1‑149 (MQ=255)
                                               ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCa                                                                                                      >  4:1019695/1‑149 (MQ=255)
                                                cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                              <  2:902994/140‑1 (MQ=255)
                                                cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                     >  3:1329103/1‑149 (MQ=255)
                                                cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                     >  4:819367/1‑149 (MQ=255)
                                                cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                     >  3:244359/1‑149 (MQ=255)
                                                cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCAc                                                                                                     <  3:388033/149‑1 (MQ=255)
                                                    aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGtt                                                                                                 >  4:198016/1‑149 (MQ=255)
                                                         cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                            >  2:1097436/1‑149 (MQ=255)
                                                         cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtc                                                                                            >  4:106416/1‑149 (MQ=255)
                                                          aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                           >  2:848072/1‑149 (MQ=255)
                                                          aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtct                                                                                           >  2:543026/1‑149 (MQ=255)
                                                            tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgc                                                                                         <  2:672422/149‑1 (MQ=255)
                                                            tAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgc                                                                                         <  3:298245/149‑1 (MQ=255)
                                                                 tGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                              >  3:264292/1‑123 (MQ=255)
                                                                 tGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                              <  4:264292/123‑1 (MQ=255)
                                                                    aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCa                                                                                 >  2:694293/1‑149 (MQ=255)
                                                                      tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGa                                                                               <  1:681561/149‑1 (MQ=255)
                                                                       gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGAc                                                                              >  1:1007295/1‑149 (MQ=255)
                                                                            tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAg                                                                         <  1:107059/149‑1 (MQ=255)
                                                                            tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAg                                                                         <  1:483909/149‑1 (MQ=255)
                                                                             gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAgg                                                                        >  4:290535/1‑149 (MQ=255)
                                                                               aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTCTAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCg                                                                      >  3:1229420/1‑149 (MQ=255)
                                                                                 gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                    <  2:1039912/149‑1 (MQ=255)
                                                                                 gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                    >  2:470181/1‑149 (MQ=255)
                                                                                 gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                    >  2:1478966/1‑149 (MQ=255)
                                                                                 gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                    <  1:853695/149‑1 (MQ=255)
                                                                                   cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                  <  2:1097334/149‑1 (MQ=255)
                                                                                     gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGc                                                                <  4:801696/149‑1 (MQ=255)
                                                                                        aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGa                                                             <  2:386280/149‑1 (MQ=255)
                                                                                        aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGa                                                             <  1:627873/149‑1 (MQ=255)
                                                                                          gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAtt                                                                                                                                       <  1:508063/73‑1 (MQ=255)
                                                                                          gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGAtt                                                           >  3:1172140/1‑149 (MQ=255)
                                                                                            gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCa                                                         >  4:587512/1‑149 (MQ=255)
                                                                                              aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGc                                                       <  3:309639/149‑1 (MQ=255)
                                                                                               gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc                                                      <  3:977988/149‑1 (MQ=255)
                                                                                               gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc                                                      >  3:320199/1‑149 (MQ=255)
                                                                                                 gcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCcac                                                    >  2:1016921/1‑149 (MQ=255)
                                                                                                    tCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAtt                                                 >  1:1192742/1‑149 (MQ=255)
                                                                                                       ggCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgc                                              <  3:686839/149‑1 (MQ=255)
                                                                                                        gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCa                                                         >  2:1404206/1‑137 (MQ=255)
                                                                                                        gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCa                                                         <  1:1404206/137‑1 (MQ=255)
                                                                                                        gCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcg                                             >  3:822801/1‑149 (MQ=255)
                                                                                                          aaTGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCGGGCGTTTTGCAGATTCAGCCACATTTgcggc                                           <  1:848072/149‑1 (MQ=255)
                                                                                                          aaTGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcggc                                           <  3:819367/149‑1 (MQ=255)
                                                                                                                 cGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGc                                    <  3:106416/149‑1 (MQ=255)
                                                                                                                   tttGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCg                                  <  2:45679/149‑1 (MQ=255)
                                                                                                                   tttGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCg                                  <  4:1065172/149‑1 (MQ=255)
                                                                                                                       tcatcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATca                              >  4:486372/1‑149 (MQ=255)
                                                                                                                          tcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                              >  1:557876/1‑66 (MQ=255)
                                                                                                                          tcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGa                                                                                                              <  2:557876/66‑1 (MQ=255)
                                                                                                                          tcaGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcacca                           >  1:175789/1‑149 (MQ=255)
                                                                                                                           aataaaaaaaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATcaccac                          >  1:1393274/8‑149 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                  <  2:6203/106‑1 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                  >  1:6203/1‑106 (MQ=255)
                                                                                                                              ttCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGa                       <  1:479314/149‑1 (MQ=255)
                                                                                                                               aaaaaaTGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaa                      >  2:1176082/4‑149 (MQ=255)
                                                                                                                                cTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTGTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGaag                     <  2:1254190/149‑2 (MQ=255)
                                                                                                                                 tAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAg                    >  2:1401532/1‑149 (MQ=255)
                                                                                                                                  aaaTGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt                   >  1:453499/1‑149 (MQ=255)
                                                                                                                                  aaaTGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGt                   >  2:8589/1‑149 (MQ=255)
                                                                                                                                   aaTGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTAGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGtt                  >  4:480355/1‑149 (MQ=255)
                                                                                                                                     tGGCACTTTATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTa                >  2:614018/1‑149 (MQ=255)
                                                                                                                                      ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTaa               <  4:546229/149‑1 (MQ=255)
                                                                                                                                      ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTaa               >  3:410943/1‑149 (MQ=255)
                                                                                                                                      ggCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTaa               <  4:1329103/149‑1 (MQ=255)
                                                                                                                                         gcTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCg            >  4:40672/1‑149 (MQ=255)
                                                                                                                                           tttATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGcc          >  1:1154807/1‑149 (MQ=255)
                                                                                                                                           tttATAAAGCGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGAc                                                                              >  3:1061948/1‑81 (MQ=255)
                                                                                                                                               taAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                  <  1:757098/89‑1 (MQ=255)
                                                                                                                                               taAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                  >  2:757098/1‑89 (MQ=255)
                                                                                                                                                aaaGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTc     >  2:1183055/1‑149 (MQ=255)
                                                                                                                                                   gTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTgg  <  4:1322477/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGG  >  NC_000913/1552240‑1552535

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