Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,644,531 G→A 91.6% 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,5310GA91.6% 766.6 / 50.9 249intergenic (+509/‑120)arsC/yhiSarsenate reductase/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base A (108/120);  total (118/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 252 total aligned reads displayed.

AAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAA  >  CP009273/3644384‑3644678
                                                                                                                                                   |                                                                                                                                                   
    tCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAAt                                                                                                                                                 >  1:1824160/1‑148 (MQ=255)
     cACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAAtt                                                                                                                                                <  2:465937/148‑1 (MQ=255)
     cACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAAtt                                                                                                                                                >  1:784651/1‑148 (MQ=255)
     cACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAAtt                                                                                                                                                <  1:1550958/148‑1 (MQ=255)
      aCTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATta                                                                                                                                               <  1:808990/148‑1 (MQ=255)
        ttCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATtata                                                                                                                                             >  1:2122817/1‑148 (MQ=255)
        ttCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATtata                                                                                                                                             >  2:1739052/1‑148 (MQ=255)
        ttCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATtata                                                                                                                                             <  2:3081611/148‑1 (MQ=255)
         tCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTatat                                                                                                                                            <  2:2730504/148‑1 (MQ=255)
         tCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTatat                                                                                                                                            <  1:256004/148‑1 (MQ=255)
          cacacaCTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATc                                                                                                                                           <  2:922655/148‑1 (MQ=255)
           acacacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCa                                                                                                                                          >  2:636130/1‑148 (MQ=255)
           acacacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCa                                                                                                                                          >  2:1609657/1‑148 (MQ=255)
            cacacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAg                                                                                                                                         >  1:589175/1‑148 (MQ=255)
             acacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATAAAAGATAATTATATCAgg                                                                                                                                        >  2:1558962/1‑148 (MQ=255)
              cacTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGa                                                                                                                                       <  1:920383/148‑1 (MQ=255)
               acTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGaa                                                                                                                                      <  1:1316071/148‑1 (MQ=255)
                cTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAg                                                                                                                                     <  1:3177384/148‑1 (MQ=255)
                cTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAg                                                                                                                                     <  1:1535387/148‑1 (MQ=255)
                cTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAg                                                                                                                                     <  1:1358892/148‑1 (MQ=255)
                 tATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGa                                                                                                                                    <  1:188499/148‑1 (MQ=255)
                 tATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGa                                                                                                                                    >  1:302446/1‑148 (MQ=255)
                     agagAACCTAATAACTTTCTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGAtata                                                                                                                                <  1:2128665/148‑1 (MQ=255)
                      gagaACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATaa                                                                                                                               >  2:2107705/1‑148 (MQ=255)
                        gaACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACg                                                                                                                             >  1:2671041/1‑148 (MQ=255)
                        gaACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACg                                                                                                                             >  2:2915740/1‑148 (MQ=255)
                            ataataACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCAtt                                                                                                                         >  2:1311584/1‑148 (MQ=255)
                              aataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTaa                                                                                                                       >  2:3048607/1‑147 (MQ=255)
                              aataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTAc                                                                                                                       >  1:736424/1‑148 (MQ=255)
                               ataaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACa                                                                                                                      >  2:3293884/1‑148 (MQ=255)
                                taaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACAt                                                                                                                     <  1:1902179/148‑1 (MQ=255)
                                 aaCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAATAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACAtt                                                                                                                    >  2:1671465/1‑148 (MQ=255)
                                  aCTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGTGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACAttt                                                                                                                   <  1:2830746/148‑1 (MQ=255)
                                    tttttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTAt                                                                                                                 >  2:2163248/1‑148 (MQ=255)
                                    tttttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTAt                                                                                                                 >  2:2163237/1‑148 (MQ=255)
                                      tttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATtg                                                                                                               <  2:2999149/148‑1 (MQ=255)
                                       ttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATtgt                                                                                                              >  1:3312244/1‑148 (MQ=255)
                                       ttAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATtgt                                                                                                              >  2:2545426/1‑148 (MQ=255)
                                           aaTGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCGTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTgtgtgt                                                                                                          <  1:3235037/148‑1 (MQ=255)
                                              gAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTaga                                                                                                       <  2:698606/148‑1 (MQ=255)
                                                gcgcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTagaag                                                                                                     >  2:2953927/1‑148 (MQ=255)
                                                 ctcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTagaaga                                                                                                    <  2:2070094/146‑1 (MQ=255)
                                                 cgcATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTagaaga                                                                                                    >  1:905034/1‑148 (MQ=255)
                                                   cATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGaaa                                                                                                  >  2:524550/1‑148 (MQ=255)
                                                      tttGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATc                                                                                               >  2:172691/1‑148 (MQ=255)
                                                      tttGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATc                                                                                               >  2:1843927/1‑148 (MQ=255)
                                                        tGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCa                                                                                             <  1:643470/148‑1 (MQ=255)
                                                        tGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCa                                                                                             <  2:1376678/148‑1 (MQ=255)
                                                           tAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTc                                                                                          <  1:2697286/148‑1 (MQ=255)
                                                              ttgttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTg                                                                                       <  1:1609657/148‑1 (MQ=255)
                                                              ttgttAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTg                                                                                       <  1:2944802/148‑1 (MQ=255)
                                                                    aCGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTaa                                                                                 >  2:2446802/1‑148 (MQ=255)
                                                                       aTTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTa                                                                              >  1:462641/1‑148 (MQ=255)
                                                                        tttCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTAt                                                                             <  1:1403085/148‑1 (MQ=255)
                                                                            cTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCaa                                                                         <  1:3299495/148‑1 (MQ=255)
                                                                             tttACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAAt                                                                        >  2:65433/1‑148 (MQ=255)
                                                                                 cGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTaa                                                                    <  1:2974375/148‑1 (MQ=255)
                                                                                    tACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACa                                                                 >  2:1064611/1‑148 (MQ=255)
                                                                                        aTCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCtt                                                             <  1:1859113/148‑1 (MQ=255)
                                                                                        aTCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCtt                                                             >  2:2636258/1‑148 (MQ=255)
                                                                                        aTCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCtt                                                             >  2:2022803/1‑148 (MQ=255)
                                                                                           aCTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGc                                                          <  1:2065892/148‑1 (MQ=255)
                                                                                            cTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCa                                                         >  2:1438463/1‑148 (MQ=255)
                                                                                             tACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCAt                                                        <  1:52664/148‑1 (MQ=255)
                                                                                             tACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCAt                                                        <  1:3034721/148‑1 (MQ=255)
                                                                                              aCCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATa                                                       >  2:2958524/1‑148 (MQ=255)
                                                                                               cccTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaa                                                      <  1:1771570/148‑1 (MQ=255)
                                                                                                ccTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaaa                                                     >  1:1703825/1‑148 (MQ=255)
                                                                                                ccTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaaa                                                     >  1:3126959/1‑148 (MQ=255)
                                                                                                 cTTTTTTTACACACAAATATTACCGTGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaaaa                                                    <  2:2000575/148‑1 (MQ=255)
                                                                                                 cTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATaaaa                                                    >  1:1001025/1‑148 (MQ=255)
                                                                                                  tttttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAAc                                                   >  2:2386685/1‑148 (MQ=255)
                                                                                                    tttttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACaa                                                 >  2:1480078/1‑148 (MQ=255)
                                                                                                       ttACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAatg                                              <  1:2291300/148‑1 (MQ=255)
                                                                                                            cacaAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGaa                                         <  2:2164972/148‑1 (MQ=255)
                                                                                                             acaAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAAt                                        <  1:1036703/148‑1 (MQ=255)
                                                                                                              caAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATa                                       >  2:2738634/1‑148 (MQ=255)
                                                                                                              caAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATa                                       >  2:1504731/1‑148 (MQ=255)
                                                                                                                aaTATTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAg                                     <  1:1311584/148‑1 (MQ=255)
                                                                                                                  tatTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGa                                   <  1:2590350/148‑1 (MQ=255)
                                                                                                                  tatTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGa                                   >  1:500927/1‑148 (MQ=255)
                                                                                                                  tatTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGa                                   <  1:1236708/148‑1 (MQ=255)
                                                                                                                   atTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGaa                                  <  2:1222232/148‑1 (MQ=255)
                                                                                                                   atTACCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGaa                                  <  2:2870207/148‑1 (MQ=255)
                                                                                                                     tACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAAtt                                <  1:636130/148‑1 (MQ=255)
                                                                                                                      aCCGGGTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAAttt                               <  1:102274/148‑1 (MQ=255)
                                                                                                                          ggTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATg                           >  2:1541935/1‑148 (MQ=255)
                                                                                                                           gTGGAGATAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGt                          <  1:2002853/148‑1 (MQ=255)
                                                                                                                              gagaTAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCaa                       >  2:2819776/1‑148 (MQ=255)
                                                                                                                               agaTAAAAGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAAt                      <  2:890219/148‑1 (MQ=255)
                                                                                                                                   aaaaGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGACTTTTATGTCAATTGAc                  >  1:382980/1‑148 (MQ=255)
                                                                                                                                    aaaGGGAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACt                 <  2:784651/148‑1 (MQ=255)
                                                                                                                                         gAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTAcc            <  1:2260250/148‑1 (MQ=255)
                                                                                                                                         gAAATCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTAcc            <  1:833073/148‑1 (MQ=255)
                                                                                                                                          aaaTCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTAccc           >  1:2458752/1‑148 (MQ=255)
                                                                                                                                             tCAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAg        >  2:2197847/1‑148 (MQ=255)
                                                                                                                                              cAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAgg       <  1:2915740/148‑1 (MQ=255)
                                                                                                                                              cAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAgg       >  2:812948/1‑148 (MQ=255)
                                                                                                                                              cAAAGATAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAgg       <  1:172691/148‑1 (MQ=255)
                                                                                                                                                  ggTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGtata   >  2:2375425/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAATCACTTCACACACTATCAGAGAACATAATAACTTTTTAAAATGAGCGCATTTTGATAGTTGTTAACGATTTCCTTTACGATACTATCACTACCCTTTTTTTACACACAAATATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAA  >  CP009273/3644384‑3644678

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: