Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,591,139 C→A 100% D66Y (GAT→TAT)  dnaT ← DNA biosynthesis protein (primosomal protein I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,591,1390CA98.8% 589.0 / ‑1.7 186D66Y (GAT→TAT) dnaTDNA biosynthesis protein (primosomal protein I)
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (91/92);  total (91/95)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 187 total aligned reads displayed.

TCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTAC  >  CP009273/4590994‑4591286
                                                                                                                                                 |                                                                                                                                                   
tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGc                                                                                                                                                   <  2:3064611/148‑1 (MQ=255)
 ctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCt                                                                                                                                                  >  1:3626424/1‑148 (MQ=255)
 ctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCt                                                                                                                                                  >  2:1918690/1‑148 (MQ=255)
 ctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCt                                                                                                                                                  >  1:3138879/1‑148 (MQ=255)
 ctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCt                                                                                                                                                  >  1:2041621/1‑148 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                                <  1:1636790/148‑1 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                                <  1:1637768/148‑1 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                                >  1:720947/1‑148 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                                >  1:2361556/1‑148 (MQ=255)
    ttAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTcc                                                                                                                                               <  1:2412582/148‑1 (MQ=255)
       gcgcCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCgg                                                                                                                                            <  2:865266/148‑1 (MQ=255)
           cACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCCTGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGAcgc                                                                                                                                        >  2:3208393/1‑148 (MQ=255)
             cGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACgcgc                                                                                                                                      >  2:2913173/1‑148 (MQ=255)
               ccccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCa                                                                                                                                    >  2:3544385/1‑148 (MQ=255)
                cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAg                                                                                                                                   <  1:560566/148‑1 (MQ=255)
                      gcgcCGCCAGGCGGATAAAATCTGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttc                                                                                                                             >  1:300250/1‑148 (MQ=255)
                      gcgcCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttc                                                                                                                             >  1:300272/1‑148 (MQ=255)
                         ccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttcttc                                                                                                                          >  2:1978478/1‑148 (MQ=255)
                          cgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTcc                                                                                                                         <  2:3377716/148‑1 (MQ=255)
                            ccAGGCGGATAAACTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAg                                                                                                                       <  2:727547/148‑1 (MQ=255)
                              aGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcg                                                                                                                     >  1:304235/1‑148 (MQ=255)
                                gCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGcc                                                                                                                   >  2:59004/1‑148 (MQ=255)
                                   gATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagc                                                                                                                <  1:3097216/148‑1 (MQ=255)
                                    aTAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagca                                                                                                               <  2:1858994/148‑1 (MQ=255)
                                    aTAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagca                                                                                                               <  2:1347699/148‑1 (MQ=255)
                                      aaaaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGt                                                                                                             >  2:2050740/1‑148 (MQ=255)
                                             gCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                      <  2:1858486/148‑1 (MQ=255)
                                               aTCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTGCAGCGCCAGCAGTTCAGCcagg                                                                                                    <  1:2359127/148‑1 (MQ=255)
                                                tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggc                                                                                                   >  1:1054706/1‑148 (MQ=255)
                                                 cGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcg                                                                                                  >  2:2228224/1‑148 (MQ=255)
                                                     ttGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgca                                                                                              <  1:3672362/148‑1 (MQ=255)
                                                       gCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcagg                                                                                            <  2:1232155/148‑1 (MQ=255)
                                                           gTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcgtg                                                                                        <  1:3470639/148‑1 (MQ=255)
                                                           gTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcgtg                                                                                        <  1:2245503/148‑1 (MQ=255)
                                                              cGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACg                                                                                     <  2:3411688/148‑1 (MQ=255)
                                                                gATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGc                                                                                   >  1:3498269/1‑148 (MQ=255)
                                                                 aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCa                                                                                  <  2:1588792/148‑1 (MQ=255)
                                                                 aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCa                                                                                  >  2:3055161/1‑148 (MQ=255)
                                                                   aCATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATa                                                                                >  1:830576/1‑148 (MQ=255)
                                                                     aTGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaa                                                                              >  2:649201/1‑148 (MQ=255)
                                                                         cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGc                                                                          >  2:2800580/1‑148 (MQ=255)
                                                                           aaTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCg                                                                        >  2:1789125/1‑148 (MQ=255)
                                                                           aaTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCg                                                                        <  2:2284029/148‑1 (MQ=255)
                                                                             tttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCggg                                                                      >  2:2622984/1‑148 (MQ=255)
                                                                              ttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgggg                                                                     >  1:1660282/1‑148 (MQ=255)
                                                                                    aTGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGt                                                               <  2:2860666/148‑1 (MQ=255)
                                                                                     tGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGtt                                                              <  2:105255/148‑1 (MQ=255)
                                                                                           cAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCaaa                                                        <  1:3017343/148‑1 (MQ=255)
                                                                                            aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                       <  2:376431/148‑1 (MQ=255)
                                                                                            aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                       <  2:376443/148‑1 (MQ=255)
                                                                                            aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                       <  2:912290/148‑1 (MQ=255)
                                                                                            aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                       <  2:912307/148‑1 (MQ=255)
                                                                                              cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                     >  1:1711072/1‑148 (MQ=255)
                                                                                              cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATGGTTAGCAAATAc                                                     >  2:2967924/1‑148 (MQ=255)
                                                                                              cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACAAAGCTTCCCGGACGCGCCAGCTTTTCTTCCACCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                     >  1:1317795/1‑148 (MQ=255)
                                                                                               aaCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACg                                                    <  1:168265/148‑1 (MQ=255)
                                                                                               aaCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACg                                                    <  1:2087605/148‑1 (MQ=255)
                                                                                                aCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACgg                                                   >  2:3080617/1‑148 (MQ=255)
                                                                                                 cGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGc                                                  >  1:847372/1‑148 (MQ=255)
                                                                                                   ggAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGATCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTTCAGGCGTGACGGCATAAAACGCCGGGGCATTTTTAGCAAAAACGGCaa                                                >  1:2237327/1‑148 (MQ=255)
                                                                                                   ggAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCaa                                                >  1:309204/1‑148 (MQ=255)
                                                                                                    gAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCCTAAAACGCCGGGGCATTGTTAGCAAATACGGCaaa                                               >  1:2950765/1‑147 (MQ=255)
                                                                                                      gCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACca                                             >  2:1833223/1‑148 (MQ=255)
                                                                                                        ggCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcaca                                           <  1:2744582/148‑1 (MQ=255)
                                                                                                          cTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACAcc                                         >  2:2855155/1‑148 (MQ=255)
                                                                                                          cTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACAcc                                         >  2:2855350/1‑148 (MQ=255)
                                                                                                            tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                       <  1:2481736/148‑1 (MQ=255)
                                                                                                            tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                       >  1:2538671/1‑148 (MQ=255)
                                                                                                             gCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGcc                                      <  1:612307/148‑1 (MQ=255)
                                                                                                             gCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGcc                                      <  2:1756431/148‑1 (MQ=255)
                                                                                                               ggTTCCTTATAGTGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCTGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGAATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCtt                                    <  1:3201121/148‑1 (MQ=255)
                                                                                                               ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCCGCGCCAGCAGTTCAGCCAGGCGTGCAGGGGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCtt                                    <  1:1001381/148‑1 (MQ=255)
                                                                                                               ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCtt                                    <  1:1001600/148‑1 (MQ=255)
                                                                                                                gTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                   <  1:3136243/148‑1 (MQ=255)
                                                                                                                gTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                   <  2:1054706/148‑1 (MQ=255)
                                                                                                                gTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                   <  1:1554399/148‑1 (MQ=255)
                                                                                                                gTTCCTGATAGAGTGGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAACACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                   <  1:448287/148‑1 (MQ=255)
                                                                                                                  tCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAg                                 <  2:1149705/148‑1 (MQ=255)
                                                                                                                   ccTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGc                                >  1:3269173/1‑148 (MQ=255)
                                                                                                                   ccTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGc                                >  1:237476/1‑148 (MQ=255)
                                                                                                                   ccTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGc                                <  2:3647809/148‑1 (MQ=255)
                                                                                                                    cTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCt                               <  1:3055161/148‑1 (MQ=255)
                                                                                                                     tGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtt                              >  2:1568370/1‑148 (MQ=255)
                                                                                                                      gATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                             <  2:3138879/148‑1 (MQ=255)
                                                                                                                      gATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                             <  2:3626424/148‑1 (MQ=255)
                                                                                                                       aTAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtttt                            >  1:942602/1‑148 (MQ=255)
                                                                                                                       aTAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtttt                            >  1:588736/1‑148 (MQ=255)
                                                                                                                        tAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTg                           >  1:1483498/1‑148 (MQ=255)
                                                                                                                        tAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTg                           >  1:1193716/1‑148 (MQ=255)
                                                                                                                                  cGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGt                 >  2:1912556/1‑148 (MQ=255)
                                                                                                                                  cGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGt                 >  1:489513/1‑148 (MQ=255)
                                                                                                                                    tCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAACACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGttt               <  2:2363716/148‑1 (MQ=255)
                                                                                                                                    tCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGttt               <  1:2729872/148‑1 (MQ=255)
                                                                                                                                     ccAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg              >  2:3027626/1‑148 (MQ=255)
                                                                                                                                     ccAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg              >  1:3568081/1‑148 (MQ=255)
                                                                                                                                        gagCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg           >  1:1846380/1‑148 (MQ=255)
                                                                                                                                        gagCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg           >  1:2986302/1‑148 (MQ=255)
                                                                                                                                           cGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATc        >  1:3458622/1‑148 (MQ=255)
                                                                                                                                             aCATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt      >  2:2112160/1‑148 (MQ=255)
                                                                                                                                              cATAGCGTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg     <  2:121306/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTAC  >  CP009273/4590994‑4591286

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