Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 385,815 G→T 92.7% 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,8150GT92.7% 727.7 / 50.0 246pseudogene (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 (7/11);  new base T (120/108);  total (127/119)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 248 total aligned reads displayed.

TTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGG  >  CP009273/385666‑385960
                                                                                                                                                     |                                                                                                                                                 
   aTGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCtata                                                                                                                                                <  2:3106543/150‑1 (MQ=255)
      gATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATc                                                                                                                                             <  2:704647/150‑1 (MQ=255)
      gATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATc                                                                                                                                             <  1:3532826/150‑1 (MQ=255)
          aaCTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGt                                                                                                                                         <  2:4630458/150‑1 (MQ=255)
          aaCTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGt                                                                                                                                         <  2:4072909/150‑1 (MQ=255)
           aCTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGtt                                                                                                                                        <  2:4123501/150‑1 (MQ=255)
            cTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTa                                                                                                                                       <  1:3199886/150‑1 (MQ=255)
            cTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTa                                                                                                                                       <  1:1703253/150‑1 (MQ=255)
               attaCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAg                                                                                                                                    >  1:2308573/1‑150 (MQ=255)
               attaCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAg                                                                                                                                    >  1:1974326/1‑150 (MQ=255)
               attaCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAg                                                                                                                                    <  1:3826582/150‑1 (MQ=255)
               attaCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACAATGATACCTGGACTTATATCTATAATCCCGTTACAg                                                                                                                                    >  1:4044552/1‑150 (MQ=255)
                ttaCGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGc                                                                                                                                   >  1:3599052/1‑150 (MQ=255)
                    gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACtt                                                                                                                               >  1:162248/1‑150 (MQ=255)
                    gTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACtt                                                                                                                               >  1:3845221/1‑150 (MQ=255)
                           aggaggGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTggg                                                                                                                        <  2:361167/150‑1 (MQ=255)
                           aggaggGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTggg                                                                                                                        >  2:323881/1‑150 (MQ=255)
                           aggaggGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTggg                                                                                                                        <  2:3850580/150‑1 (MQ=255)
                            ggaggGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGa                                                                                                                       <  1:100459/150‑1 (MQ=255)
                               gggATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCa                                                                                                                    >  2:2717142/1‑150 (MQ=255)
                                ggATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCat                                                                                                                   <  2:1556638/150‑1 (MQ=255)
                                ggATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTAAAGCACTTGATTGGGATCat                                                                                                                   >  2:4054921/1‑150 (MQ=255)
                                 gATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCata                                                                                                                  >  1:1999719/1‑150 (MQ=255)
                                 gATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCata                                                                                                                  >  1:469259/1‑150 (MQ=255)
                                    gACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATgg                                                                                                               >  2:1880856/1‑150 (MQ=255)
                                     aCAAATACACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGa                                                                                                              >  2:1513010/1‑150 (MQ=255)
                                       aaaTGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATa                                                                                                            >  1:2818102/1‑150 (MQ=255)
                                        aaTGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTGGATTGGGATCATATGGATaa                                                                                                           >  1:143337/1‑150 (MQ=255)
                                          tGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAagg                                                                                                         <  1:1018791/150‑1 (MQ=255)
                                             aCTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAaggcag                                                                                                      <  1:2165396/150‑1 (MQ=255)
                                                ccTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAggcaggca                                                                                                   >  2:261658/1‑150 (MQ=255)
                                                 cTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCAc                                                                                                  <  1:3852649/150‑1 (MQ=255)
                                                  tGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACt                                                                                                 >  1:323402/1‑150 (MQ=255)
                                                  tGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACt                                                                                                 <  2:2968253/150‑1 (MQ=255)
                                                  tGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACt                                                                                                 <  1:1754063/150‑1 (MQ=255)
                                                   ggAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTg                                                                                                <  1:3394044/150‑1 (MQ=255)
                                                     aaaaaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTa                                                                                              >  1:4474990/1‑150 (MQ=255)
                                                     aaaaaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCGATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATCATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTTTa                                                                                              <  2:2936273/150‑1 (MQ=255)
                                                      aaaaGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAg                                                                                             <  1:1592960/150‑1 (MQ=255)
                                                             aCTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCaa                                                                                      >  1:2706491/1‑150 (MQ=255)
                                                               ttGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACt                                                                                    <  1:4268813/150‑1 (MQ=255)
                                                                tGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTa                                                                                   <  1:1178965/150‑1 (MQ=255)
                                                                  tagtaTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATg                                                                                 <  2:2308573/150‑1 (MQ=255)
                                                                  tagtaTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATg                                                                                 >  1:1233873/1‑150 (MQ=255)
                                                                   agtaTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGt                                                                                >  2:1225673/1‑150 (MQ=255)
                                                                    gtaTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTa                                                                               >  2:4513135/1‑150 (MQ=255)
                                                                      aTCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTaaa                                                                             >  2:3456756/1‑150 (MQ=255)
                                                                        cacaATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACt                                                                           <  2:2643983/150‑1 (MQ=255)
                                                                        cacaATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACt                                                                           >  1:3931128/1‑150 (MQ=255)
                                                                                gCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGt                                                                   >  2:3568749/1‑150 (MQ=255)
                                                                                   ttGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGtttt                                                                <  1:256280/150‑1 (MQ=255)
                                                                                   ttGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGtttt                                                                <  1:256104/150‑1 (MQ=255)
                                                                                    tGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTg                                                               >  1:213837/1‑150 (MQ=255)
                                                                                     ggCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGt                                                              >  2:2749503/1‑150 (MQ=255)
                                                                                     ggCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGt                                                              <  2:2203753/150‑1 (MQ=255)
                                                                                      gCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTc                                                             <  2:3965316/150‑1 (MQ=255)
                                                                                      gCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTc                                                             <  2:1640354/150‑1 (MQ=255)
                                                                                        ttGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTa                                                           <  1:2725500/150‑1 (MQ=255)
                                                                                         tGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTAt                                                          <  1:2176442/150‑1 (MQ=255)
                                                                                         tGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTAt                                                          >  1:4393842/1‑150 (MQ=255)
                                                                                           cACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTa                                                        >  1:4614752/1‑150 (MQ=255)
                                                                                           cACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTa                                                        <  2:1905069/150‑1 (MQ=255)
                                                                                             ctctGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTaca                                                      >  1:1799657/1‑150 (MQ=255)
                                                                                                  gCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTacaacaac                                                 >  1:4320804/1‑150 (MQ=255)
                                                                                                        aCCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAAt                                           <  2:610014/150‑1 (MQ=255)
                                                                                                         ccTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATg                                          >  2:585301/1‑150 (MQ=255)
                                                                                                          ctctCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATgg                                         <  1:4119489/150‑1 (MQ=255)
                                                                                                             tcAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGa                                      >  1:2868194/1‑150 (MQ=255)
                                                                                                             tcAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGa                                      >  1:2018484/1‑150 (MQ=255)
                                                                                                              cAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAt                                     <  2:3892993/150‑1 (MQ=255)
                                                                                                               aGGCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAtt                                    <  2:4273810/150‑1 (MQ=255)
                                                                                                                ggCTAATGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGAttt                                   <  2:3410017/150‑1 (MQ=255)
                                                                                                                   tAATGCATCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGa                                <  2:945367/150‑1 (MQ=255)
                                                                                                                    aaTGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGAt                               >  1:302498/1‑150 (MQ=255)
                                                                                                                    aaTGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGAt                               >  1:3463821/1‑150 (MQ=255)
                                                                                                                    aaTGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGAt                               >  1:1864191/1‑150 (MQ=255)
                                                                                                                      tGCAGCAAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCa                             <  2:1529226/150‑1 (MQ=255)
                                                                                                                        cagcaAACTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGt                           <  2:2484002/150‑1 (MQ=255)
                                                                                                                            aaaCTATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCtt                       >  2:20318/1‑150 (MQ=255)
                                                                                                                                tATGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAAc                   <  1:432369/150‑1 (MQ=255)
                                                                                                                                  tGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACgg                 >  1:2545439/1‑150 (MQ=255)
                                                                                                                                  tGATACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACgg                 >  1:4265636/1‑150 (MQ=255)
                                                                                                                                    aTACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCg               <  1:3526213/150‑1 (MQ=255)
                                                                                                                                    aTACCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCg               <  2:4386080/150‑1 (MQ=255)
                                                                                                                                      aCCTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGAt             <  2:2573319/150‑1 (MQ=255)
                                                                                                                                        cTGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAc           >  2:3723393/1‑150 (MQ=255)
                                                                                                                                         tGGACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAcc          <  2:3450169/150‑1 (MQ=255)
                                                                                                                                          ggACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCg         >  2:890845/1‑150 (MQ=255)
                                                                                                                                          ggACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCg         <  1:2337586/150‑1 (MQ=255)
                                                                                                                                          ggACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCg         >  1:1194988/1‑150 (MQ=255)
                                                                                                                                           gACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGt        >  1:4388154/1‑150 (MQ=255)
                                                                                                                                           gACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGt        <  1:3626930/150‑1 (MQ=255)
                                                                                                                                           gACTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGt        >  1:725357/1‑150 (MQ=255)
                                                                                                                                            aCTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGtt       >  2:528823/1‑150 (MQ=255)
                                                                                                                                            aCTTATATATATAATCCCGTTAGAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAAATTTTGTGGTTTTGTATATTACATCAACACCAATGTTGATTTCGATCAGTCCTTTAACGGCGATACCGtt       <  2:4268690/150‑1 (MQ=255)
                                                                                                                                             cTTATATCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTa      <  1:243765/150‑1 (MQ=255)
                                                                                                                                               tataTCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAAc    <  1:1493994/150‑1 (MQ=255)
                                                                                                                                               tataTCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAAc    <  2:2074206/150‑1 (MQ=255)
                                                                                                                                                atatCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACg   <  2:816845/150‑1 (MQ=255)
                                                                                                                                                 tatCTATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACgg  <  2:416913/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCGTTAACGG  >  CP009273/385666‑385960

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: