Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,110,391 C→A 17.2% R73L (CGC→CTC)  wcaL ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,110,3910CA17.2% 399.1 / 49.8 140R73L (CGC→CTC) wcaLputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (59/56);  new base A (6/18);  total (65/75)
Fisher's exact test for biased strand distribution p-value = 2.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 143 total aligned reads displayed.

ACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAA  >  CP009273/2110245‑2110536
                                                                                                                                                  |                                                                                                                                                 
acacATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGgcg                                                                                                                                                  <  1:1902418/148‑1 (MQ=255)
 cacaTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGATTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGgagc                                                                                                                                                 <  2:700351/148‑1 (MQ=255)
  acaTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCa                                                                                                                                                <  1:2198762/148‑1 (MQ=255)
   caTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAg                                                                                                                                               >  2:1266858/1‑148 (MQ=255)
   caTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAg                                                                                                                                               >  1:2603040/1‑148 (MQ=255)
     tCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGtt                                                                                                                                             >  1:3187772/1‑148 (MQ=255)
      cGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGttt                                                                                                                                            >  1:2723501/1‑148 (MQ=255)
        gcgcGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCg                                                                                                                                          <  2:1135406/148‑1 (MQ=255)
           cgAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCa                                                                                                                                       >  1:1415258/1‑148 (MQ=255)
            gAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCAc                                                                                                                                      <  2:1841936/148‑1 (MQ=255)
             aaaCGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACt                                                                                                                                     >  1:2897684/1‑148 (MQ=255)
              aaCGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGATGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACtt                                                                                                                                    <  2:2406079/148‑1 (MQ=255)
                 ggTGTTGCGACCTGGGCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCTCCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGAGCAAGGTCTGTCTGGATAGGTGGCGCAGTTTCGCCACTTTGc                                                                                                                                 <  2:144312/148‑1 (MQ=255)
                 ggTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGc                                                                                                                                 <  2:1296568/148‑1 (MQ=255)
                 ggTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGc                                                                                                                                 <  2:2373497/148‑1 (MQ=255)
                 ggTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGc                                                                                                                                 <  2:144332/148‑1 (MQ=255)
                     ttGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGt                                                                                                                             <  2:1279541/148‑1 (MQ=255)
                     ttGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGt                                                                                                                             >  2:2701375/1‑148 (MQ=255)
                            cTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCg                                                                                                                      >  1:2154939/1‑148 (MQ=255)
                             tGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt                                                                                                                     <  2:2813001/148‑1 (MQ=255)
                             tGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt                                                                                                                     <  1:2718090/148‑1 (MQ=255)
                                 cGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTg                                                                                                                 <  1:1400957/148‑1 (MQ=255)
                                    aaaTGGCAGACAAAATCAGGGTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTaa                                                                                                              <  2:2783576/148‑1 (MQ=255)
                                      aTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAcc                                                                                                            >  1:2222609/1‑148 (MQ=255)
                                         gCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGc                                                                                                         >  1:2074134/1‑148 (MQ=255)
                                         gCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGc                                                                                                         <  2:287929/148‑1 (MQ=255)
                                            gACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCggg                                                                                                      <  2:1483845/148‑1 (MQ=255)
                                             aCAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGt                                                                                                     >  1:2290305/1‑148 (MQ=255)
                                               aaaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTc                                                                                                   <  1:1942417/148‑1 (MQ=255)
                                               aaaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTc                                                                                                   <  1:1942415/148‑1 (MQ=255)
                                               aaaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTc                                                                                                   <  1:459079/148‑1 (MQ=255)
                                               aaaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTc                                                                                                   <  1:2433436/148‑1 (MQ=255)
                                               aaaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTc                                                                                                   <  1:1352020/148‑1 (MQ=255)
                                                 aaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGTGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTg                                                                                                 <  1:721160/148‑1 (MQ=255)
                                                  aTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTgg                                                                                                >  2:1612635/1‑148 (MQ=255)
                                                   tCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgc                                                                                               <  2:2381430/148‑1 (MQ=255)
                                                   tCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgc                                                                                               >  1:1035858/1‑148 (MQ=255)
                                                        ttCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCaa                                                                                          <  1:3150636/148‑1 (MQ=255)
                                                        ttCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCaa                                                                                          <  1:1509090/148‑1 (MQ=255)
                                                            gcgACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTGACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGtt                                                                                      <  2:2319178/148‑1 (MQ=255)
                                                             cgACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTg                                                                                     <  1:1266858/148‑1 (MQ=255)
                                                             cgACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTg                                                                                     <  1:1775366/148‑1 (MQ=255)
                                                              gACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGt                                                                                    <  2:566713/148‑1 (MQ=255)
                                                              gACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGt                                                                                    <  2:565869/148‑1 (MQ=255)
                                                               aCTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTa                                                                                   >  1:2534830/1‑148 (MQ=255)
                                                               aCTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTa                                                                                   <  1:746116/148‑1 (MQ=255)
                                                                  cGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTAttt                                                                                <  2:1485880/148‑1 (MQ=255)
                                                                     cACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTCTTTCGt                                                                             >  2:4937/1‑148 (MQ=255)
                                                                     cACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGt                                                                             >  2:1679808/1‑148 (MQ=255)
                                                                       ccATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTcc                                                                           >  2:1474598/1‑148 (MQ=255)
                                                                       ccATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTcc                                                                           >  2:677257/1‑148 (MQ=255)
                                                                        cATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCa                                                                          >  2:2555461/1‑148 (MQ=255)
                                                                          tAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATg                                                                        <  1:3144796/148‑1 (MQ=255)
                                                                          tAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATg                                                                        <  1:2949477/148‑1 (MQ=255)
                                                                                   aGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTggg                                                               <  1:2375291/148‑1 (MQ=255)
                                                                                         aGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGtttt                                                         >  1:2851257/1‑148 (MQ=255)
                                                                                           cgcCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtgt                                                       >  1:1051873/1‑148 (MQ=255)
                                                                                               tGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTgtgtgt                                                   >  1:50193/1‑148 (MQ=255)
                                                                                                gCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTc                                                  >  2:3313361/1‑148 (MQ=255)
                                                                                                gCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTc                                                  >  2:3313374/1‑148 (MQ=255)
                                                                                                   agatATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGCTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGcc                                               >  2:2396056/4‑148 (MQ=255)
                                                                                                    gatATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCt                                              >  2:1235383/3‑148 (MQ=255)
                                                                                                    gatATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCt                                              >  2:3220078/3‑148 (MQ=255)
                                                                                                    gatATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCt                                              <  2:2952713/146‑1 (MQ=255)
                                                                                                    gatATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCt                                              <  2:2952453/146‑1 (MQ=255)
                                                                                                     atATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCtt                                             >  2:1890929/2‑148 (MQ=255)
                                                                                                        tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCttttt                                          <  2:2094334/148‑1 (MQ=255)
                                                                                                           ttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCa                                       >  1:1439903/1‑148 (MQ=255)
                                                                                                           ttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCa                                       >  1:2194909/1‑148 (MQ=255)
                                                                                                           ttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCa                                       >  1:189564/1‑148 (MQ=255)
                                                                                                              cGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAgcg                                    <  1:1510020/148‑1 (MQ=255)
                                                                                                                aTGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGcgcg                                  <  2:15691/148‑1 (MQ=255)
                                                                                                                aTGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGcgcg                                  >  1:178536/1‑148 (MQ=255)
                                                                                                                 tGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCCAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCga                                 >  2:518267/1‑148 (MQ=255)
                                                                                                                 tGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCga                                 >  2:518242/1‑148 (MQ=255)
                                                                                                                 tGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCga                                 >  2:518272/1‑148 (MQ=255)
                                                                                                                  gAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCgag                                >  1:114079/1‑148 (MQ=255)
                                                                                                                    aTGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGaa                              <  1:963944/148‑1 (MQ=255)
                                                                                                                     tGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAAt                             >  2:2390391/1‑148 (MQ=255)
                                                                                                                     tGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAAt                             >  2:3230886/1‑148 (MQ=255)
                                                                                                                     tGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAAt                             <  2:2222609/148‑1 (MQ=255)
                                                                                                                     tGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTCGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAAt                             <  2:2742569/148‑1 (MQ=255)
                                                                                                                       ccGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAAtct                           <  1:1716543/148‑1 (MQ=255)
                                                                                                                         gcgcAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATctct                         <  2:2089093/148‑1 (MQ=255)
                                                                                                                          cgcAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTa                        <  2:360483/148‑1 (MQ=255)
                                                                                                                            cAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTAcc                      >  1:1433283/1‑148 (MQ=255)
                                                                                                                              aGGTCTGGCTGGCTCGGTGGCGCAGTTACGCCACTTTGCCCGAAGGTTCGTCCTGAAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTc                    >  2:1355538/1‑148 (MQ=255)
                                                                                                                               ggTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCa                   <  1:144689/148‑1 (MQ=255)
                                                                                                                                gTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCaa                  <  1:451893/148‑1 (MQ=255)
                                                                                                                                 tCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCaaa                 >  1:2783409/1‑148 (MQ=255)
                                                                                                                                  ctggctggCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAAt                <  2:2428278/148‑1 (MQ=255)
                                                                                                                                    ggctggctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATcc              <  1:2701375/148‑1 (MQ=255)
                                                                                                                                       tggctCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCat           <  1:717140/148‑1 (MQ=255)
                                                                                                                                        ggctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCata          >  2:1940561/1‑148 (MQ=255)
                                                                                                                                         gctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCatat         >  1:455973/1‑148 (MQ=255)
                                                                                                                                         gctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCatat         >  1:1973887/1‑148 (MQ=255)
                                                                                                                                         gctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGATCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTATTCAGCGCGAGAATCTCTAACTCAAATCCCatat         >  1:832952/1‑148 (MQ=255)
                                                                                                                                          ctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATc        >  2:2328223/1‑148 (MQ=255)
                                                                                                                                               tGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATaa   <  1:2202828/148‑1 (MQ=255)
                                                                                                                                                ggCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATaaa  <  2:1051873/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAA  >  CP009273/2110245‑2110536

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: