Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,644,531 G→A 100% intergenic (+509/‑120) arsC → / → yhiS arsenate reductase/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,644,5310GA98.4% 1251.7 / 5.3 383intergenic (+509/‑120)arsC/yhiSarsenate reductase/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (164/213);  total (166/217)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 383 total aligned reads displayed.

AATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGT  >  CP009273/3644386‑3644674
                                                                                                                                                 |                                                                                                                                               
 aTCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATaa                                                                                                                                              <  1:1909206/148‑1 (MQ=255)
  tCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAAt                                                                                                                                             >  2:2604039/1‑148 (MQ=255)
     cTTCACACACTATAAGAGAACATAATAACTATATAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATtat                                                                                                                                          <  2:3561904/148‑1 (MQ=255)
       tCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTatat                                                                                                                                        <  1:145539/148‑1 (MQ=255)
           acacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAgg                                                                                                                                    >  2:1536603/1‑148 (MQ=255)
            cacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGa                                                                                                                                   >  2:3535562/1‑148 (MQ=255)
            cacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGa                                                                                                                                   <  2:631887/148‑1 (MQ=255)
                aTCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAATGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGat                                                                                                                               <  1:2815258/148‑1 (MQ=255)
                aTCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGat                                                                                                                               <  1:2668449/148‑1 (MQ=255)
                  cAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGatat                                                                                                                             >  1:1252151/1‑148 (MQ=255)
                  cAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGatat                                                                                                                             <  2:2389437/148‑1 (MQ=255)
                     agaACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAAc                                                                                                                          >  1:2333523/1‑148 (MQ=255)
                      gaACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACg                                                                                                                         >  2:3699407/1‑148 (MQ=255)
                        aCATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCa                                                                                                                       <  1:10918/148‑1 (MQ=255)
                        aCATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCa                                                                                                                       <  1:3282015/148‑1 (MQ=255)
                        aCATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCa                                                                                                                       >  2:1074775/1‑148 (MQ=255)
                         cATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCAt                                                                                                                      <  1:3116374/148‑1 (MQ=255)
                         cATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCAt                                                                                                                      >  1:2060068/1‑148 (MQ=255)
                         cATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCAt                                                                                                                      >  2:2426575/1‑148 (MQ=255)
                         cATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCAt                                                                                                                      <  1:1327228/148‑1 (MQ=255)
                           taataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTa                                                                                                                    >  1:2387844/1‑148 (MQ=255)
                            aataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTAc                                                                                                                   >  1:480898/1‑148 (MQ=255)
                            aataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTAc                                                                                                                   <  2:3653826/148‑1 (MQ=255)
                            aataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTAc                                                                                                                   >  1:2978397/1‑148 (MQ=255)
                             ataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACa                                                                                                                  >  1:3447220/1‑148 (MQ=255)
                             ataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACa                                                                                                                  >  1:3388292/1‑148 (MQ=255)
                              taaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACAt                                                                                                                 <  1:10083/148‑1 (MQ=255)
                               aaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACAtt                                                                                                                <  1:763364/148‑1 (MQ=255)
                                 cTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTa                                                                                                              <  2:3388631/148‑1 (MQ=255)
                                  tttttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTAt                                                                                                             <  1:135809/148‑1 (MQ=255)
                                   ttttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTAtt                                                                                                            <  2:2993174/148‑1 (MQ=255)
                                        aaaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATtgtgtg                                                                                                       <  1:1145083/148‑1 (MQ=255)
                                         aaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTgtgtgt                                                                                                      >  1:197433/1‑148 (MQ=255)
                                          aTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTa                                                                                                     >  1:2747590/1‑148 (MQ=255)
                                            gAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTaga                                                                                                   >  1:2446800/1‑148 (MQ=255)
                                            gAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTaga                                                                                                   <  1:1718672/148‑1 (MQ=255)
                                              gcgcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTagaag                                                                                                 <  1:1790286/148‑1 (MQ=255)
                                               cgcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTagaaga                                                                                                >  1:317236/1‑148 (MQ=255)
                                               cgcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTagaaga                                                                                                >  2:560438/1‑148 (MQ=255)
                                                  aTTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGaaaa                                                                                             <  1:1825506/148‑1 (MQ=255)
                                                   ttttGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAAt                                                                                            >  2:917774/1‑148 (MQ=255)
                                                    tttGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATc                                                                                           >  1:3356909/1‑148 (MQ=255)
                                                    tttGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATc                                                                                           >  2:1818889/1‑148 (MQ=255)
                                                        aTAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCAtt                                                                                       <  2:1598435/148‑1 (MQ=255)
                                                         tAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTc                                                                                      >  2:449580/1‑148 (MQ=255)
                                                         tAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTc                                                                                      >  2:2956351/1‑148 (MQ=255)
                                                         tAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTc                                                                                      >  2:2670126/1‑148 (MQ=255)
                                                         tAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTc                                                                                      >  2:1419971/1‑148 (MQ=255)
                                                            ttgttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTg                                                                                   >  2:3752394/1‑148 (MQ=255)
                                                             tgttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGt                                                                                  >  1:2595232/1‑148 (MQ=255)
                                                              gttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGtt                                                                                 <  1:836477/148‑1 (MQ=255)
                                                                      tttCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTAt                                                                         >  2:3011390/1‑148 (MQ=255)
                                                                      tttCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTAt                                                                         <  1:2923460/148‑1 (MQ=255)
                                                                       ttCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATc                                                                        <  2:1935268/148‑1 (MQ=255)
                                                                        tCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATcc                                                                       >  1:601939/1‑148 (MQ=255)
                                                                          cTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCaa                                                                     <  2:1336221/148‑1 (MQ=255)
                                                                           tttACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAAt                                                                    <  1:1199872/148‑1 (MQ=255)
                                                                            ttACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAAtt                                                                   >  2:2111537/1‑148 (MQ=255)
                                                                             tACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAAttt                                                                  >  2:2875897/1‑148 (MQ=255)
                                                                               cGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTaa                                                                <  2:2841864/148‑1 (MQ=255)
                                                                                 aTACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAAc                                                              >  2:2739438/1‑148 (MQ=255)
                                                                                    cTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATc                                                           <  2:2012319/148‑1 (MQ=255)
                                                                                      aTCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCtt                                                         <  1:223101/148‑1 (MQ=255)
                                                                                        cACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAg                                                       >  2:3225046/1‑148 (MQ=255)
                                                                                          cTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCa                                                     >  1:3047831/1‑148 (MQ=255)
                                                                                          cTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCa                                                     <  1:2109449/148‑1 (MQ=255)
                                                                                             cccTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaa                                                  >  1:872393/1‑148 (MQ=255)
                                                                                              ccTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaaa                                                 <  2:3320755/148‑1 (MQ=255)
                                                                                                 ttttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACTGCTTAGCATAAATCa                                              >  2:2250419/1‑148 (MQ=255)
                                                                                                 ttttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACa                                              <  1:3529617/148‑1 (MQ=255)
                                                                                                 ttttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACa                                              <  2:2152181/148‑1 (MQ=255)
                                                                                                  tttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACaa                                             <  1:2899725/148‑1 (MQ=255)
                                                                                                   ttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACaaa                                            <  1:3402940/148‑1 (MQ=255)
                                                                                                    tttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAat                                           >  1:529744/1‑148 (MQ=255)
                                                                                                     ttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAatg                                          <  2:421300/148‑1 (MQ=255)
                                                                                                      tACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAatga                                         <  2:884032/148‑1 (MQ=255)
                                                                                                      tACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAatga                                         <  1:984115/148‑1 (MQ=255)
                                                                                                          cacaAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGaa                                     >  2:2504973/1‑148 (MQ=255)
                                                                                                           acaAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAAt                                    <  1:2187557/148‑1 (MQ=255)
                                                                                                            caAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATa                                   <  2:3302583/148‑1 (MQ=255)
                                                                                                               atatTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAgg                                <  2:2936447/148‑1 (MQ=255)
                                                                                                                tatTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGa                               >  2:745596/1‑148 (MQ=255)
                                                                                                                tatTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGa                               <  2:197433/148‑1 (MQ=255)
                                                                                                                  ttACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAAt                             <  2:3482834/148‑1 (MQ=255)
                                                                                                                   tACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAAtt                            >  1:3579833/1‑148 (MQ=255)
                                                                                                                    aCCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAAttt                           >  1:16994/1‑148 (MQ=255)
                                                                                                                      ctgttggatataaaatggatataaaaTATAATTTTATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATGATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTa                         <  2:3244672/144‑1 (MQ=255)
                                                                                                                        ggTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATg                       >  2:815907/1‑148 (MQ=255)
                                                                                                                        ggTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATg                       <  1:3265641/148‑1 (MQ=255)
                                                                                                                            gagaTAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCaa                   <  2:1034653/148‑1 (MQ=255)
                                                                                                                            gagaTAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCaa                   <  2:2333523/148‑1 (MQ=255)
                                                                                                                                tAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGa               <  1:1419971/148‑1 (MQ=255)
                                                                                                                                 aaaaGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGAc              >  1:416780/1‑148 (MQ=255)
                                                                                                                                   aaGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACtt            >  1:694440/1‑148 (MQ=255)
                                                                                                                                     gggAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTa          >  2:3063424/1‑148 (MQ=255)
                                                                                                                                     gggAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTa          >  2:3111286/1‑148 (MQ=255)
                                                                                                                                      ggAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTAc         >  2:671866/1‑148 (MQ=255)
                                                                                                                                           tCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGACTTTTATGTCAATTGACTTTACCCCAg    >  1:730229/1‑148 (MQ=255)
                                                                                                                                           tCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAg    >  2:1231425/1‑148 (MQ=255)
                                                                                                                                             aaaGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGt  <  2:1123076/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGT  >  CP009273/3644386‑3644674

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: