Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 248,263 G→A 5.1% intergenic (+129/+95) rayT → / ← lfhA RAYT REP element‑mobilizing transposase; TnpA(REP)/pseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,2630GA 5.1% 248.3 / 2.2 80intergenic (+129/+95)rayT/lfhARAYT REP element‑mobilizing transposase; TnpA(REP)/pseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base A (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

GAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCG  >  NC_000913/248117‑248406
                                                                                                                                                  |                                                                                                                                               
gAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGgc                                                                                                                                               <  3:240808/149‑1 (MQ=255)
         aTTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                      <  4:334975/149‑1 (MQ=255)
         aTTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                      <  8:131096/149‑1 (MQ=255)
         aTTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCAGCATCCGGCAGGTATGCGCAGATGCCTGATGCGACTCTGGCGCGTCTTa                                                                                                                                       >  4:186802/1‑148 (MQ=255)
         aTTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGGTTATGCCGGATGGGGAAGTTATGCGGAGATGCCTGATGCGACGGGTGCGCGtggtct                                                                                                                                      >  4:328495/1‑144 (MQ=21)
          tttCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGACGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                      >  8:335561/1‑148 (MQ=255)
             cATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTGTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                                                                                                   <  2:342926/148‑1 (MQ=255)
              aTAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                                                                                                 >  5:168930/1‑149 (MQ=255)
              aTAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                                                                                                 <  4:287634/149‑1 (MQ=255)
                  ttGTTCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGTATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAc                                                                                                                             <  3:186802/149‑1 (MQ=255)
                     tgCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaaa                                                                                                                          <  2:57177/149‑1 (MQ=255)
                          aGATGCCTGACGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTtgt                                                                                                                     <  7:94051/149‑1 (MQ=255)
                                   aTGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACg                                                                                                            >  7:404029/1‑149 (MQ=255)
                                    tGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACgg                                                                                                           <  8:404029/149‑1 (MQ=255)
                                      cGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                          <  7:335561/148‑1 (MQ=255)
                                       gACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTGATCATGCCTACAAACTTGTGCCTGAACGGTAg                                                                                                        >  3:225526/1‑149 (MQ=255)
                                               cgcgTCTTATCATGCCTACAAACTTGTGCCGTATCGGTAGGCCGGATAAGGCGTTAATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATa                                                                                                <  3:328495/149‑1 (MQ=255)
                                                  gTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAgg                                                                                             >  5:41556/1‑149 (MQ=255)
                                                   tCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTAATGCCGCATCCGGCACTTATGCGCAGCTGCCTGATGCGACGCTGGAGCGTCTTATCATTACTACAAACTGGTGCCTGAACGGTAGGACGGAGAAGGc                                                                                            >  4:234290/1‑149 (MQ=21)
                                                     ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaa                                                                                                                           <  2:291232/116‑1 (MQ=255)
                                                     ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaa                                                                                                                           >  1:291232/1‑116 (MQ=255)
                                                     ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGt                                                                                          >  5:173556/1‑149 (MQ=255)
                                                       aTCCTGCCTACCCACTTGCGCCGGATCTGTAGGCCGTATATGGCGTTTATGCCGCATCCGTCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGt                                                                                          <  6:80877/147‑1 (MQ=21)
                                                       aTCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTAAGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTAATCATGCCTACAAACTTGTGGCTGAACGGTAGGACGGATAAGGCGtt                                                                                         >  7:332914/1‑148 (MQ=255)
                                                        tCATGCCTACAACCTTGTGCCGGATCGGTTGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGAAGATTCCTGATCCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCa                                                                                       <  3:186200/149‑1 (MQ=37)
                                                         cATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTAATGCATCATCCGGCAGTTATGCTCAGATGCCTGATGAAACGCTGGCGCGAATTATCATGACGACAAACTTGTGACTGAACGGTAGGACGGATAAGGGGTTCAc                                                                                      >  4:186200/1‑149 (MQ=37)
                                                                 cAAACTTGTGCCGTATCGTGAGGCCGGATAAGGCGTTTATTCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCa                                                                                <  6:168930/147‑1 (MQ=255)
                                                                 cAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCa                                                                                <  5:161155/147‑1 (MQ=255)
                                                                 cAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATc                                                                              <  3:234290/149‑1 (MQ=255)
                                                                         tgCGGGAGCGGTAGGTCGGAGATGGCGTTTCTGCCGCATCCGGAATTTATGCGCAGATGCCTGATGCGACGGTGGCGCGTCTTATCATGGCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGtt                                                                      <  1:107749/149‑1 (MQ=18)
                                                                          gCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTtg                                                                     >  7:389199/1‑149 (MQ=255)
                                                                           ccGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCAGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTtgt                                                                    <  2:142031/149‑1 (MQ=255)
                                                                            cGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTtgtg                                                                   <  2:3091/149‑1 (MQ=255)
                                                                               aTCGGTAGGCCGGATCAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCAc                                                                <  5:353770/149‑1 (MQ=255)
                                                                                   gTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGCCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGAGCACCAAt                                                            >  6:372169/1‑149 (MQ=255)
                                                                                     aGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTCGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGc                                                          >  8:49394/1‑149 (MQ=255)
                                                                                         cGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGa                                                      >  7:85110/1‑149 (MQ=255)
                                                                                          ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                                                                           <  4:277873/127‑1 (MQ=255)
                                                                                          ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                                                                           >  3:277873/1‑127 (MQ=255)
                                                                                           gATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATa                                                    <  8:85110/149‑1 (MQ=255)
                                                                                             tAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACt                                                  >  7:133823/1‑149 (MQ=255)
                                                                                                ggCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTa                                                >  8:296289/1‑148 (MQ=255)
                                                                                                 gCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACa                                              >  2:421559/1‑149 (MQ=255)
                                                                                                 gCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTGATCATGCCGACAAACTTGTGCCTGAACGGTAGGAAGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACa                                              >  2:292862/1‑149 (MQ=255)
                                                                                                   gTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                                                       <  1:145961/58‑1 (MQ=255)
                                                                                                   gTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                                                       >  2:145961/1‑58 (MQ=255)
                                                                                                    tttATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATc                                            >  8:115/1‑148 (MQ=255)
                                                                                                      tATGCCGCATACGGCAGTTATGCTCAGATGCCTTATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCaa                                          <  6:41556/148‑1 (MQ=255)
                                                                                                       aTGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACg                                        >  5:177988/1‑149 (MQ=255)
                                                                                                           cGCATCCTTCAGTTATGCGCAGATGCCTTATGCGACGCTTGCGCGTCTTATCATGCCTACAAACTTGTGCCTGATCGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTTGATACTTACATCAACGcccc                                    <  8:332914/149‑1 (MQ=21)
                                                                                                                 cGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTAGCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATaa                               <  1:311727/148‑1 (MQ=255)
                                                                                                                   gCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCt                            >  5:70553/1‑149 (MQ=255)
                                                                                                                   gCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCt                            >  5:177750/1‑149 (MQ=255)
                                                                                                                       ttgtGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                          <  3:304403/64‑1 (MQ=0)
                                                                                                                       ttgtGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                          <  3:309848/64‑1 (MQ=0)
                                                                                                                       ttgtGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                          >  4:309848/4‑67 (MQ=0)
                                                                                                                       ttgtGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                          >  4:304403/4‑67 (MQ=0)
                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGa                                                      >  7:267861/1‑119 (MQ=255)
                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGa                                                      <  8:267861/119‑1 (MQ=255)
                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCaa                        >  4:89505/1‑149 (MQ=255)
                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCaa                        <  5:372169/149‑1 (MQ=255)
                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCaa                        <  2:408210/149‑1 (MQ=255)
                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCaa                        >  1:198474/1‑149 (MQ=255)
                                                                                                                           gcgcAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTc                     >  6:143978/1‑148 (MQ=255)
                                                                                                                            cgcAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCaa                                                             <  7:40483/107‑1 (MQ=255)
                                                                                                                            cgcAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCaa                                                             >  8:40483/1‑107 (MQ=255)
                                                                                                                            cgcAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCa                   <  8:133823/149‑1 (MQ=255)
                                                                                                                             gcAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAg                  >  2:200041/1‑149 (MQ=255)
                                                                                                                                 aTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGCACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTc               >  1:294294/1‑148 (MQ=255)
                                                                                                                                    ccTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATc            <  7:115/148‑1 (MQ=255)
                                                                                                                                    ccTGATGCGACGCTGGCGCGTCGTATCATGCCTACAAACGTGTGCCTGCCCGGTAGGACGGAGAAGGCGTTCACGCCGCCGCCGGCAGTTGTGCCCCAATGCTGGATACTTACCTCAACGCCCCCATAATCTTCACCTCCCGCTCATcc           >  5:38587/1‑149 (MQ=37)
                                                                                                                                     cTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGTTCATCCg          >  3:405919/1‑149 (MQ=255)
                                                                                                                                     cTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCg          >  3:221712/1‑149 (MQ=255)
                                                                                                                                        aTGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGc        <  5:81877/148‑1 (MQ=255)
                                                                                                                                        aTGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCa       >  1:67664/1‑149 (MQ=255)
                                                                                                                                           cGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGccc                                     >  1:173864/1‑116 (MQ=255)
                                                                                                                                           cGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGccc                                     <  2:173864/116‑1 (MQ=255)
                                                                                                                                            gACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTc   >  4:224139/1‑149 (MQ=255)
                                                                                                                                             aCGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCg  >  8:3809/1‑149 (MQ=255)
                                                                                                                                             aCGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCg  >  4:77241/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCG  >  NC_000913/248117‑248406

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑