Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 385,815 G→T 89.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,8150GT89.9% 589.2 / 62.8 216pseudogene (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 (8/14);  new base T (112/82);  total (120/96)
Fisher's exact test for biased strand distribution p-value = 7.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 216 total aligned reads displayed.

TGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGGC  >  CP009273/385667‑385961
                                                                                                                                                    |                                                                                                                                                  
     gATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATc                                                                                                                                              <  1:1463759/150‑1 (MQ=255)
     gATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATc                                                                                                                                              >  1:2320371/1‑150 (MQ=255)
       tCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATccc                                                                                                                                            >  2:2023770/1‑150 (MQ=255)
       tCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATccc                                                                                                                                            >  2:2539587/1‑150 (MQ=255)
        cAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCg                                                                                                                                           >  2:3826535/1‑150 (MQ=255)
        cAAATTATTACGTCACTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGATATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCg                                                                                                                                           <  2:69057/150‑1 (MQ=255)
          aCTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGtt                                                                                                                                         >  1:2427519/1‑150 (MQ=255)
               ttaCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGc                                                                                                                                    >  2:1674695/1‑150 (MQ=255)
                taCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCa                                                                                                                                   >  2:40445/1‑150 (MQ=255)
                   gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACtt                                                                                                                                >  2:1167927/1‑150 (MQ=255)
                   gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACtt                                                                                                                                <  1:4154276/150‑1 (MQ=255)
                   gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACtt                                                                                                                                >  1:3129077/1‑150 (MQ=255)
                   gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACtt                                                                                                                                >  1:805834/1‑150 (MQ=255)
                    tCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTg                                                                                                                               >  2:82206/1‑150 (MQ=255)
                     cccTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGa                                                                                                                              >  1:787158/1‑150 (MQ=255)
                      ccTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGAt                                                                                                                             >  2:2866635/1‑150 (MQ=255)
                        tGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTg                                                                                                                           <  2:600716/150‑1 (MQ=255)
                        tGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTg                                                                                                                           <  2:600811/150‑1 (MQ=255)
                         gaggagGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTgg                                                                                                                          >  2:3521390/1‑150 (MQ=255)
                         gaggagGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTgg                                                                                                                          >  2:402202/1‑150 (MQ=255)
                            gaggGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGAt                                                                                                                       <  1:3318253/150‑1 (MQ=255)
                              gggATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACGTGATTGGGATCa                                                                                                                     >  2:4100266/1‑150 (MQ=255)
                                 aTGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCatat                                                                                                                  >  2:891805/1‑150 (MQ=255)
                                   gACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATgg                                                                                                                >  2:1425412/1‑150 (MQ=255)
                                     cAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGAt                                                                                                              <  1:4319232/150‑1 (MQ=255)
                                      aaaTGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATa                                                                                                             <  2:2393268/150‑1 (MQ=255)
                                        aTGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAag                                                                                                           <  2:121885/150‑1 (MQ=255)
                                         tGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAagg                                                                                                          <  2:3198721/150‑1 (MQ=255)
                                         tGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAagg                                                                                                          <  1:919138/150‑1 (MQ=255)
                                           cACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAaggca                                                                                                        >  1:3852470/1‑150 (MQ=255)
                                                 tGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACt                                                                                                  >  2:1718682/1‑150 (MQ=255)
                                                 tGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACt                                                                                                  >  2:4453222/1‑150 (MQ=255)
                                                  ggAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTg                                                                                                 <  1:2466071/150‑1 (MQ=255)
                                                  ggAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTg                                                                                                 <  1:3735366/150‑1 (MQ=255)
                                                  ggAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTCATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTg                                                                                                 >  2:465901/1‑150 (MQ=255)
                                                    aaaaaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTa                                                                                               >  1:2442151/1‑150 (MQ=255)
                                                      aaaGAAACTTGTAGTATCACAATTAGCATTGGCGTGCACTCTGGCTATCACCTCTCAGGCTAAAGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGAGTGTGATCATATGGATAAGGCAGGCACAGTAGa                                                                                             >  1:86729/1‑150 (MQ=255)
                                                        aGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATg                                                                                           <  1:1389645/150‑1 (MQ=255)
                                                        aGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATg                                                                                           >  2:3362141/1‑150 (MQ=255)
                                                        aGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATg                                                                                           <  2:3596762/150‑1 (MQ=255)
                                                        aGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATg                                                                                           >  1:2961202/1‑150 (MQ=255)
                                                        aGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATg                                                                                           >  1:2020401/1‑150 (MQ=255)
                                                         gAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATgg                                                                                          >  1:1108287/1‑150 (MQ=255)
                                                          aaaCTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGc                                                                                         >  1:3673319/1‑150 (MQ=255)
                                                          aaaCTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGc                                                                                         <  1:2866635/150‑1 (MQ=255)
                                                             cTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAAc                                                                                      <  1:1406815/150‑1 (MQ=255)
                                                             cTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAAc                                                                                      <  1:4204634/150‑1 (MQ=255)
                                                              ttGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACt                                                                                     <  2:2178151/150‑1 (MQ=255)
                                                                 tagtaTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATg                                                                                  >  1:1885996/1‑150 (MQ=255)
                                                                         caATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACtat                                                                          <  1:3376648/150‑1 (MQ=255)
                                                                         caATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACtat                                                                          <  2:4013509/150‑1 (MQ=255)
                                                                         caATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACtat                                                                          <  2:4337492/150‑1 (MQ=255)
                                                                                cATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGtt                                                                   >  2:2493763/1‑150 (MQ=255)
                                                                                 aTTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGttt                                                                  <  2:1240472/150‑1 (MQ=255)
                                                                                 aTTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGttt                                                                  <  2:1240491/150‑1 (MQ=255)
                                                                                   tGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTg                                                                <  1:2532455/150‑1 (MQ=255)
                                                                                    ggCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGt                                                               >  1:963127/1‑150 (MQ=255)
                                                                                       ttGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTa                                                            <  2:1959026/150‑1 (MQ=255)
                                                                                        tGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTAt                                                           <  1:600436/150‑1 (MQ=255)
                                                                                         gCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTAtt                                                          <  1:4379338/150‑1 (MQ=255)
                                                                                           aCTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTac                                                        >  2:1671774/1‑150 (MQ=255)
                                                                                             tctGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGTATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaa                                                      >  1:3951790/1‑150 (MQ=255)
                                                                                              ctGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaac                                                     <  1:3075351/150‑1 (MQ=255)
                                                                                              ctGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaac                                                     <  1:2158560/150‑1 (MQ=255)
                                                                                               tGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaaca                                                    >  2:3237424/1‑150 (MQ=255)
                                                                                                   tATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAacac                                                >  1:936237/1‑150 (MQ=255)
                                                                                                   tATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAacac                                                >  1:454941/1‑150 (MQ=255)
                                                                                                   tATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAacac                                                <  2:4350351/150‑1 (MQ=255)
                                                                                                   tATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAacac                                                <  2:3022821/150‑1 (MQ=255)
                                                                                                    aTCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACAcc                                               >  1:2873212/1‑150 (MQ=255)
                                                                                                          tctcAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGt                                         <  2:407671/150‑1 (MQ=255)
                                                                                                           ctcAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTg                                        >  2:396053/1‑150 (MQ=255)
                                                                                                            tcAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGa                                       >  1:4590906/1‑150 (MQ=255)
                                                                                                              aGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAtt                                     <  2:2228987/150‑1 (MQ=255)
                                                                                                              aGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAtt                                     >  1:3617938/1‑150 (MQ=255)
                                                                                                              aGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAtt                                     <  1:503773/150‑1 (MQ=255)
                                                                                                               ggCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAttt                                    <  1:1212080/150‑1 (MQ=255)
                                                                                                                gCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTc                                   <  2:2367232/150‑1 (MQ=255)
                                                                                                                 cTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCg                                  >  1:1922430/1‑150 (MQ=255)
                                                                                                                  tAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGa                                 >  2:3255416/1‑150 (MQ=255)
                                                                                                                   aaTGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGAt                                >  2:3344347/1‑150 (MQ=255)
                                                                                                                     tGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCa                              >  2:3943717/1‑150 (MQ=255)
                                                                                                                       cagcaAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGt                            <  1:4484030/150‑1 (MQ=255)
                                                                                                                         gcaAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTcc                          <  1:45430/150‑1 (MQ=255)
                                                                                                                          caAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCt                         >  2:856044/1‑150 (MQ=255)
                                                                                                                           aaaCTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCtt                        >  1:1446413/1‑150 (MQ=255)
                                                                                                                            aaCTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCttt                       <  1:2539587/150‑1 (MQ=255)
                                                                                                                               tATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAAc                    >  2:3755076/1‑150 (MQ=255)
                                                                                                                                  gATACATTGACTTATATCTAGAATCGCGTTAAAGCACTTGATTGGGTTCATATGGATAAGTCATGCACTGTAGATGTCAACTATGTAAACTATAGTGGTTTTGTCTATTTCAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGc                 <  1:1744804/150‑1 (MQ=255)
                                                                                                                                   aTACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCg                >  2:2520970/1‑150 (MQ=255)
                                                                                                                                     aCCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGAt              >  2:1765019/1‑150 (MQ=255)
                                                                                                                                      ccTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATa             >  1:166625/1‑150 (MQ=255)
                                                                                                                                        tGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAcc           >  2:4240779/1‑150 (MQ=255)
                                                                                                                                        tGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAcc           <  1:4236842/150‑1 (MQ=255)
                                                                                                                                         ggACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCg          <  1:1295172/150‑1 (MQ=255)
                                                                                                                                          gACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGt         <  1:3821642/150‑1 (MQ=255)
                                                                                                                                          gACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGt         >  1:792180/1‑150 (MQ=255)
                                                                                                                                              tataTCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAAc     >  1:4224614/1‑150 (MQ=255)
                                                                                                                                               atatCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACg    <  2:787158/150‑1 (MQ=255)
                                                                                                                                                 atCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGGc  >  1:826570/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGGC  >  CP009273/385667‑385961

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: