Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 5.2% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC 5.2% 551.5 / 16.1 155intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (78/69);  new base C (4/4);  total (82/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 159 total aligned reads displayed.

TGGCGCAGGCCAAAGAAGCGGATAAAGGCGGCATGAATTTCTCGTTTATTAATTCGGCAGGCCAGTATCAACTCGAGGCTAAAAAATACGTTCGCCGCATCAGAGATAAAGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGG  >  NC_000913/4087819‑4088347
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                
tGGCGCAGGCCAAAGAAGCGGATAAAGGCGGCATGAATTTCTCGTTTATTAATTCGGCAGGCCAGTATCAACTCGAGGCTAAAAAATACGTTCGCCGCATCAGAGATAAAGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                                                                                                                                                                                 <  4:67525/274‑1 (MQ=255)
         ccAAAGAAGCGGATAAAGGCGGCATGAATTTCTCGTTTATTAATTCGGCAGGCAAGTATCAACTCGAGGCTAAAAAATACGTTCGCCGCATCAGAGATAAAGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                                                                                                                                                                                       <  1:4867/275‑1 (MQ=255)
         acAAAAAAGCGGAAAAAGGTGGCAAGATNTNCTCTTTTATTAATTCGGCAGGCCAGTATCAACTCGAGGCTAAAAANTACGTTCGCCGCATCAGAGATAAAGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACNCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                                                                                                                                                                                       <  3:37675/274‑1 (MQ=255)
                                              tcttaATTCGGCAGGCAAGAACAAAATGAGGGCTAAAAAAAACGTTCGCCGCATCAGAGGTAAAAGCCCGTACTCTGACTGGGATAAAGGGCAACTGCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCt                                                                                                                                                                                                                   <  4:176306/272‑1 (MQ=255)
                                                                           aGGCTAAAAAATACGTTCGCCGCATCCGAGATAAAGTCCCGTACTCTGACTGTGATAAAGAGCAACTTCAGGACGCAAACCCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCCAGTACAACGAAATGGTTGATGACTGTAATAGCCTGCGTTAAACGTTATTGATGATACATGAAGTGTCGGCGTCGGCGCATCCTACACATTTACCCCCTTGGAAATCAACGCCAAAAAAGCTCCATATCCTTTTATCCCGGCACTATtcccg                                                                                                                                                                                      >  4:31146/1‑271 (MQ=255)
                                                                                      tACGTTCGCCGCATCAGAGATAAAGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg                                                                                                                                                                           >  2:254769/1‑274 (MQ=255)
                                                                                                           aaaGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCNGTGTNGGATGCGGCGCATCCGACACAGTTCCCCCCTtgtg                                                                                                                                                                                                                                       >  4:356401/1‑193 (MQ=255)
                                                                                                               tCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                                                                                                                                                                                   <  2:356401/193‑1 (MQ=255)
                                                                                                                            cTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgg                                                                                                                                                                                                                                                                  >  5:730305/1‑149 (MQ=255)
                                                                                                                                gATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                                                                                                                                                                              >  6:1206746/1‑149 (MQ=255)
                                                                                                                                gATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTc                                                                                                                                                                                                                                                 >  4:443/1‑162 (MQ=255)
                                                                                                                                 aTAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAt                                                                                                                                                                                                                                                             >  5:647816/1‑149 (MQ=255)
                                                                                                                                  tAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                                                                                                                            >  8:186422/1‑149 (MQ=255)
                                                                                                                                  tAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                                                                                                                            >  5:1283848/1‑149 (MQ=255)
                                                                                                                                  tAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                                                                                                                                                                                            >  6:202937/1‑149 (MQ=255)
                                                                                                                                     agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGa                                                                                                                                                                                                                                                         <  8:719832/149‑1 (MQ=255)
                                                                                                                                          aaCTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                                                                                                                                                                                                    <  8:913496/149‑1 (MQ=255)
                                                                                                                                             ttCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaa                                                                                                                                                                                                                         <  2:220770/173‑1 (MQ=255)
                                                                                                                                               cAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                                                                                                                                                                               <  7:1180641/149‑1 (MQ=255)
                                                                                                                                                    cGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                                                                                                                                                                                          <  5:80677/149‑1 (MQ=255)
                                                                                                                                                    cGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                                                                                                                                                                                          <  5:206454/149‑1 (MQ=255)
                                                                                                                                                     gCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                                                                                                                                                                                         >  6:96089/1‑149 (MQ=255)
                                                                                                                                                       aaaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTGTAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACAAAATATGCTTGCCAGATAGCTAAATATTTTTAGCTAATGACCCGCATCGTCGGCGGCCCAGGc                                                                                                         >  4:193960/1‑275 (MQ=255)
                                                                                                                                                       aaaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATATCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgtg                                                                                                                                                                                                                                       >  7:992289/1‑149 (MQ=255)
                                                                                                                                                         aCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                                                                                                                                     <  8:710897/149‑1 (MQ=255)
                                                                                                                                                           tCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                                                                                                                                                                                       <  7:623030/129‑1 (MQ=255)
                                                                                                                                                               gCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACt                                                                                                                                                                                                                               >  7:1008800/1‑149 (MQ=255)
                                                                                                                                                                cTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                                                                                                                              <  5:1173415/149‑1 (MQ=255)
                                                                                                                                                                cTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                                                                                                                              <  6:1092986/149‑1 (MQ=255)
                                                                                                                                                                        tCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                                                                                                                                                                                      >  6:473730/1‑149 (MQ=255)
                                                                                                                                                                        tCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                                                                                                                                                                                      >  5:446533/1‑149 (MQ=255)
                                                                                                                                                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                                                                                                                                    >  6:197738/1‑149 (MQ=255)
                                                                                                                                                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                                                                                                                                    >  6:1049444/1‑149 (MQ=255)
                                                                                                                                                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                                                                                                                                    >  7:242840/1‑149 (MQ=255)
                                                                                                                                                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                                                                                                                                    >  5:734229/1‑149 (MQ=255)
                                                                                                                                                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                                                                                                                                    >  8:591497/1‑149 (MQ=255)
                                                                                                                                                                            aGACTCTTCCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGATTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                                                                                                                                                                                  >  6:1247422/1‑149 (MQ=255)
                                                                                                                                                                             gACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                                                                                                                                                                                 >  5:460748/1‑149 (MQ=255)
                                                                                                                                                                               ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                                                                                                                                                                                                <  7:855349/116‑1 (MQ=255)
                                                                                                                                                                                tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCatat                                                                                                                                                                                                              >  6:1334451/1‑149 (MQ=255)
                                                                                                                                                                                tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCatat                                                                                                                                                                                                              >  5:1038917/1‑149 (MQ=255)
                                                                                                                                                                                 ctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                                                                                                                                                                                             <  8:460759/149‑1 (MQ=255)
                                                                                                                                                                                  tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCt                                                                                                                                                                                                                   >  3:22902/1‑142 (MQ=255)
                                                                                                                                                                                  tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                                                                                                                                                                                            <  6:1251173/149‑1 (MQ=255)
                                                                                                                                                                                   ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                                                                                                                                           <  6:9060/149‑1 (MQ=255)
                                                                                                                                                                                   ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                                                                                                                                           <  7:485790/149‑1 (MQ=255)
                                                                                                                                                                                      ccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCata                                                                                                                                                                                                               <  1:22902/142‑1 (MQ=255)
                                                                                                                                                                                      ccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                                                                                                                                                                                        >  6:500556/1‑149 (MQ=255)
                                                                                                                                                                                                   gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                                                                                                                                                                                           >  7:1440574/1‑149 (MQ=255)
                                                                                                                                                                                                   gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                                                                                                                                                                                           >  8:1355618/1‑149 (MQ=255)
                                                                                                                                                                                                     gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                                                                                                                                                                                         <  6:730305/149‑1 (MQ=255)
                                                                                                                                                                                                        tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                                                                                                                                                                      >  7:490082/1‑149 (MQ=255)
                                                                                                                                                                                                        tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                                                                                                                                                                      <  8:265406/149‑1 (MQ=255)
                                                                                                                                                                                                           aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                                                                                                                                                                   >  6:330862/1‑149 (MQ=255)
                                                                                                                                                                                                              gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                                                                                                                                                                >  5:234817/1‑149 (MQ=255)
                                                                                                                                                                                                                 aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc                                                                                                                                                                             <  8:1334697/149‑1 (MQ=255)
                                                                                                                                                                                                                         tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACaa                                                                                                                                                                     <  6:647816/149‑1 (MQ=255)
                                                                                                                                                                                                                          gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                                                                                                                                                                    >  5:368597/1‑149 (MQ=255)
                                                                                                                                                                                                                          gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                                                                                                                                                                    <  6:470811/149‑1 (MQ=255)
                                                                                                                                                                                                                             tttaataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                                                                                                                                                                 >  6:274238/3‑149 (MQ=255)
                                                                                                                                                                                                                             tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                                                                                                                                                                 >  8:1394368/1‑149 (MQ=255)
                                                                                                                                                                                                                               taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt                                                                                                                                                               >  5:325759/1‑149 (MQ=255)
                                                                                                                                                                                                                                aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt                                                                                                                                                              <  6:368597/149‑1 (MQ=255)
                                                                                                                                                                                                                                    gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCa                                                                                                                                                          >  6:401154/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                                                                                                                                                                                             >  4:92355/1‑63 (MQ=255)
                                                                                                                                                                                                                                       tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGa                                                                                                                                                        >  5:1247664/1‑148 (MQ=255)
                                                                                                                                                                                                                                       tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt                                                                                                                                                       <  7:1394368/149‑1 (MQ=255)
                                                                                                                                                                                                                                        gCGTTAAACGTTTTGGATAATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTc                        <  2:288263/275‑1 (MQ=255)
                                                                                                                                                                                                                                         cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATTCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCCGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCTCAGGCCATCACACCATGGTCAATGACCTCGTTCTTCCCCTCAAATGGTGGGGCGGGGGCAATAAAAAATCCACGAAACCGTCGCTcc                       >  4:274072/1‑275 (MQ=255)
                                                                                                                                                                                                                                         cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACCCAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGCCATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAACCGTCGCTcc                       >  2:9722/1‑275 (MQ=255)
                                                                                                                                                                                                                                         cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAAGCGCATCGGCGGCGGCGCAGGCGATCACGCCATGGTTAATGACTTCGTTATTCACTTCAACTGTTCGGGCGAGGGTAATAAAAAATCAACGAAAAccttgatca                       >  1:48894/1‑267 (MQ=255)
                                                                                                                                                                                                                                          gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGa                                                                                                                                                        <  1:25325/145‑1 (MQ=255)
                                                                                                                                                                                                                                          gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc                                                                                                                                                    >  7:598487/1‑149 (MQ=255)
                                                                                                                                                                                                                                           ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                                                                                                                                                                                         <  2:92355/63‑1 (MQ=255)
                                                                                                                                                                                                                                             aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                                                                                                                                <  2:363643/150‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                                                                                                                                <  5:180629/149‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                                                                                                                                >  5:891793/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                                                                                                                                >  6:1416718/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                                                                                                                                                                                         >  7:376769/1‑91 (MQ=255)
                                                                                                                                                                                                                                                cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAata                                                                                                                                              <  5:1423603/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttt                                                                                                                                           <  6:128260/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAtttt                                                                                                                                          <  5:1180870/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAg           <  4:271770/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTGCACCCTTGTGAAGTAACTCCATAAATGCCTCATAGCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCATATAGCTAAATATGTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCa         <  2:276967/274‑1 (MQ=255)
                                                                                                                                                                                                                                                       gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAg                                                                                                                                       <  7:1355618/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                                                                                                                                   >  7:1331618/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                                                                                                                                   <  6:1209560/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                                                                                                                                   >  6:1430237/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                                                                                                                                   <  6:234817/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCATTAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCAt                                                                                                                                  >  8:625401/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGa                                                                                                                                <  7:432262/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAcgc                                                                                                                             <  8:502106/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCg                                                                                >  4:236619/1‑190 (MQ=255)
                                                                                                                                                                                                                                                                     tgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCt               >  4:232966/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                     tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                                                                                                                                         <  6:1008284/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg                                                                                                                       >  6:791207/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCAt                                                                                                                                  <  5:650210/136‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAg           <  2:232966/255‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcg                                                                                                                    >  7:939073/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCggcggcggc                                                                                                               >  8:827510/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                
TGGCGCAGGCCAAAGAAGCGGATAAAGGCGGCATGAATTTCTCGTTTATTAATTCGGCAGGCCAGTATCAACTCGAGGCTAAAAAATACGTTCGCCGCATCAGAGATAAAGTCCCGTACTCTGACTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCGCAGGCGATCAAGACATGGTTAATGACTTCGTTATTCACCTCAACTGGTTGGGCGAGGGTAATAAAAAATCCACGAAAACGTCGCTCCTGTTCGCTCCAGCAATGTGGG  >  NC_000913/4087819‑4088347

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