Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,872,520 G→T 100% P183P (CCG→CCT yeaP → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,872,5200GT99.6% 747.5 / ‑4.4 227P183P (CCG→CCTyeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (109/117);  total (110/117)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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 227 total aligned reads displayed.

GGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTT  >  CP009273/1872375‑1872663
                                                                                                                                                 |                                                                                                                                               
ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTaa                                                                                                                                               <  1:2199305/148‑1 (MQ=255)
ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTaa                                                                                                                                               <  2:853451/148‑1 (MQ=255)
   gTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcg                                                                                                                                            <  2:1410213/148‑1 (MQ=255)
    tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcgg                                                                                                                                           >  2:616088/1‑148 (MQ=255)
            cgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAt                                                                                                                                   <  2:3097775/148‑1 (MQ=255)
            cgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAt                                                                                                                                   <  2:1615418/148‑1 (MQ=255)
             ggAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAtt                                                                                                                                  <  1:3167435/148‑1 (MQ=255)
              gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttt                                                                                                                                 <  2:70244/148‑1 (MQ=255)
              gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttt                                                                                                                                 >  2:1053583/1‑148 (MQ=255)
               aaCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAtttt                                                                                                                                >  1:1462648/1‑148 (MQ=255)
                aCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                               >  2:1019302/1‑148 (MQ=255)
                aCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                               <  2:1122236/148‑1 (MQ=255)
                 cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTg                                                                                                                              >  2:351187/1‑148 (MQ=255)
                  aGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGa                                                                                                                             <  1:2557276/148‑1 (MQ=255)
                   ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAANCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGaa                                                                                                                            >  1:1808226/1‑148 (MQ=255)
                   ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaa                                                                                                                            <  1:2087163/148‑1 (MQ=255)
                    gTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaa                                                                                                                           >  1:3507276/1‑148 (MQ=255)
                     ttttACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGACTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaaa                                                                                                                          >  2:1321895/1‑148 (MQ=255)
                               agtCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACGGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGaca                                                                                                                <  1:3548942/146‑1 (MQ=255)
                                ttCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACCGCTGCGCGAAGCCATTGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGacac                                                                                                               <  2:2467790/148‑1 (MQ=255)
                                ttCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGacac                                                                                                               <  2:3068541/148‑1 (MQ=255)
                                 tCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACt                                                                                                              >  2:1797884/1‑148 (MQ=255)
                                  cGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTg                                                                                                             >  1:1604881/1‑148 (MQ=255)
                                  cGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTg                                                                                                             <  1:1011306/148‑1 (MQ=255)
                                    cgggAGTGTTTGCACAATATTTTGAATTAGTTGCAGTTTTTTAACAGCTGCGGGATGAAATTGATGGTCTTATTGCTCATTGTTATACCGACTCTTTTACCGGGCTACGTATTCGGCGGGCGATTTTTGAAAATCTGACGACACTGtt                                                                                                           <  1:466749/146‑1 (MQ=255)
                                      tgACTGATTGCACATTATATTCAAAAAGAGGCACTGGTTGAACAGCTGGGCGTTGCCATTGCTGAGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTTCCTAATCGGCGGGCGATTTTTGATAATCTGACGACAATGtttt                                                                                                         <  2:122487/147‑1 (MQ=255)
                                      ggACTGATTGAACTATTTTTTCAAAAAGAGGCAATGGTTGAACAGCTGCGCGAAGCCAATGGTGCGATGTTTGGGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTTAAAATCTGAAGACACTGtttt                                                                                                         <  2:1511515/148‑1 (MQ=255)
                                          tGATTGCACAATATATTCAAAAAGAGTCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCTTTAACCGGGCTACCTAATCGGAGGGCGATTTTTGAAAATCTGACGACACTTTGTTACCt                                                                                                     >  2:2799176/1‑148 (MQ=255)
                                           gATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                    >  2:703693/1‑148 (MQ=255)
                                            aTTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCg                                                                                                   <  2:2982198/148‑1 (MQ=255)
                                            aTTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCg                                                                                                   <  2:1765957/148‑1 (MQ=255)
                                              tGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAAACGGGCTACCTAATCGGCGGGCGATTTTTGAAAATATGACGACACTGTTTTCCCTCgac                                                                                                 >  1:1061200/1‑146 (MQ=255)
                                                  caATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGc                                                                                             <  2:1650825/148‑1 (MQ=255)
                                                     tatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATc                                                                                          <  1:3238048/148‑1 (MQ=255)
                                                      atatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                         >  2:786654/1‑148 (MQ=255)
                                                       tatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCtt                                                                                        >  2:364671/1‑148 (MQ=255)
                                                         ttCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTaa                                                                                      <  2:868375/148‑1 (MQ=255)
                                                         ttCAAAAAGAGGCACTGGTGGATCAGCGGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCTTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGTTTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTaa                                                                                      <  2:2276783/148‑1 (MQ=255)
                                                             aaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                  >  2:688803/1‑148 (MQ=255)
                                                             aaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                  >  2:688874/1‑148 (MQ=255)
                                                             aaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                  >  1:862658/1‑148 (MQ=255)
                                                              aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                 >  2:1460605/1‑148 (MQ=255)
                                                              aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                 >  2:1460619/1‑148 (MQ=255)
                                                              aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                 >  1:2395923/1‑148 (MQ=255)
                                                              aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                 >  2:2969205/1‑148 (MQ=255)
                                                              aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                 >  2:2929946/1‑148 (MQ=255)
                                                                  agGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAt                                                                             >  1:421228/1‑148 (MQ=255)
                                                                   ggCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATa                                                                            >  2:106077/1‑148 (MQ=255)
                                                                           ttGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                    >  1:310053/1‑148 (MQ=255)
                                                                            tGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgcg                                                                   <  2:1082518/148‑1 (MQ=255)
                                                                            tGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgcg                                                                   >  2:721459/1‑148 (MQ=255)
                                                                             gAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGt                                                                  <  1:3451025/148‑1 (MQ=255)
                                                                              aaCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGtt                                                                 >  1:338982/1‑148 (MQ=255)
                                                                               aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                >  1:938900/1‑148 (MQ=255)
                                                                               aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                <  1:1768310/148‑1 (MQ=255)
                                                                               aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                <  2:1057407/148‑1 (MQ=255)
                                                                               aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTAGACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                <  2:1006743/148‑1 (MQ=255)
                                                                                 aGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTAt                                                              <  2:3128186/148‑1 (MQ=255)
                                                                                 aGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACAGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTAt                                                              >  2:3408554/1‑148 (MQ=255)
                                                                                   cTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCg                                                            <  1:2072163/148‑1 (MQ=255)
                                                                                     gcgcGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGAt                                                          <  2:2038305/148‑1 (MQ=255)
                                                                                      cgcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATc                                                         <  1:2863201/148‑1 (MQ=255)
                                                                                       gcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCt                                                        <  2:2455846/148‑1 (MQ=255)
                                                                                       gcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCt                                                        <  2:515190/148‑1 (MQ=255)
                                                                                          aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGa                                                     >  1:139058/1‑148 (MQ=255)
                                                                                          aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGa                                                     <  2:1569793/148‑1 (MQ=255)
                                                                                          aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGa                                                     <  1:3340618/148‑1 (MQ=255)
                                                                                           aGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGAt                                                    >  1:1984303/1‑148 (MQ=255)
                                                                                           aGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGAt                                                    >  2:2194761/1‑148 (MQ=255)
                                                                                             ccAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAGGATAATGATCGCGTTTATCGATCTGGATaa                                                  <  2:1020334/148‑1 (MQ=255)
                                                                                               aaTGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACt                                                >  1:1654851/1‑148 (MQ=255)
                                                                                                   ctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCaa                                            <  2:3449129/148‑1 (MQ=255)
                                                                                                     gcgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAt                                          <  2:945654/148‑1 (MQ=255)
                                                                                                     gcgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAt                                          <  2:2728116/148‑1 (MQ=255)
                                                                                                     gcgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACAATGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAt                                          >  1:2302913/1‑148 (MQ=255)
                                                                                                          gATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATc                                     >  1:2884039/1‑148 (MQ=255)
                                                                                                               cgcAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGa                                <  2:1871791/148‑1 (MQ=255)
                                                                                                                gcAATCGTATACCGACTCGTTAACCGGTCTACCTAATCGGCGGTCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAAATTCAAATTAATCAATGAt                               <  2:3025637/148‑1 (MQ=255)
                                                                                                                    tCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                           >  1:95781/1‑148 (MQ=255)
                                                                                                                    tCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGAGCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                           >  1:1741571/1‑148 (MQ=255)
                                                                                                                     cGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACGTCAAATTAATCAATGATCGttt                          >  2:1972912/1‑148 (MQ=255)
                                                                                                                        ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTgg                       >  2:1806694/1‑148 (MQ=255)
                                                                                                                          aCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                     <  2:457968/148‑1 (MQ=255)
                                                                                                                           ccGTCTGGTTATCCGGGCTACCTTATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCTTAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGTTCa                    <  2:2949405/148‑1 (MQ=255)
                                                                                                                             gACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCata                  <  1:1019302/148‑1 (MQ=255)
                                                                                                                               cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataat                <  1:721296/148‑1 (MQ=255)
                                                                                                                               cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataat                >  2:676841/1‑148 (MQ=255)
                                                                                                                               cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataat                >  1:2270281/1‑148 (MQ=255)
                                                                                                                                tCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataata               >  2:3429837/1‑148 (MQ=255)
                                                                                                                                  gTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt             <  1:29696/148‑1 (MQ=255)
                                                                                                                                  gTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt             <  1:1834991/148‑1 (MQ=255)
                                                                                                                                   ttAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTg            <  1:707591/148‑1 (MQ=255)
                                                                                                                                    tAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTgg           <  2:2408406/148‑1 (MQ=255)
                                                                                                                                      aCCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCg         >  1:1163338/1‑148 (MQ=255)
                                                                                                                                       ccGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGa        <  1:616088/148‑1 (MQ=255)
                                                                                                                                       ccGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGa        >  2:1168344/1‑148 (MQ=255)
                                                                                                                                        cGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGAt       <  1:702774/148‑1 (MQ=255)
                                                                                                                                          ggCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCt     >  1:683888/1‑148 (MQ=255)
                                                                                                                                           gCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTGTTGGTGATAATAGTGGCGATCTg    >  2:2774613/1‑148 (MQ=255)
                                                                                                                                             tACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGtt  <  2:1749088/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTT  >  CP009273/1872375‑1872663

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: