Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,618,693 G→T 100% E74* (GAA→TAA)  ppx → exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,618,6930GT95.7% 696.1 / 20.2 233E74* (GAA→TAA) ppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base T (94/129);  total (99/134)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 233 total aligned reads displayed.

GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAG  >  CP009273/2618548‑2618837
                                                                                                                                                 |                                                                                                                                                
ggTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTaa                                                                                                                                                >  2:1609951/1‑148 (MQ=255)
 gTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAAc                                                                                                                                               <  1:516210/148‑1 (MQ=255)
      aGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTa                                                                                                                                          >  2:171860/1‑148 (MQ=255)
      aGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTa                                                                                                                                          >  1:574290/1‑148 (MQ=255)
      aGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTa                                                                                                                                          >  1:1586029/1‑148 (MQ=255)
      aGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGAAGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTa                                                                                                                                          <  1:1035104/148‑1 (MQ=255)
       gCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTAc                                                                                                                                         >  1:2510660/1‑148 (MQ=255)
           gtgtGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAgg                                                                                                                                     >  1:2628129/1‑148 (MQ=255)
            tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAggg                                                                                                                                    <  1:3507115/148‑1 (MQ=255)
            tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAggg                                                                                                                                    <  2:1485736/148‑1 (MQ=255)
             gtgGTAGATGGTGCCATGCAGATTATTGGCCGACTGAAACAGCGGCTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGTCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGt                                                                                                                                   >  1:1855929/1‑148 (MQ=255)
              tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtt                                                                                                                                  <  2:3242766/148‑1 (MQ=255)
                gTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtttt                                                                                                                                >  1:1066049/1‑148 (MQ=255)
                   gATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtctc                                                                                                                             <  2:408050/148‑1 (MQ=255)
                     tGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCt                                                                                                                           >  2:2619845/1‑148 (MQ=255)
                            aTGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCg                                                                                                                    >  1:3145625/1‑148 (MQ=255)
                             tGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGt                                                                                                                   <  1:1819359/148‑1 (MQ=255)
                             tGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGt                                                                                                                   >  1:3311659/1‑148 (MQ=255)
                              gCAGATTAATGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGAGAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                                                                                  >  2:65482/1‑148 (MQ=255)
                               cAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCt                                                                                                                 <  1:3299386/148‑1 (MQ=255)
                                  attattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGta                                                                                                              <  1:1618287/148‑1 (MQ=255)
                                   ttattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                                                                                             <  1:334998/148‑1 (MQ=255)
                                    tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                                                            <  2:111009/148‑1 (MQ=255)
                                    tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                                                            >  1:766592/1‑148 (MQ=255)
                                    tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                                                            <  2:600366/148‑1 (MQ=255)
                                     attGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                                                                                           >  2:1632220/1‑148 (MQ=255)
                                      ttGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                          >  1:3106624/1‑148 (MQ=255)
                                        ggCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTg                                                                                                        <  2:804185/148‑1 (MQ=255)
                                         gccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                                                                                       <  2:1848750/148‑1 (MQ=255)
                                           cgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTa                                                                                                     <  1:566885/148‑1 (MQ=255)
                                                gAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACGAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAt                                                                                                <  1:3135286/148‑1 (MQ=255)
                                                               tCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcg                                                                                 <  2:1549298/148‑1 (MQ=255)
                                                               tCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcg                                                                                 <  2:103972/148‑1 (MQ=255)
                                                                cTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcgc                                                                                <  2:1570112/148‑1 (MQ=255)
                                                                  ggCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTg                                                                              <  2:1183321/148‑1 (MQ=255)
                                                                       aCGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc                                                                         <  1:2169612/148‑1 (MQ=255)
                                                                             tGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa                                                                   <  1:2994940/148‑1 (MQ=255)
                                                                               gggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACt                                                                 >  1:1515668/1‑148 (MQ=255)
                                                                                 gCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACttt                                                               <  2:2911282/148‑1 (MQ=255)
                                                                                 gCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACttt                                                               >  2:1537916/1‑148 (MQ=255)
                                                                                  ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTc                                                              <  2:3357902/148‑1 (MQ=255)
                                                                                  ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTc                                                              <  2:3357948/148‑1 (MQ=255)
                                                                                   cAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCt                                                             <  2:2781012/148‑1 (MQ=255)
                                                                                    aGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTg                                                            >  1:2582365/1‑148 (MQ=255)
                                                                                     gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGa                                                           >  2:1975837/1‑148 (MQ=255)
                                                                                     gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGa                                                           <  2:2195796/148‑1 (MQ=255)
                                                                                         atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                       <  2:994204/148‑1 (MQ=255)
                                                                                         atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                       >  1:2019048/1‑148 (MQ=255)
                                                                                           atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                     >  2:1264725/1‑148 (MQ=255)
                                                                                           atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                     >  2:3446142/1‑148 (MQ=255)
                                                                                              ttGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCgg                                                  >  1:1108652/1‑148 (MQ=255)
                                                                                              ttGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCgg                                                  >  1:2483316/1‑148 (MQ=255)
                                                                                                 aGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaaa                                               <  1:1248912/148‑1 (MQ=255)
                                                                                                  gTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaaaa                                              <  2:3477144/148‑1 (MQ=255)
                                                                                                   tGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAg                                             <  2:2628137/148‑1 (MQ=255)
                                                                                                     aaGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGt                                           <  2:2238038/148‑1 (MQ=255)
                                                                                                       gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTCCAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCa                                         >  2:2168218/1‑148 (MQ=255)
                                                                                                       gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCa                                         >  2:902455/1‑148 (MQ=255)
                                                                                                         ggCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCAtt                                       <  2:2385572/148‑1 (MQ=255)
                                                                                                           cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                     <  1:545033/148‑1 (MQ=255)
                                                                                                           cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                     <  2:91269/148‑1 (MQ=255)
                                                                                                           cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                     >  2:2004/1‑148 (MQ=255)
                                                                                                            aaTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTccc                                    <  2:3515275/148‑1 (MQ=255)
                                                                                                            aaTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTccc                                    <  2:2510660/148‑1 (MQ=255)
                                                                                                            aaTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTccc                                    <  1:859668/148‑1 (MQ=255)
                                                                                                              tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                  >  1:1750703/1‑148 (MQ=255)
                                                                                                              tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                  <  1:1971332/148‑1 (MQ=255)
                                                                                                              tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                  >  1:3445905/1‑148 (MQ=255)
                                                                                                              tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACCCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                  >  1:2724308/1‑148 (MQ=255)
                                                                                                               gACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAc                                 <  2:2280490/148‑1 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                                >  1:618227/1‑148 (MQ=255)
                                                                                                                   cgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAACAGGTCATTCCCTACCCGa                             >  2:138089/1‑148 (MQ=255)
                                                                                                                   cgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGa                             <  1:73584/148‑1 (MQ=255)
                                                                                                                   cgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGa                             >  1:2126881/1‑148 (MQ=255)
                                                                                                                    gcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGAt                            <  2:55203/148‑1 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCGATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                        >  1:2097140/1‑148 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                        <  1:660148/148‑1 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                        <  2:4297/148‑1 (MQ=255)
                                                                                                                         ttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaaa                       <  2:3285115/148‑1 (MQ=255)
                                                                                                                          tAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAat                      >  1:3437810/1‑148 (MQ=255)
                                                                                                                          tAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAat                      >  1:1668749/1‑148 (MQ=255)
                                                                                                                           aaaCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAatt                     <  2:2458201/148‑1 (MQ=255)
                                                                                                                             aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattat                   >  1:1300197/1‑148 (MQ=255)
                                                                                                                             aCTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattat                   >  2:2364778/1‑148 (MQ=255)
                                                                                                                              ctgtctgtCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                  >  2:1687479/1‑148 (MQ=255)
                                                                                                                              ctgtctgtCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                  >  1:2362619/1‑148 (MQ=255)
                                                                                                                               tgtctgtcGCTGTTTGCCTAACGGCTACTAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATAAGCTGCGTCAGGCGCTGATCGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCAGATTGAAATTAttt                 <  1:916383/148‑1 (MQ=255)
                                                                                                                                 tctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc               >  2:2908880/1‑148 (MQ=255)
                                                                                                                                 tctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc               <  2:449289/148‑1 (MQ=255)
                                                                                                                                   tgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCgg             <  1:1537916/148‑1 (MQ=255)
                                                                                                                                   tgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCgg             <  1:712422/148‑1 (MQ=255)
                                                                                                                                    gtcGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGt            <  2:3311659/148‑1 (MQ=255)
                                                                                                                                     tcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTa           >  2:2664968/1‑148 (MQ=255)
                                                                                                                                        cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATg        >  2:708130/1‑148 (MQ=255)
                                                                                                                                        cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATg        >  2:2492294/1‑148 (MQ=255)
                                                                                                                                        cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGAGTGAAATTATTTCCGGTAATg        >  2:1916378/1‑148 (MQ=255)
                                                                                                                                           tttGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaag     >  2:181527/1‑148 (MQ=255)
                                                                                                                                             tGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaagaa   <  2:3133542/148‑1 (MQ=255)
                                                                                                                                             tGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaagaa   <  2:964329/148‑1 (MQ=255)
                                                                                                                                              gCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGaagaag  <  1:848108/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAG  >  CP009273/2618548‑2618837

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: ‑