Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,431 C→A 100% P90T (CCA→ACA)  yiaN → 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,738,4310CA98.2% 964.2 / 5.8 284P90T (CCA→ACA) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (162/117);  total (166/118)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
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 286 total aligned reads displayed.

CAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGG  >  CP009273/3738284‑3738577
                                                                                                                                                   |                                                                                                                                                  
 aGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                                 >  1:4093819/1‑150 (MQ=255)
 aGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                                 >  2:4148660/1‑150 (MQ=255)
 aGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                                 >  1:3514151/1‑150 (MQ=255)
 aGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                                 >  2:167973/1‑150 (MQ=255)
     cATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgg                                                                                                                                             <  2:3814096/150‑1 (MQ=255)
     cATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgg                                                                                                                                             >  1:3276125/1‑150 (MQ=255)
     cATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgg                                                                                                                                             >  1:1847023/1‑150 (MQ=255)
             aaaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggct                                                                                                                                     >  2:4483212/1‑150 (MQ=255)
             aaaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggct                                                                                                                                     <  1:4337412/150‑1 (MQ=255)
             aaaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggct                                                                                                                                     >  1:2866341/1‑150 (MQ=255)
             aaaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggct                                                                                                                                     >  2:839042/1‑150 (MQ=255)
              aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTa                                                                                                                                    >  1:4578160/1‑150 (MQ=255)
               aCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTAc                                                                                                                                   <  2:1978248/150‑1 (MQ=255)
                  cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTc                                                                                                                                >  2:2476916/1‑150 (MQ=255)
                  cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTc                                                                                                                                >  2:730369/1‑150 (MQ=255)
                   tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCg                                                                                                                               <  1:3326172/150‑1 (MQ=255)
                    ggTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCgg                                                                                                                              <  2:2141002/150‑1 (MQ=255)
                     gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGc                                                                                                                             >  1:2513138/1‑150 (MQ=255)
                      tGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCg                                                                                                                            <  2:107429/150‑1 (MQ=255)
                       gAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGTGCTGGTCTACGTCGGCGt                                                                                                                           >  2:756032/1‑150 (MQ=255)
                          cGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCt                                                                                                                        <  1:2353546/150‑1 (MQ=255)
                          cGGAGCCGATAGCTTCCCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCt                                                                                                                        <  1:3872703/150‑1 (MQ=255)
                            gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                      >  2:643581/1‑150 (MQ=255)
                             aGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                     >  2:811039/1‑150 (MQ=255)
                             aGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                     >  1:3595460/1‑150 (MQ=255)
                                 gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                 >  2:4279785/1‑150 (MQ=255)
                                 gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                 >  1:3203360/1‑150 (MQ=255)
                                 gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                 >  2:1629825/1‑150 (MQ=255)
                                   tAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTAt                                                                                                               <  2:1821705/150‑1 (MQ=255)
                                   tAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTAt                                                                                                               <  1:4172743/150‑1 (MQ=255)
                                       ttCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGAtt                                                                                                           <  2:3143555/150‑1 (MQ=255)
                                           cccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATgg                                                                                                       <  1:2727453/150‑1 (MQ=255)
                                            ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGc                                                                                                      <  2:3974408/150‑1 (MQ=255)
                                               cctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                   >  2:4096014/2‑150 (MQ=255)
                                                ctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGc                                                                                                  >  1:3781703/1‑150 (MQ=255)
                                                  ggCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCt                                                                                                >  1:258194/1‑150 (MQ=255)
                                                  ggCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCt                                                                                                >  2:3069574/1‑150 (MQ=255)
                                                    cGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCttt                                                                                              <  2:2679120/150‑1 (MQ=255)
                                                     gATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTc                                                                                             <  1:3099756/150‑1 (MQ=255)
                                                        tCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTgg                                                                                          <  2:3514462/150‑1 (MQ=255)
                                                        tCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTgg                                                                                          <  2:2579926/150‑1 (MQ=255)
                                                          cGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                        <  2:893415/150‑1 (MQ=255)
                                                          cGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                        <  2:3648951/150‑1 (MQ=255)
                                                           gTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTc                                                                                       >  2:1725033/1‑150 (MQ=255)
                                                               tttGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGcc                                                                                   <  2:3653521/150‑1 (MQ=255)
                                                               tttGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGcc                                                                                   <  1:1141271/150‑1 (MQ=255)
                                                                 tgtgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGt                                                                                 >  2:437265/1‑150 (MQ=255)
                                                                 tgtgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGt                                                                                 >  2:2968862/1‑150 (MQ=255)
                                                                  gtgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAATAGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGtt                                                                                <  1:2236532/150‑1 (MQ=255)
                                                                  gtgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTTCGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGtt                                                                                <  1:904846/150‑1 (MQ=255)
                                                                   tgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTg                                                                               <  1:3555153/150‑1 (MQ=255)
                                                                      tGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCgg                                                                            >  1:2469401/1‑150 (MQ=255)
                                                                        gCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGACGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGAt                                                                          >  2:3392106/1‑150 (MQ=255)
                                                                              gAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgcc                                                                    >  2:854352/1‑150 (MQ=255)
                                                                                aaTCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgc                                                                  >  1:4243069/1‑150 (MQ=255)
                                                                                 aTCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgcc                                                                 <  2:2745797/150‑1 (MQ=255)
                                                                                    aTGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGtt                                                              <  1:862021/150‑1 (MQ=255)
                                                                                      gAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGc                                                            >  1:2147522/1‑150 (MQ=255)
                                                                                            gggcggCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGctct                                                      >  1:1492773/1‑150 (MQ=255)
                                                                                            gggcggCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGctct                                                      >  2:1124381/1‑150 (MQ=255)
                                                                                             ggcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTa                                                     >  2:3872012/1‑150 (MQ=255)
                                                                                               cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACt                                                   <  2:3502229/150‑1 (MQ=255)
                                                                                               cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACt                                                   >  2:1487762/1‑150 (MQ=255)
                                                                                                ggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                  >  1:4384328/1‑150 (MQ=255)
                                                                                                 gcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTgg                                                 >  2:4070751/1‑150 (MQ=255)
                                                                                                  cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGt                                                <  1:38363/150‑1 (MQ=255)
                                                                                                     gTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGcc                                             >  2:1467954/1‑150 (MQ=255)
                                                                                                       cAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCga                                           <  2:475411/150‑1 (MQ=255)
                                                                                                       cAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCga                                           >  1:1872437/1‑150 (MQ=255)
                                                                                                       cAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCga                                           >  1:136723/1‑150 (MQ=255)
                                                                                                         aaaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                         >  2:1382070/1‑150 (MQ=255)
                                                                                                         aaaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                         >  2:2680936/1‑150 (MQ=255)
                                                                                                          aaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatga                                        >  1:1443599/1‑150 (MQ=255)
                                                                                                              cATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcg                                    >  1:2169944/1‑150 (MQ=255)
                                                                                                                 tgttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAg                                 >  1:3228661/1‑150 (MQ=255)
                                                                                                                  gttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgc                                <  2:3330916/150‑1 (MQ=255)
                                                                                                                  gttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgc                                <  2:2081255/150‑1 (MQ=255)
                                                                                                                   ttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgcg                               <  1:2946630/150‑1 (MQ=255)
                                                                                                                    tgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgcgc                              >  2:598668/1‑150 (MQ=255)
                                                                                                                    tgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgcgc                              >  1:3716140/1‑150 (MQ=255)
                                                                                                                     gACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGcc                             >  2:4063612/1‑150 (MQ=255)
                                                                                                                      aCCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCa                            >  2:1656146/1‑150 (MQ=255)
                                                                                                                      aCCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCa                            <  1:1600615/150‑1 (MQ=255)
                                                                                                                      aCCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCa                            <  1:3644305/150‑1 (MQ=255)
                                                                                                                       ccTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCaa                           >  2:4430962/1‑150 (MQ=255)
                                                                                                                        cTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAAc                          >  1:838588/1‑150 (MQ=255)
                                                                                                                         tGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACt                         <  2:3726944/150‑1 (MQ=255)
                                                                                                                             gATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATcc                     >  1:851414/1‑150 (MQ=255)
                                                                                                                             gATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATcc                     >  1:1025398/1‑150 (MQ=255)
                                                                                                                                gAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGt                  >  2:3183905/1‑150 (MQ=255)
                                                                                                                                gAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGt                  <  1:4006277/150‑1 (MQ=255)
                                                                                                                                   aCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCaa               <  1:1071109/150‑1 (MQ=255)
                                                                                                                                    cTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAAc              >  1:3132328/1‑150 (MQ=255)
                                                                                                                                    cTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAAc              <  1:4270675/150‑1 (MQ=255)
                                                                                                                                     tggtggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAAcc             >  2:2385325/1‑150 (MQ=255)
                                                                                                                                      ggtggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCg            <  2:3595460/150‑1 (MQ=255)
                                                                                                                                        tggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCggg          <  2:4204007/150‑1 (MQ=255)
                                                                                                                                        tggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCggg          <  2:2980016/150‑1 (MQ=255)
                                                                                                                                              aTAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGcggcgg    <  1:2973036/150‑1 (MQ=255)
                                                                                                                                                aaaCCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCgggg  >  2:2791674/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGG  >  CP009273/3738284‑3738577

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: