Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,210,753 C→A 100% E555* (GAA→TAA)  arpA ← ankyrin repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,210,7530CA98.7% 730.5 / 0.4 228E555* (GAA→TAA) arpAankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (111/114);  total (112/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 228 total aligned reads displayed.

GTATTAAACAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTCCTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACATGT  >  CP009273/4210608‑4210897
                                                                                                                                                 |                                                                                                                                                
 tATTAAACAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCCTAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTg                                                                                                                                               >  1:326846/1‑148 (MQ=255)
  aTTAAACAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTgg                                                                                                                                              >  1:1269161/1‑148 (MQ=255)
   ttAAACAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTggg                                                                                                                                             >  2:1221234/1‑148 (MQ=255)
   ttAAACAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTggg                                                                                                                                             >  2:531154/1‑148 (MQ=255)
        cAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAgg                                                                                                                                        <  2:2425802/148‑1 (MQ=255)
             tGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTggg                                                                                                                                   <  2:3011326/148‑1 (MQ=255)
              gggTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGc                                                                                                                                  <  1:2746012/148‑1 (MQ=255)
              gggTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGc                                                                                                                                  >  1:2536212/1‑148 (MQ=255)
               ggTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCa                                                                                                                                 >  1:1727907/1‑148 (MQ=255)
               ggTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCa                                                                                                                                 >  2:2627267/1‑148 (MQ=255)
                gTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCaa                                                                                                                                <  2:1921038/148‑1 (MQ=255)
                gTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCaa                                                                                                                                <  2:3171366/148‑1 (MQ=255)
                   aTGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGc                                                                                                                             >  2:852756/1‑148 (MQ=255)
                   aTGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGc                                                                                                                             <  2:1099387/148‑1 (MQ=255)
                    tGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCt                                                                                                                            >  1:791898/1‑148 (MQ=255)
                       gTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCa                                                                                                                         <  1:741371/148‑1 (MQ=255)
                          aaaaaTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGa                                                                                                                      >  1:1860470/1‑148 (MQ=255)
                           aaaaTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGaa                                                                                                                     <  1:2922772/148‑1 (MQ=255)
                             aaTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATa                                                                                                                   <  1:462493/148‑1 (MQ=255)
                             aaTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATa                                                                                                                   >  1:3277734/1‑148 (MQ=255)
                             aaTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATa                                                                                                                   <  1:2650058/148‑1 (MQ=255)
                             aaTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATa                                                                                                                   <  2:76962/148‑1 (MQ=255)
                             aaTAGTATTAATAACGTTAATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATa                                                                                                                   >  2:362110/1‑148 (MQ=255)
                              aTAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATaa                                                                                                                  >  1:3150242/1‑148 (MQ=255)
                               tAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAAc                                                                                                                 <  2:16907/148‑1 (MQ=255)
                               tAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAAc                                                                                                                 <  1:723708/148‑1 (MQ=255)
                                 gTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCATGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACtt                                                                                                               <  2:1568182/148‑1 (MQ=255)
                                 gTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACtt                                                                                                               <  2:2768277/148‑1 (MQ=255)
                                  tATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACttt                                                                                                              <  1:2198279/148‑1 (MQ=255)
                                    ttAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCa                                                                                                            >  1:1870516/1‑148 (MQ=255)
                                     taataaCGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCAc                                                                                                           >  2:2774060/1‑148 (MQ=255)
                                      aataaCGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACg                                                                                                          <  1:2227281/148‑1 (MQ=255)
                                       ataaCGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGa                                                                                                         <  1:3388780/148‑1 (MQ=255)
                                           cGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATc                                                                                                     <  2:451662/148‑1 (MQ=255)
                                           cGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATc                                                                                                     >  2:2022318/1‑148 (MQ=255)
                                            gTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCg                                                                                                    >  2:1056251/1‑148 (MQ=255)
                                            gTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCg                                                                                                    <  2:1860470/148‑1 (MQ=255)
                                             ttCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGc                                                                                                   >  2:438704/1‑148 (MQ=255)
                                             ttCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGc                                                                                                   >  2:3018907/1‑148 (MQ=255)
                                              tCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCt                                                                                                  <  2:3463064/148‑1 (MQ=255)
                                              tCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCt                                                                                                  >  2:384577/1‑148 (MQ=255)
                                                   tgtCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTc                                                                                             <  2:469184/148‑1 (MQ=255)
                                                      ccGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGATATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGATTAACTTTCACGATATCGCTATGTCCAt                                                                                          <  1:2501131/148‑1 (MQ=255)
                                                       cgcgCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCAtt                                                                                         >  2:1413858/1‑148 (MQ=255)
                                                       cgcgCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCAtt                                                                                         <  2:333669/148‑1 (MQ=255)
                                                       cgcgCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCAtt                                                                                         >  1:2059160/1‑148 (MQ=255)
                                                        gcgcTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCAttt                                                                                        <  1:362110/148‑1 (MQ=255)
                                                             gCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCAt                                                                                   <  2:2897557/148‑1 (MQ=255)
                                                                 ggCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGcc                                                                               <  2:1580435/148‑1 (MQ=255)
                                                                    cATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGa                                                                            >  1:851268/1‑148 (MQ=255)
                                                                     aTAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGat                                                                           >  1:2558462/1‑148 (MQ=255)
                                                                      tAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGata                                                                          >  1:1397389/1‑148 (MQ=255)
                                                                       aaacaaacCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGatat                                                                         >  1:1149243/1‑148 (MQ=255)
                                                                       aaacaaacCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGatat                                                                         <  2:1377093/148‑1 (MQ=255)
                                                                         acaaacCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAg                                                                       >  2:2554816/1‑148 (MQ=255)
                                                                           aaacCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTc                                                                     <  1:1565289/148‑1 (MQ=255)
                                                                           aaacCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTc                                                                     <  1:195166/148‑1 (MQ=255)
                                                                                   gtCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCaa                                                             >  2:1014355/1‑148 (MQ=255)
                                                                                    tCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAAc                                                            <  1:1065215/148‑1 (MQ=255)
                                                                                    tCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAAc                                                            >  2:1401335/1‑148 (MQ=255)
                                                                                    tCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAAc                                                            >  2:1633235/1‑148 (MQ=255)
                                                                                      gcgcgcAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCa                                                          <  1:1534635/148‑1 (MQ=255)
                                                                                       cgcgcAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCAt                                                         <  2:1269161/148‑1 (MQ=255)
                                                                                        gcgcAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATa                                                        <  2:3205729/148‑1 (MQ=255)
                                                                                         cgcAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATaa                                                       <  1:2022318/148‑1 (MQ=255)
                                                                                          gcAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATaaa                                                      <  1:3516770/148‑1 (MQ=255)
                                                                                              gACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATc                                                  >  2:1973212/1‑148 (MQ=255)
                                                                                              gACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATc                                                  <  1:2698940/148‑1 (MQ=255)
                                                                                                 ttttAGCGGTCAGAAGATCTAGAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCtt                                               >  2:2945773/1‑148 (MQ=255)
                                                                                                   ttAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTg                                             <  2:2249375/148‑1 (MQ=255)
                                                                                                       cGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATCTAGTCCTGGGCAACCATAAAAATCCTTTGCCtt                                         >  2:2119968/1‑148 (MQ=255)
                                                                                                       cGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCtt                                         <  2:3076298/148‑1 (MQ=255)
                                                                                                        ggTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTc                                        <  2:3173240/148‑1 (MQ=255)
                                                                                                        ggTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTc                                        <  2:2391842/148‑1 (MQ=255)
                                                                                                        ggTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTc                                        <  1:9115/148‑1 (MQ=255)
                                                                                                             gaagaTCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGgag                                   >  1:3082265/1‑148 (MQ=255)
                                                                                                              aataTCTACAATATAGGAATCTGAAATGTTTATTTACTGTGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAATATCCTTTGCCTTCAGgaga                                  <  2:184009/145‑1 (MQ=255)
                                                                                                              aagaTCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGgaga                                  <  1:1640382/148‑1 (MQ=255)
                                                                                                              aagaTCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGgaga                                  >  1:2960225/1‑148 (MQ=255)
                                                                                                                  tCTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGa                              <  2:3442624/148‑1 (MQ=255)
                                                                                                                   cTACAATATCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAg                             <  2:796047/148‑1 (MQ=255)
                                                                                                                       aaTATCGGAAGCTGAAATGTTAATTTCCTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTAcc                         <  1:2867649/148‑1 (MQ=255)
                                                                                                                        atatCGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCt                        >  1:2403115/1‑148 (MQ=255)
                                                                                                                            cGGAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTc                    >  1:3479270/1‑148 (MQ=255)
                                                                                                                              gAAGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCtt                  <  1:2878028/148‑1 (MQ=255)
                                                                                                                               aaGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCttt                 >  1:2273097/1‑148 (MQ=255)
                                                                                                                               aaGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCttt                 >  1:2273491/1‑148 (MQ=255)
                                                                                                                                aGCTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTa                >  1:540117/1‑148 (MQ=255)
                                                                                                                                  cTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGt              <  1:2683877/148‑1 (MQ=255)
                                                                                                                                  cTGAAATGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGt              >  2:2489425/1‑148 (MQ=255)
                                                                                                                                      aaTGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCaaa          <  2:3146224/148‑1 (MQ=255)
                                                                                                                                      aaTGTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCaaa          <  1:1183930/148‑1 (MQ=255)
                                                                                                                                         gTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGa       >  2:3255630/1‑148 (MQ=255)
                                                                                                                                         gTTAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGa       <  2:2296657/148‑1 (MQ=255)
                                                                                                                                          ttAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGAc      >  2:2814562/1‑148 (MQ=255)
                                                                                                                                           tAATTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACa     <  2:2511385/148‑1 (MQ=255)
                                                                                                                                            aaTTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACAt    <  2:2665766/148‑1 (MQ=255)
                                                                                                                                             aTTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACAtg   <  1:2774060/148‑1 (MQ=255)
                                                                                                                                             aTTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACAtg   >  2:1802634/1‑148 (MQ=255)
                                                                                                                                             aTTTACTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACAtg   <  2:2705200/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTATTAAACAGAGTGGGTAATGCGTTAAAAATAGTATTAATAACGTTCATGTGTCCGCGCTGCATGGCCATAAACAAACCCGTGTCGCGCGCAAGACTTTTAGCGGTCAGAAGATCTACAATATCGGAAGCTGAAATGTTAATTTCCTGGGCTAGGCTGGGCAATGCTTCAAGAATAACTTTCACGATATCGCTATGTCCATTTTGCATCGCCAGATATAGTCCTGGGCAACCATAAAAATCCTTTGCCTTCAGGAGATCGAGTACCTGTTCTTTAGTCAAATGACATGT  >  CP009273/4210608‑4210897

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