Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,837,177 C→T 30.2% intergenic (‑28/‑147) yicO ← / → ade putative adenine permease/cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,837,1770CT30.2% 587.2 / inf 368intergenic (‑28/‑147)yicO/adeputative adenine permease/cryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (106/151);  new base T (57/54);  total (163/205)
Fisher's exact test for biased strand distribution p-value = 8.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 375 total aligned reads displayed.

AAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAATAT  >  CP009273/3837028‑3837325
                                                                                                                                                     |                                                                                                                                                    
 aaaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATcc                                                                                                                                                     <  2:3325245/150‑1 (MQ=255)
  aaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATtcc                                                                                                                                                    <  1:2515016/150‑1 (MQ=255)
    aaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCtt                                                                                                                                                  <  2:1086497/150‑1 (MQ=255)
     aGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTg                                                                                                                                                 >  1:4381246/1‑150 (MQ=255)
     aGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTg                                                                                                                                                 >  2:482705/1‑150 (MQ=255)
         gTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCttt                                                                                                                                             <  1:3129807/150‑1 (MQ=255)
          tCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCtttt                                                                                                                                            <  1:914759/150‑1 (MQ=255)
          tCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCtttt                                                                                                                                            >  2:938648/1‑150 (MQ=255)
          tCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCACTTGCtttt                                                                                                                                            >  2:2318046/1‑150 (MQ=255)
            aTCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGt                                                                                                                                          >  1:1929165/1‑150 (MQ=255)
              ccccGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTaa                                                                                                                                        <  1:3842120/150‑1 (MQ=255)
              ccccGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCATTGCTTTTGTaa                                                                                                                                        >  2:3105368/1‑150 (MQ=255)
                cctcaataaatTCTTTTCGGACGGTTGTCCCATTCTGATGTTGTTTAAATAATCTAGAAAGAGTCCCTTATTCATTAATAACGTAATCTTTATTTTCATTATTCATTTTTTTGTCCATTGGAGGATATTTAATCACTTTCTTTTGTaaaa                                                                                                                                      <  2:762670/147‑1 (MQ=255)
                  gCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACa                                                                                                                                    >  1:4696995/1‑150 (MQ=255)
                  gCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACa                                                                                                                                    >  2:286982/1‑150 (MQ=255)
                     atcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAAt                                                                                                                                 <  1:1000397/150‑1 (MQ=255)
                     atcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAAt                                                                                                                                 >  2:3786468/1‑150 (MQ=255)
                      tcaatTCTGTGCGGACGGTGGGCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCCTTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATg                                                                                                                                <  2:3685673/150‑1 (MQ=255)
                       caatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCAAGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGa                                                                                                                               <  1:2864059/150‑1 (MQ=255)
                        aatTCTGTGCGGACGGTGGTCCCATGCTGATGTAGTTTAAATAATCGCTAAAGCGTCCCTGATACAGTACTAACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGac                                                                                                                              <  2:2705329/150‑1 (MQ=255)
                        aatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGac                                                                                                                              >  1:852175/1‑150 (MQ=255)
                         atTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGaca                                                                                                                             <  1:3332795/150‑1 (MQ=255)
                          ttCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGacaa                                                                                                                            >  1:1410766/1‑150 (MQ=255)
                             tgtgCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGAcaacaa                                                                                                                         <  1:3591713/150‑1 (MQ=255)
                                  ggACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGtt                                                                                                                    >  2:3402615/1‑150 (MQ=255)
                                  ggACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGtt                                                                                                                    <  1:1797328/150‑1 (MQ=255)
                                    aCGGTGGTCCCATGCTGAGGTATTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCa                                                                                                                  <  2:297073/150‑1 (MQ=255)
                                        tggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGa                                                                                                              <  2:4883518/150‑1 (MQ=255)
                                        tggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGa                                                                                                              >  1:4040246/1‑150 (MQ=255)
                                        tggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATGCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGa                                                                                                              <  2:4196740/150‑1 (MQ=255)
                                            cccATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGAGTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAAt                                                                                                          <  2:3090695/150‑1 (MQ=255)
                                             ccATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAAtt                                                                                                         >  1:1668770/1‑150 (MQ=255)
                                             ccATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAAtt                                                                                                         <  1:5106603/150‑1 (MQ=255)
                                               aTGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTAt                                                                                                       <  1:1199214/150‑1 (MQ=255)
                                               aTGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTAt                                                                                                       >  2:4016427/1‑150 (MQ=255)
                                                  cTGAGGTAGTTTAAATACTCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATcac                                                                                                    <  1:2673630/150‑1 (MQ=255)
                                                  cTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATcac                                                                                                    >  1:3627165/1‑150 (MQ=255)
                                                     aGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACAtt                                                                                                 >  1:458276/1‑150 (MQ=255)
                                                               aataatCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGcac                                                                                       >  1:3714142/1‑150 (MQ=255)
                                                                 taatCTCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTAATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCACTTGCACAt                                                                                     <  1:1078290/150‑1 (MQ=255)
                                                                 taatCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACAt                                                                                     >  1:19617/1‑150 (MQ=255)
                                                                   atCGCGATAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTa                                                                                   <  2:912262/150‑1 (MQ=255)
                                                                   atCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTa                                                                                   >  1:586853/1‑150 (MQ=255)
                                                                     cgctAAAGCGTCCATGATTCAAAACTCACGTAATCGGTATTGTCATAATTCAATTTTTTGTCCATTGTAGGATATTTAATCCATTGCTTTTGTAAAACAAATGACAACAACCGTTAACTGATAATTATCAAATTACAATTTCACaataat                                                                                 >  2:4156780/1‑150 (MQ=255)
                                                                      gcgAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATg                                                                                <  1:4631391/150‑1 (MQ=255)
                                                                       cgAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATgg                                                                               <  1:2830254/150‑1 (MQ=255)
                                                                       cgAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATGGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATgg                                                                               <  2:557339/150‑1 (MQ=255)
                                                                        gAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGa                                                                              <  2:1515052/150‑1 (MQ=255)
                                                                         aaaGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGat                                                                             >  1:4244762/1‑150 (MQ=255)
                                                                          aaGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGata                                                                            <  1:4893606/150‑1 (MQ=255)
                                                                          aaGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGata                                                                            <  2:1320550/150‑1 (MQ=255)
                                                                              gTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCTCATTAATGGATATTCt                                                                        >  2:2142778/1‑150 (MQ=255)
                                                                              gTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCt                                                                        >  2:1294832/1‑150 (MQ=255)
                                                                               tCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCtt                                                                       <  2:259195/150‑1 (MQ=255)
                                                                                cccTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCttt                                                                      <  1:2921817/150‑1 (MQ=255)
                                                                                     aTTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTtaataa                                                                 >  2:318807/1‑150 (MQ=255)
                                                                                      ttCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGAGTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTaataat                                                                >  2:1045130/1‑150 (MQ=255)
                                                                                       tCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtc                                                               <  2:2453042/150‑1 (MQ=255)
                                                                                        cATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct                                                              >  1:2473015/1‑150 (MQ=255)
                                                                                        cATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct                                                              <  2:1988542/150‑1 (MQ=255)
                                                                                          ttACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTcg                                                            <  1:1324528/150‑1 (MQ=255)
                                                                                           tACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTcgc                                                           <  1:4960964/150‑1 (MQ=255)
                                                                                            aCTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTcgcg                                                          <  1:3713275/150‑1 (MQ=255)
                                                                                             cTCACGTAATCGGTGTTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGa                                                         >  2:502043/1‑150 (MQ=255)
                                                                                             cTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGa                                                         >  1:2938715/1‑150 (MQ=255)
                                                                                             cTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGa                                                         >  2:685968/1‑150 (MQ=255)
                                                                                               cACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACg                                                       <  1:2115867/150‑1 (MQ=255)
                                                                                                 cGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGtt                                                     <  1:4153718/150‑1 (MQ=255)
                                                                                                 cGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGtt                                                     <  2:1359162/150‑1 (MQ=255)
                                                                                                  gTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACAGTAATGGATATTCTTTAATAATATCTCGACGttt                                                    >  2:3347854/1‑150 (MQ=255)
                                                                                                    aaTCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCt                                                  <  1:4302146/150‑1 (MQ=255)
                                                                                                    aaTCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCt                                                  >  2:446101/1‑150 (MQ=255)
                                                                                                     aTCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCtt                                                 <  2:1548095/150‑1 (MQ=255)
                                                                                                       cGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTa                                               <  1:3479566/150‑1 (MQ=255)
                                                                                                        gttATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTaa                                              <  2:2174967/148‑2 (MQ=255)
                                                                                                        ggTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTAt                                              <  1:2385830/150‑1 (MQ=255)
                                                                                                        ggTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTAt                                              <  1:65710/150‑1 (MQ=255)
                                                                                                         gTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATg                                             <  2:429767/150‑1 (MQ=255)
                                                                                                         gTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATg                                             >  1:4086301/1‑150 (MQ=255)
                                                                                                          tATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGa                                            >  2:1179538/1‑150 (MQ=255)
                                                                                                             tGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataa                                         >  1:3203476/1‑150 (MQ=255)
                                                                                                              gTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaa                                        >  2:2004615/1‑150 (MQ=255)
                                                                                                                   tattCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAaataata                                   <  1:4371201/150‑1 (MQ=255)
                                                                                                                      tCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATc                                <  2:3512990/150‑1 (MQ=255)
                                                                                                                       cATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCa                               >  2:1001084/1‑150 (MQ=255)
                                                                                                                        aTTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCaa                              >  1:1368977/1‑150 (MQ=255)
                                                                                                                        aTTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCaa                              >  2:2573979/1‑150 (MQ=255)
                                                                                                                        aTTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCaa                              >  1:752484/1‑150 (MQ=255)
                                                                                                                              ttGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGag                        >  2:2356375/1‑149 (MQ=255)
                                                                                                                                gTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCTAATTGATaa                      >  2:3613858/1‑150 (MQ=255)
                                                                                                                                 tCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATaaa                     <  1:341087/150‑1 (MQ=255)
                                                                                                                                   cATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTATTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAAt                   <  2:4654693/150‑1 (MQ=255)
                                                                                                                                       ggaggaGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCaaa               >  2:1913934/1‑150 (MQ=255)
                                                                                                                                          ggagATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGAGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATg            <  2:4386978/150‑1 (MQ=255)
                                                                                                                                           gagaTTTAATCCCTTGCTTTTGTAAAACAATTGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATga           <  1:5164467/150‑1 (MQ=255)
                                                                                                                                            agaTTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATgag          >  2:1768916/1‑150 (MQ=255)
                                                                                                                                             gaTTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATgaga         >  2:2403905/1‑150 (MQ=255)
                                                                                                                                              aTTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGaa        >  2:4732810/1‑150 (MQ=255)
                                                                                                                                              aTTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGaa        >  2:4565689/1‑150 (MQ=255)
                                                                                                                                                 tAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGaaaaa     <  1:4471307/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAATAT  >  CP009273/3837028‑3837325

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