Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 248,868 C→G 5.2% pseudogene (1203/1713 nt) lfhA ← 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,8680CG 5.2% 535.0 / 15.8 153pseudogene (1203/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (65/80);  new base G (5/3);  total (70/83)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 158 total aligned reads displayed.

GGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCT  >  NC_000913/248720‑249016
                                                                                                                                                    |                                                                                                                                                    
ggCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGcgcc                                                                                                                                                      >  4:188878/1‑149 (MQ=255)
        tCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCaa                                                                                                                                              >  3:91172/1‑149 (MQ=255)
          ggCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAAt                                                                                                                                            >  2:1338263/1‑149 (MQ=255)
             cgccAGCAGACTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCACCTGGCTGTAATTGAGCGCCGCGCTCAAATCtt                                                                                                                                         <  4:1102523/149‑1 (MQ=255)
             cgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCTGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaattaca                                                                                                                                         <  4:153844/149‑6 (MQ=255)
             cgccAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACAAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCGGCGCTCaattaca                                                                                                                                         <  1:133800/149‑6 (MQ=255)
                    agAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAg                                                                                                                                  >  3:300875/1‑149 (MQ=255)
                    agAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAg                                                                                                                                  >  1:1173219/1‑149 (MQ=255)
                      aaTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAc                                                                                                                                <  2:1283272/149‑1 (MQ=255)
                      aaTATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAc                                                                                                                                >  4:1140168/1‑149 (MQ=255)
                        tatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                                                                                                                              >  1:287151/1‑149 (MQ=255)
                        tatGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgc                                                                                                                              >  4:270236/1‑149 (MQ=255)
                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgcg                                                                                                                             <  4:91172/149‑1 (MQ=255)
                         atGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGAcgcg                                                                                                                             <  3:1350395/149‑1 (MQ=255)
                          tGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCgg                                                                                                                            >  4:202797/1‑149 (MQ=255)
                           gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                           <  1:415936/149‑1 (MQ=255)
                           gATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGt                                                                                                                           >  2:1009209/1‑149 (MQ=255)
                             tCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGc                                                                                                                         <  3:82650/149‑1 (MQ=255)
                               tttCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCg                                                                                                                       >  2:532156/1‑149 (MQ=255)
                                ttCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGt                                                                                                                      <  1:551302/149‑1 (MQ=255)
                                      aCCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTcgacga                                                                                                                <  4:192600/149‑1 (MQ=255)
                                           gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAAc                                                                                                           <  3:1661629/149‑1 (MQ=255)
                                           gCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAAc                                                                                                           <  4:815561/149‑1 (MQ=255)
                                            cTACTGGCGACCAGCACGGTGGCGGTGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg                                                                                                          <  4:1609305/149‑1 (MQ=255)
                                            cTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACg                                                                                                          >  1:1549317/1‑149 (MQ=255)
                                             tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGcc                                                                                                                        <  1:1127301/134‑1 (MQ=255)
                                             tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGcc                                                                                                                        >  2:1127301/1‑134 (MQ=255)
                                             tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACgg                                                                                                         <  1:152945/149‑1 (MQ=255)
                                             tACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACgg                                                                                                         >  3:1223951/1‑149 (MQ=255)
                                              aCTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGt                                                                                                        <  4:308897/149‑1 (MQ=255)
                                               cTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGtt                                                                                                       <  1:532156/149‑1 (MQ=255)
                                                tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                                                                      >  1:8903/1‑149 (MQ=255)
                                                tGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTa                                                                                                      >  4:1339795/1‑149 (MQ=255)
                                                 ggagACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTAt                                                                                                     <  4:1206535/146‑1 (MQ=255)
                                                 ggCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTAt                                                                                                     <  4:1506248/149‑1 (MQ=255)
                                                       cAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACt                                                                                               >  2:253804/1‑149 (MQ=255)
                                                       cAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACt                                                                                               >  1:462073/1‑149 (MQ=255)
                                                         gCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTTAGCAGCGCTCGGTACACATTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAACCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGc                                                                                             <  4:1030119/149‑1 (MQ=255)
                                                         gCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGc                                                                                             <  1:281256/149‑1 (MQ=255)
                                                         gCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGc                                                                                             >  2:890862/1‑149 (MQ=255)
                                                           aCGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGt                                                                                           <  3:202797/149‑1 (MQ=255)
                                                                  cGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatca                                                                                    <  4:47780/149‑1 (MQ=255)
                                                                   gATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcat                                                                                   >  1:712564/1‑149 (MQ=255)
                                                                   gATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcat                                                                                   >  1:784191/1‑149 (MQ=255)
                                                                    aTGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                                                                                  <  1:684644/149‑1 (MQ=255)
                                                                    aTGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                                                                                  <  4:149156/149‑1 (MQ=255)
                                                                    aTGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                                                                                  <  3:1670686/149‑1 (MQ=255)
                                                                    aTGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATaccatc                                                                                  <  2:650498/149‑1 (MQ=255)
                                                                       gTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATa                                                                               <  2:382825/149‑1 (MQ=255)
                                                                        tGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCCCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAg                                                                              >  4:1421135/1‑149 (MQ=255)
                                                                         gACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGt                                                                             >  4:1149272/1‑149 (MQ=255)
                                                                          aCGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGtt                                                                            <  1:610135/149‑1 (MQ=255)
                                                                             atatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaaa                                                                         <  2:124794/149‑1 (MQ=255)
                                                                             atatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaaa                                                                         <  4:739531/149‑1 (MQ=255)
                                                                             atatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTaaa                                                                         <  2:517414/149‑1 (MQ=255)
                                                                              tatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAAc                                                                        >  2:713773/1‑149 (MQ=255)
                                                                              tatCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAAc                                                                        <  4:898575/149‑1 (MQ=255)
                                                                               atCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGtt                                                                            <  4:1529491/144‑1 (MQ=255)
                                                                                tCGCGTAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaa                                                                      >  2:1420500/1‑149 (MQ=255)
                                                                                tCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaa                                                                      >  3:403453/1‑149 (MQ=255)
                                                                                    gcAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCa                                                                  <  3:270236/149‑1 (MQ=255)
                                                                                    gcAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCa                                                                  >  3:924455/1‑149 (MQ=255)
                                                                                     cAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAg                                                                 <  4:1223951/149‑1 (MQ=255)
                                                                                     cAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAg                                                                 <  4:924455/149‑1 (MQ=255)
                                                                                     cAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAg                                                                 <  3:47725/149‑1 (MQ=255)
                                                                                       gggAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGa                                                               <  3:696680/149‑1 (MQ=255)
                                                                                        ggAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGaa                                                              <  3:1555286/149‑1 (MQ=255)
                                                                                            aCGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAAtata                                                          <  2:784191/149‑1 (MQ=255)
                                                                                               ccTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCt                                                       >  2:613311/1‑149 (MQ=255)
                                                                                                cTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACaaa                                                                     <  4:473420/134‑1 (MQ=255)
                                                                                                cTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTg                                                      >  2:811067/1‑149 (MQ=255)
                                                                                                 ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc                                                     >  3:688614/1‑149 (MQ=255)
                                                                                                 ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc                                                     >  3:1589106/1‑149 (MQ=255)
                                                                                                 ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc                                                     <  3:1493672/149‑1 (MQ=255)
                                                                                                 ttCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTgc                                                     <  3:1339795/149‑1 (MQ=255)
                                                                                                         gcagcGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATg                                                                                                    <  4:1254289/94‑1 (MQ=255)
                                                                                                         gcagcGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATg                                                                                                    >  3:1254289/1‑94 (MQ=255)
                                                                                                           agcGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCtt                                           <  3:188878/149‑1 (MQ=255)
                                                                                                            gcgcCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTa                                          <  2:1157205/149‑1 (MQ=255)
                                                                                                              gcCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTAtt                                        <  2:462073/149‑1 (MQ=255)
                                                                                                              gcCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTAtt                                        <  2:210395/149‑1 (MQ=255)
                                                                                                                       acTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACaa                                               <  2:1079690/133‑1 (MQ=255)
                                                                                                                       acTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACaa                                               >  1:1079690/1‑133 (MQ=255)
                                                                                                                        cTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTa                              <  3:569352/149‑1 (MQ=255)
                                                                                                                         tttCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCaa                                                                                                 >  2:53606/1‑81 (MQ=255)
                                                                                                                         tttCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATggg                                >  2:720068/1‑146 (MQ=255)
                                                                                                                            cAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGa                          <  4:283015/149‑1 (MQ=255)
                                                                                                                            cAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGa                          <  1:1285147/149‑1 (MQ=255)
                                                                                                                                aaCTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTa                                          <  2:1175944/129‑1 (MQ=255)
                                                                                                                                aaCTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTa                                          >  1:1175944/1‑129 (MQ=255)
                                                                                                                                aaCTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCAc                      <  4:837913/149‑1 (MQ=255)
                                                                                                                                aaCTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCAc                      >  4:200031/1‑149 (MQ=255)
                                                                                                                                  cTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATatcatc                                                                                  >  4:648257/1‑87 (MQ=255)
                                                                                                                                       tGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGgcg               >  4:1415582/1‑149 (MQ=255)
                                                                                                                                       tGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGgcg               >  4:1327023/1‑149 (MQ=255)
                                                                                                                                          atttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTg            >  1:173236/3‑149 (MQ=255)
                                                                                                                                           tttGAGCGCGGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGATTATGCAACTGCGTAATATCATCATAGTTAAACAAATAAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGt           >  2:1716476/2‑149 (MQ=255)
                                                                                                                                               aGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGAt       >  4:523268/1‑149 (MQ=255)
                                                                                                                                                   ccgcGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCAcc   <  2:1195978/149‑1 (MQ=255)
                                                                                                                                                    cgcgCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCACCTGCGTAATATCATCATAGTTAAACAAATCAGGAAGATAGCTGCGCACAATCTTATTTGCATGGGTAGCGATCACGCTGGCTCTGTCGATCACCt  <  4:1362935/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGCACATCGGCCGCCAGCAGAATATGATCTTTCGTCACCGTGCTACTGGCGACCAGCACGGTGGCGATGGTGACGATATCGCGCAGGGAAACGCCTTCGGTCAGCAGCGCCCGGTACACTTTCAGCAACTGGCTGTAATTGAGCGCCGCGCTCAAATCTTCCGCCAGACGCGGTGCCGTCGACGATAAACGGTTATGCAACTGCGTAATATCATCATAGTTAAACAAATCAGGAATATAGCTGCGCACAATCTTATTTACATGGGTAGCGATCACGCTGGCGCTGTCGATCACCT  >  NC_000913/248720‑249016

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