Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,755,723 G→T 100% E61* (GAA→TAA)  rhsA → Rhs family protein, putative polymorphic toxin; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,755,7230GT97.3% 173.5 / ‑0.7 75E61* (GAA→TAA) rhsARhs family protein, putative polymorphic toxin; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (25/48);  total (25/50)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-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.

GACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCG  >  CP009273/3755578‑3755769
                                                                                                                                                 |                                              
gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaa                                              >  2:2689319/1‑148 (MQ=21)
gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaa                                              >  1:987634/1‑148 (MQ=21)
gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaa                                              <  1:1750538/148‑1 (MQ=21)
gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCGGTTCGGGGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGACCTACCCGGTaaa                                              >  2:2689320/1‑148 (MQ=11)
gACGCAGTATGGAGGTAGCATTGTTCAGGGTTAAGCCGGGGTGCGCATTGGTGCCCCAACCGGCGTGGCCTGTTCGGTGTGACCCGGCGGAGTGACGTCAGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaa                                              >  2:3322536/1‑148 (MQ=11)
 aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaaa                                             <  1:1163987/148‑1 (MQ=21)
 aCGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTaaa                                             >  2:2438568/1‑148 (MQ=21)
    cAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCg                                          <  2:1778976/148‑1 (MQ=21)
    cAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCg                                          <  2:1779007/148‑1 (MQ=21)
       tATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACa                                       >  1:1763084/1‑148 (MQ=21)
        aTGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                      >  1:673213/1‑148 (MQ=21)
        aTGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACAt                                      >  1:673161/1‑148 (MQ=21)
          ggCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCg                                    <  2:2027288/148‑1 (MQ=21)
           gCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGc                                   <  1:1680491/148‑1 (MQ=21)
            cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGcc                                  <  2:2521061/148‑1 (MQ=21)
            cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGcc                                  <  1:836560/148‑1 (MQ=21)
            cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGcc                                  <  1:2528182/148‑1 (MQ=21)
            cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGcc                                  <  2:569457/148‑1 (MQ=21)
            cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGcc                                  <  1:862998/148‑1 (MQ=21)
             ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGccc                                 <  2:174323/148‑1 (MQ=21)
             ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGccc                                 <  2:190215/148‑1 (MQ=21)
             ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGccc                                 <  2:174563/148‑1 (MQ=21)
             ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGccc                                 <  1:2675803/148‑1 (MQ=34)
             ggTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGccc                                 <  1:2675936/148‑1 (MQ=34)
              gTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCATGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGccat                                >  2:3196383/1‑146 (MQ=19)
                aGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgc                              >  2:178162/1‑148 (MQ=21)
                aGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgc                              <  1:666207/148‑1 (MQ=21)
                  cATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCCTCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgccc                            >  1:2486687/1‑148 (MQ=21)
                   aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccg                           >  2:549555/1‑148 (MQ=21)
                   aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccg                           >  2:549605/1‑148 (MQ=21)
                   aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccg                           >  2:549682/1‑148 (MQ=21)
                    ttGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccgg                          >  1:1579960/1‑148 (MQ=21)
                    ttGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccgg                          <  2:1465301/148‑1 (MQ=21)
                     tGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggc                         <  2:777065/148‑1 (MQ=21)
                     tGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggc                         <  2:776895/148‑1 (MQ=21)
                      gTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggcc                        <  1:1374283/148‑1 (MQ=21)
                       ttCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggccc                       >  1:2298738/1‑148 (MQ=21)
                        tCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGGTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTgcccggcccg                      >  2:3007145/1‑148 (MQ=21)
                          aGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGGTCGGGGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCAGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCCCt                    >  1:736817/1‑148 (MQ=11)
                           gggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTg                   <  2:1399953/148‑1 (MQ=21)
                            ggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGc                  >  2:447374/1‑148 (MQ=21)
                            ggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGc                  <  1:2099613/148‑1 (MQ=21)
                            ggTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGc                  >  2:1912160/1‑148 (MQ=21)
                              ttCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                <  1:1171474/148‑1 (MQ=21)
                              ttCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                <  1:1173387/148‑1 (MQ=21)
                              ttCAGCCGGGGTGCGCATGGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                <  1:713870/148‑1 (MQ=21)
                                cAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGtt              <  2:3354064/148‑1 (MQ=21)
                                cAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGtt              >  2:2572999/1‑148 (MQ=21)
                                 aGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc             <  1:2636315/148‑1 (MQ=21)
                                 aGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc             <  2:1280274/148‑1 (MQ=21)
                                 aGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc             >  2:2087012/1‑148 (MQ=21)
                                   ccGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAt           >  1:1834661/1‑148 (MQ=21)
                                    cGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATc          <  2:456769/148‑1 (MQ=21)
                                    cGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATc          <  2:1335421/148‑1 (MQ=21)
                                    cGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATc          <  2:458108/148‑1 (MQ=21)
                                    cGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATc          <  1:398974/148‑1 (MQ=21)
                                      gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct        <  2:2378650/148‑1 (MQ=21)
                                      gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct        <  2:1109428/148‑1 (MQ=21)
                                      gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct        >  2:3461260/1‑148 (MQ=21)
                                      gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct        <  2:2380927/148‑1 (MQ=21)
                                      gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTCATCct        >  2:2987458/1‑148 (MQ=21)
                                        gTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct      <  1:2094105/148‑1 (MQ=21)
                                        gTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct      <  1:2094076/148‑1 (MQ=21)
                                        gTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCctct      <  1:2094080/148‑1 (MQ=21)
                                         tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCTCGCTGCCGTTCATCCtctc     <  2:673161/148‑1 (MQ=18)
                                         tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCTCGCTGCCGTTCATCCtctc     <  2:673213/148‑1 (MQ=18)
                                         tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCtctc     <  1:957072/148‑1 (MQ=18)
                                         tGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCtctc     <  2:466499/148‑1 (MQ=18)
                                          gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAGAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc    <  2:617075/148‑1 (MQ=17)
                                          gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc    <  2:180735/148‑1 (MQ=17)
                                          gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc    <  2:254730/148‑1 (MQ=17)
                                          gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc    <  2:617015/148‑1 (MQ=17)
                                          gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGGTAAACCGACCTCGCCCTGCCCGGCCCGCGGCCGTTCATCCTCGcc    >  1:3118661/1‑148 (MQ=11)
                                           cgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTccc   <  1:1663522/148‑1 (MQ=17)
                                            gcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCg  <  1:1653467/148‑1 (MQ=14)
                                            gcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTTAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCg  <  1:1653464/148‑1 (MQ=14)
                                                                                                                                                 |                                              
GACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCG  >  CP009273/3755578‑3755769

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