Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,869,780 G→A 100% intergenic (‑92/‑160) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,869,7800GA96.5% 747.1 / 12.9 230intergenic (‑92/‑160)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base A (110/112);  total (112/118)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 230 total aligned reads displayed.

GTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGT  >  CP009273/2869635‑2869924
                                                                                                                                                 |                                                                                                                                                
gtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAag                                                                                                                                                <  1:1410881/148‑1 (MQ=255)
       cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTcc                                                                                                                                         >  2:420818/1‑148 (MQ=255)
       cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTcc                                                                                                                                         <  2:3568093/148‑1 (MQ=255)
         ttCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCtt                                                                                                                                       >  1:3790997/1‑148 (MQ=255)
          tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTa                                                                                                                                      <  1:414811/148‑1 (MQ=255)
                tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGaa                                                                                                                                <  2:1011967/148‑1 (MQ=255)
                   aGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATaa                                                                                                                             <  2:1790463/148‑1 (MQ=255)
                   aGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATaa                                                                                                                             <  2:3235361/148‑1 (MQ=255)
                     ttGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACg                                                                                                                           >  1:2325308/1‑148 (MQ=255)
                      tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                          <  2:368445/148‑1 (MQ=255)
                      tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                          <  1:3281810/148‑1 (MQ=255)
                      tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                          <  2:3153591/148‑1 (MQ=255)
                      tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                          <  1:2504172/148‑1 (MQ=255)
                       gCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGAt                                                                                                                         <  1:535311/148‑1 (MQ=255)
                          gCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTg                                                                                                                      <  1:3798865/148‑1 (MQ=255)
                             ggTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                                                                                   >  1:3223971/1‑148 (MQ=255)
                              gTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                                                                  >  2:2148408/1‑148 (MQ=255)
                                   tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTa                                                                                                             >  2:1557842/1‑148 (MQ=255)
                                    aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                                                                            >  1:220187/1‑148 (MQ=255)
                                     aGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAat                                                                                                           <  2:2975932/148‑1 (MQ=255)
                                      ttcttATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAata                                                                                                          <  2:2962661/144‑1 (MQ=255)
                                            tttGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaa                                                                                                    >  1:920445/1‑148 (MQ=255)
                                            tttGATCCATAACCGTTCCTTTGCAAGACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaa                                                                                                    <  2:337300/148‑1 (MQ=255)
                                                 tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGc                                                                                               >  2:3175943/1‑148 (MQ=255)
                                                 tCCATAACAGTTCCTTTGCACTACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGc                                                                                               <  2:1257048/148‑1 (MQ=255)
                                                  ccATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCt                                                                                              <  2:301881/148‑1 (MQ=255)
                                                    aTAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTa                                                                                            <  2:466718/148‑1 (MQ=255)
                                                      aaCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                                                                          >  2:2998773/1‑148 (MQ=255)
                                                       aCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                                                                         <  2:1306358/148‑1 (MQ=255)
                                                       aCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                                                                         >  2:3714502/1‑148 (MQ=255)
                                                        ccGATCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATCAAATTACCAATTAGAATAATTAGTAACTAAACGGAAAAACGAATTGGCAAAGCTAATATCAAAAAGAGCTTAAgcc                                                                                        >  1:2819759/1‑146 (MQ=255)
                                                         cGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                                                                       <  2:2916427/148‑1 (MQ=255)
                                                          gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                                                                      <  2:2319041/148‑1 (MQ=255)
                                                          gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                                                                      >  1:1567593/1‑148 (MQ=255)
                                                            tCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAcag                                                                                    >  2:2808578/1‑146 (MQ=255)
                                                             ccTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                                                                   >  2:493502/1‑148 (MQ=255)
                                                              cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGa                                                                                  <  2:3555619/148‑1 (MQ=255)
                                                                ttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAtt                                                                                >  2:579045/1‑148 (MQ=255)
                                                                 tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAttt                                                                               <  2:972344/148‑1 (MQ=255)
                                                                   cAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                             <  1:3137772/148‑1 (MQ=255)
                                                                      tACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGc                                                                          >  1:1169558/1‑148 (MQ=255)
                                                                       aCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCg                                                                         <  2:2053496/148‑1 (MQ=255)
                                                                       aCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCg                                                                         >  1:2408213/1‑148 (MQ=255)
                                                                        ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCACTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGt                                                                        <  2:2497973/148‑1 (MQ=255)
                                                                        ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGt                                                                        <  2:2030216/148‑1 (MQ=255)
                                                                        ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGt                                                                        <  2:568391/148‑1 (MQ=255)
                                                                         cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGtt                                                                       >  2:3051740/1‑148 (MQ=255)
                                                                          gCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGttt                                                                      <  2:657310/148‑1 (MQ=255)
                                                                             aTTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAgg                                                                   >  1:902529/1‑148 (MQ=255)
                                                                             aTTTTCTTGCCATCAGATGTTTCGACTATACGGAGCGTAAGAGAACGAATTAAATAACCAATTAGAATAAGTATTTCCTTAACGTAATAACTATTTGGCAAATCTAATATCAAAAAGTGCTTAAGTCACCCGATTACTAGCGTGTAgg                                                                   >  1:2883753/1‑148 (MQ=255)
                                                                              ttttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGa                                                                  <  1:1909022/148‑1 (MQ=255)
                                                                              ttttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGa                                                                  >  1:3417949/1‑148 (MQ=255)
                                                                               tttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGaa                                                                 >  2:1928523/1‑148 (MQ=255)
                                                                               tttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGaa                                                                 >  1:3310472/1‑148 (MQ=255)
                                                                               tttCTTGCCATAAGATGATTCGACTATAGGGAGCGTAAGAGAAAGAATGAAATTACCAATTAGAATAATTAGTTCCTTAACGTAATAACGATTTGGCAAATCTAATATCAAAAAGTGCTTAAGGCACCGTATTTCGTGCGTTTAGGaa                                                                 >  2:3814798/1‑148 (MQ=255)
                                                                                  cTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAt                                                              <  1:660296/148‑1 (MQ=255)
                                                                                  cTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAt                                                              <  2:118848/148‑1 (MQ=255)
                                                                                   ttGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAtt                                                             <  1:1648120/148‑1 (MQ=255)
                                                                                    tGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAttt                                                            >  2:2400888/1‑148 (MQ=255)
                                                                                      ccATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGa                                                          <  1:303374/148‑1 (MQ=255)
                                                                                       cATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaa                                                         <  2:2641603/148‑1 (MQ=255)
                                                                                        aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa                                                        <  2:882214/148‑1 (MQ=255)
                                                                                         tCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCATAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                       <  2:2681088/148‑1 (MQ=255)
                                                                                          cAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAtt                                                      >  2:1674955/1‑148 (MQ=255)
                                                                                           aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTg                                                     <  2:447594/148‑1 (MQ=255)
                                                                                              tGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGttt                                                  >  1:3723690/1‑148 (MQ=255)
                                                                                               gTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtttt                                                 >  1:1461230/1‑148 (MQ=255)
                                                                                               gTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtttt                                                 >  1:1834312/1‑148 (MQ=255)
                                                                                                tttCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTa                                                <  2:2289898/148‑1 (MQ=255)
                                                                                                 ttCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAg                                               <  1:3526357/148‑1 (MQ=255)
                                                                                                     aCTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCa                                           >  1:2980111/1‑148 (MQ=255)
                                                                                                      cTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAg                                          <  2:1681404/148‑1 (MQ=255)
                                                                                                        ataGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                        >  2:144082/1‑148 (MQ=255)
                                                                                                         taGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCgg                                       >  1:1052531/1‑148 (MQ=255)
                                                                                                           gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                     >  2:1857007/1‑148 (MQ=255)
                                                                                                           gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                     >  2:1711129/1‑148 (MQ=255)
                                                                                                             gAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGt                                   >  1:132875/1‑148 (MQ=255)
                                                                                                              aGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                                  >  2:3117566/1‑148 (MQ=255)
                                                                                                                cGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTc                                <  2:2249730/148‑1 (MQ=255)
                                                                                                                cGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTc                                <  2:3310472/148‑1 (MQ=255)
                                                                                                                  tAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCAt                              <  1:934882/148‑1 (MQ=255)
                                                                                                                    agagAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATAc                            >  2:1297358/1‑148 (MQ=255)
                                                                                                                    agagAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATAc                            >  1:735251/1‑148 (MQ=255)
                                                                                                                     gagaACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACt                           >  2:1956425/1‑148 (MQ=255)
                                                                                                                       gaACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTAt                         <  2:1051975/148‑1 (MQ=255)
                                                                                                                         aCGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATgg                       <  1:2140564/148‑1 (MQ=255)
                                                                                                                         aCGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATgg                       <  2:3089565/148‑1 (MQ=255)
                                                                                                                            aaTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCgg                    >  2:305966/1‑148 (MQ=255)
                                                                                                                             aTGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGt                   >  1:3819745/1‑148 (MQ=255)
                                                                                                                              tGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTa                  <  2:2192923/148‑1 (MQ=255)
                                                                                                                              tGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACAGTATTTCGGGCGTTTAGGACGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATAATATGTCGGTa                  >  1:391593/1‑148 (MQ=255)
                                                                                                                               gAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaa                 >  1:1131313/1‑148 (MQ=255)
                                                                                                                                aaaTTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaa                <  1:3555946/148‑1 (MQ=255)
                                                                                                                                   ttACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaa             >  2:1521335/1‑148 (MQ=255)
                                                                                                                                   aaaCCAATTATAATAAGTAGATCCTTAACGGAAAAAAGATTTGTCAAAGCAAATAACAAAAAGTGCTTAAGGCAACGTATTTCGGGCGTTTAGTAAGATTTGAAATTGTTTTAGCTCAGCGGCAGTTTCATACTATGGCGGTaaacaa             >  2:1604355/3‑148 (MQ=255)
                                                                                                                                    tACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaaa            >  2:1875045/1‑148 (MQ=255)
                                                                                                                                     aCCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAt           <  1:2706041/148‑1 (MQ=255)
                                                                                                                                     aCCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAt           <  2:936556/148‑1 (MQ=255)
                                                                                                                                          ttAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa      <  1:711827/148‑1 (MQ=255)
                                                                                                                                              aaTAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGt  >  2:3502182/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGT  >  CP009273/2869635‑2869924

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: