Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 385,815 G→T 59.9% pseudogene (109/1458 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273385,8150GT59.9% 86.2 / inf 257pseudogene (109/1458 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base G (45/58);  new base T (80/74);  total (125/132)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 259 total aligned reads displayed.

TTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGGCACG  >  CP009273/385666‑385964
                                                                                                                                                     |                                                                                                                                                     
ttGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCg                                                                                                                                                       <  1:4058765/150‑1 (MQ=255)
     ggATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAAt                                                                                                                                                  <  2:435716/150‑1 (MQ=255)
     ggATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAAt                                                                                                                                                  <  2:711623/150‑1 (MQ=255)
       aTCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATcc                                                                                                                                                <  2:2260441/150‑1 (MQ=255)
       aTCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATcc                                                                                                                                                <  2:3123542/150‑1 (MQ=255)
        tCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATccc                                                                                                                                               <  2:1734421/150‑1 (MQ=255)
         cAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCg                                                                                                                                              >  1:1682134/1‑150 (MQ=255)
                 taCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCa                                                                                                                                      <  1:1472800/150‑1 (MQ=255)
                 taCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCa                                                                                                                                      <  2:1456847/150‑1 (MQ=255)
                   cGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACt                                                                                                                                    <  1:3939069/150‑1 (MQ=255)
                   cGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACt                                                                                                                                    >  2:3409139/1‑150 (MQ=255)
                    gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACtt                                                                                                                                   <  1:2387196/150‑1 (MQ=255)
                    gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACtt                                                                                                                                   >  2:2061341/1‑150 (MQ=255)
                      cccTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGa                                                                                                                                 >  2:607749/1‑150 (MQ=255)
                       ccTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGAt                                                                                                                                <  1:170618/150‑1 (MQ=255)
                       ccTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGAt                                                                                                                                <  2:1628216/150‑1 (MQ=255)
                        cTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGAtt                                                                                                                               <  1:670710/150‑1 (MQ=255)
                        cTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGAtt                                                                                                                               <  2:1686453/150‑1 (MQ=255)
                        cTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGAtt                                                                                                                               <  1:2248911/150‑1 (MQ=255)
                        cTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGAtt                                                                                                                               <  1:2249006/150‑1 (MQ=255)
                         tGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTg                                                                                                                              >  2:2183688/1‑150 (MQ=255)
                               gggATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCa                                                                                                                        >  1:1319555/1‑150 (MQ=255)
                                ggATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCat                                                                                                                       >  1:912142/1‑150 (MQ=255)
                                  aTGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCatat                                                                                                                     <  2:2106891/150‑1 (MQ=255)
                                         aTGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAag                                                                                                              <  1:866288/150‑1 (MQ=255)
                                             aCTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAaggcag                                                                                                          >  2:1770467/1‑150 (MQ=255)
                                                 cTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCAc                                                                                                      <  1:2299557/150‑1 (MQ=255)
                                                 cTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCAc                                                                                                      <  2:2045794/150‑1 (MQ=255)
                                                  tGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACt                                                                                                     >  2:126778/1‑150 (MQ=255)
                                                    gAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGt                                                                                                   <  2:820153/150‑1 (MQ=255)
                                                    gAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGt                                                                                                   >  1:1016243/1‑150 (MQ=255)
                                                    gAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGt                                                                                                   >  1:644141/1‑150 (MQ=255)
                                                     aaaaaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTa                                                                                                  <  2:1044941/150‑1 (MQ=255)
                                                     aaaaaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAATGCAGGCACTGTa                                                                                                  <  2:390658/150‑1 (MQ=255)
                                                        aaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGAt                                                                                               <  1:585714/150‑1 (MQ=255)
                                                             aCTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCaa                                                                                          <  1:3140220/150‑1 (MQ=255)
                                                              cTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAAc                                                                                         >  2:4019063/1‑150 (MQ=255)
                                                                tGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGTAACTa                                                                                       <  2:2162503/150‑1 (MQ=255)
                                                                tGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTa                                                                                       >  1:61186/1‑150 (MQ=255)
                                                                   agtaTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGt                                                                                    <  1:161928/150‑1 (MQ=255)
                                                                        cacaATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACt                                                                               >  2:3531701/1‑150 (MQ=255)
                                                                         acaATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACta                                                                              >  1:196439/1‑150 (MQ=255)
                                                                          caATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACtat                                                                             <  1:2206689/150‑1 (MQ=255)
                                                                             ttAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGt                                                                          >  2:1649449/1‑150 (MQ=255)
                                                                              tAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTg                                                                         >  2:58573/1‑150 (MQ=255)
                                                                               ggCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTgg                                                                        >  2:3176598/2‑150 (MQ=255)
                                                                               aGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTgg                                                                        <  1:1847568/150‑1 (MQ=255)
                                                                                   ttGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATATATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGtttt                                                                    <  1:3675729/150‑1 (MQ=255)
                                                                                      gCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTTGTTTTGTc                                                                 <  2:840781/150‑1 (MQ=255)
                                                                                       cTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCt                                                                <  2:1086961/150‑1 (MQ=255)
                                                                                        ttGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTa                                                               <  2:2920885/150‑1 (MQ=255)
                                                                                        ttGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTa                                                               <  2:2448491/150‑1 (MQ=255)
                                                                                          gCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTAtt                                                             <  1:3465851/150‑1 (MQ=255)
                                                                                            aCTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTac                                                           <  1:3022909/150‑1 (MQ=255)
                                                                                               ctGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaac                                                        <  2:2836674/150‑1 (MQ=255)
                                                                                               ctGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaac                                                        <  1:436481/150‑1 (MQ=255)
                                                                                                 ggCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaacaa                                                      >  2:3256272/1‑150 (MQ=255)
                                                                                                    tATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAacac                                                   >  1:385232/1‑150 (MQ=255)
                                                                                                    tATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACATATATCTATAATCCCGTTACAGAACTTGATTGGGATCATATGGATAATGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGGTTTGTCTATTACAACAacac                                                   >  2:95182/1‑150 (MQ=255)
                                                                                                        aCCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAAt                                               >  1:1798729/1‑150 (MQ=255)
                                                                                                        aCCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAAt                                               >  1:3897071/1‑150 (MQ=255)
                                                                                                         ccTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATg                                              >  1:4170082/1‑150 (MQ=255)
                                                                                                         ccTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATg                                              >  2:4108759/1‑150 (MQ=255)
                                                                                                         ccTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATg                                              <  2:2153940/150‑1 (MQ=255)
                                                                                                         ccTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATg                                              <  2:1063922/150‑1 (MQ=255)
                                                                                                             tcAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGa                                          >  2:3510197/1‑150 (MQ=255)
                                                                                                                ggCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAttt                                       <  2:4002706/150‑1 (MQ=255)
                                                                                                                ggCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAttt                                       >  1:291158/1‑150 (MQ=255)
                                                                                                                ggCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAttt                                       >  1:3724370/1‑150 (MQ=255)
                                                                                                                  cTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCg                                     >  1:3729241/1‑150 (MQ=255)
                                                                                                                   tAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGa                                    >  1:2541302/1‑150 (MQ=255)
                                                                                                                   tAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGa                                    <  2:923564/150‑1 (MQ=255)
                                                                                                                   tAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGa                                    >  2:1897760/1‑150 (MQ=255)
                                                                                                                    aaTGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGAt                                   >  2:2396112/1‑150 (MQ=255)
                                                                                                                    aaTGCAGCAAACTATGATACCTGGACTTATATCTAAAATCCCGTTACAGCAATTGATTGGGATCATATGTATAAGGCAGGCACGGAAGATGGCAACGATGTAAACAATAGTGGTTTTGTCTATTAAAACAACAACAATGGTGATTTCGAt                                   >  2:1983905/1‑150 (MQ=255)
                                                                                                                     aTGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATc                                  <  2:2800042/150‑1 (MQ=255)
                                                                                                                     aTGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATc                                  >  1:1497948/1‑150 (MQ=255)
                                                                                                                      tGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCa                                 <  1:482934/150‑1 (MQ=255)
                                                                                                                         agcaAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTc                              <  2:1395826/150‑1 (MQ=255)
                                                                                                                         agcaAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATGACAACAACACCAATGGTGATCTCGATCAGTc                              >  2:135596/1‑150 (MQ=255)
                                                                                                                          gcaAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTcc                             >  2:3828099/1‑150 (MQ=255)
                                                                                                                            aaaCTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCtt                           <  2:3349196/150‑1 (MQ=255)
                                                                                                                              aCTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTa                         >  1:2692422/1‑150 (MQ=255)
                                                                                                                              aCTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTa                         >  2:4121133/1‑150 (MQ=255)
                                                                                                                               cTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTaa                        <  1:1171927/150‑1 (MQ=255)
                                                                                                                                tATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAAc                       <  1:2821056/150‑1 (MQ=255)
                                                                                                                                 aTGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACg                      <  2:1557963/150‑1 (MQ=255)
                                                                                                                                     tACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGa                  <  1:2026557/150‑1 (MQ=255)
                                                                                                                                     tACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGa                  <  1:4126370/150‑1 (MQ=255)
                                                                                                                                        cTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAc               >  1:2467184/1‑150 (MQ=255)
                                                                                                                                        cTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAc               >  2:1515074/1‑150 (MQ=255)
                                                                                                                                         tGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAcc              >  1:1129334/1‑150 (MQ=255)
                                                                                                                                          ggACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCg             >  2:414586/1‑150 (MQ=255)
                                                                                                                                           gACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGt            >  1:3298824/1‑150 (MQ=255)
                                                                                                                                            aCTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGtt           <  1:3584384/150‑1 (MQ=255)
                                                                                                                                            aCTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGtt           >  2:1822025/1‑150 (MQ=255)
                                                                                                                                             cTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTa          >  2:2558871/1‑150 (MQ=255)
                                                                                                                                             cTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTa          >  2:2672691/1‑150 (MQ=255)
                                                                                                                                                atatCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACg       >  2:2663465/1‑150 (MQ=255)
                                                                                                                                                  atCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGGc     >  2:1527681/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGGCACG  >  CP009273/385666‑385964

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