Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,767,000 G→T 100% D454Y (GAC→TAC)  mtlA → fused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,767,0000GT99.1% 1063.8 / ‑0.1 331D454Y (GAC→TAC) mtlAfused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (166/162);  total (167/164)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 338 total aligned reads displayed.

GCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTC  >  CP009273/3766853‑3767146
                                                                                                                                                   |                                                                                                                                                  
 cAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACa                                                                                                                                                 <  1:896945/150‑1 (MQ=255)
  aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAg                                                                                                                                                >  1:1265506/1‑150 (MQ=255)
  aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAg                                                                                                                                                <  2:1468252/150‑1 (MQ=255)
       tGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGcc                                                                                                                                           <  1:1473102/150‑1 (MQ=255)
           gCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTa                                                                                                                                       <  2:1697077/150‑1 (MQ=255)
            cAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTac                                                                                                                                      >  2:1965802/1‑150 (MQ=255)
             aGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTaca                                                                                                                                     <  1:1674223/150‑1 (MQ=255)
             aGATTTCTGTTACTAACAGCGCGATCAACAAACTGCCTCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACTAGCAATGCGCCAGGTTCCGCAGGCACATCATATTTCTCTTACCAACATCCTCTACAGCGGCCTGTaca                                                                                                                                     >  2:892289/1‑150 (MQ=255)
               aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                   >  2:3100161/1‑150 (MQ=255)
               aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                   <  2:3011988/150‑1 (MQ=255)
                    tGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCt                                                                                                                              >  2:2958345/1‑150 (MQ=255)
                     gTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTg                                                                                                                             <  2:3038577/150‑1 (MQ=255)
                      ttACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGa                                                                                                                            >  2:1213107/1‑150 (MQ=255)
                      ttACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGa                                                                                                                            >  2:2147640/1‑150 (MQ=255)
                        aCTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGAcc                                                                                                                          >  1:93577/1‑150 (MQ=255)
                         cTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCg                                                                                                                         >  1:4908517/1‑150 (MQ=255)
                           aaCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACTTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGaa                                                                                                                       >  2:1735774/1‑150 (MQ=255)
                             cAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACg                                                                                                                     >  1:1342198/1‑150 (MQ=255)
                              aGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGt                                                                                                                    >  1:4837557/1‑150 (MQ=255)
                              aGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGt                                                                                                                    >  2:2902509/1‑150 (MQ=255)
                               gcgcGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTc                                                                                                                   >  2:2075372/1‑150 (MQ=255)
                                   gATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                               <  1:3879956/150‑1 (MQ=255)
                                     tCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                                                                             <  2:3235267/150‑1 (MQ=255)
                                          aacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGcc                                                                                                        >  1:977089/1‑150 (MQ=255)
                                          aacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGcc                                                                                                        <  2:1858658/150‑1 (MQ=255)
                                           acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGccc                                                                                                       <  2:3608241/150‑1 (MQ=255)
                                           acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGccc                                                                                                       <  2:4044457/150‑1 (MQ=255)
                                              tGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAAc                                                                                                    >  1:812539/1‑150 (MQ=255)
                                                  gccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                                                                                >  1:4678581/1‑150 (MQ=255)
                                                      gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                            <  1:1034478/150‑1 (MQ=255)
                                                      gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                            >  2:3592591/1‑150 (MQ=255)
                                                       aTGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                                                                           >  2:1628601/1‑150 (MQ=255)
                                                        tgtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGc                                                                                          <  1:3069007/150‑1 (MQ=255)
                                                        tgtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGc                                                                                          >  1:2367955/1‑150 (MQ=255)
                                                           ggACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                       <  1:3721776/150‑1 (MQ=255)
                                                           ggACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                       >  1:4252116/1‑150 (MQ=255)
                                                            gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                      <  1:3304935/150‑1 (MQ=255)
                                                            gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                      <  1:2883704/150‑1 (MQ=255)
                                                              ccTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGa                                                                                    <  1:3966085/150‑1 (MQ=255)
                                                               cTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                                                                   >  1:3675121/1‑150 (MQ=255)
                                                                tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAag                                                                                  <  1:5208623/150‑1 (MQ=255)
                                                                  gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAagag                                                                                <  2:1803502/150‑1 (MQ=255)
                                                                    catcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaa                                                                              >  2:1545513/1‑150 (MQ=255)
                                                                     atcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaaa                                                                             <  2:4539337/150‑1 (MQ=255)
                                                                         ctcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaa                                                                         <  2:4148348/150‑1 (MQ=255)
                                                                         ctcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaa                                                                         >  2:2578679/1‑150 (MQ=255)
                                                                         ctcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaa                                                                         <  2:3947176/150‑1 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                        <  1:3713725/150‑1 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                        <  2:2426355/150‑1 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                        >  2:614211/1‑150 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                        <  1:5282779/150‑1 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                        <  1:3503630/150‑1 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                        <  2:1155602/150‑1 (MQ=255)
                                                                            acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                      >  1:1546381/1‑150 (MQ=255)
                                                                               gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                                                                   >  1:4488720/1‑150 (MQ=255)
                                                                               gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                                                                   <  1:3377484/150‑1 (MQ=255)
                                                                                tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                                  >  1:2681238/1‑150 (MQ=255)
                                                                                tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                                  <  1:229454/150‑1 (MQ=255)
                                                                                 gacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc                                                                 <  1:2242535/150‑1 (MQ=255)
                                                                                 gacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc                                                                 <  1:2270038/150‑1 (MQ=255)
                                                                                  acctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGcc                                                                >  1:5101663/1‑150 (MQ=255)
                                                                                    ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCACACGAAGAGAAAGTAAAAGACAGCCTg                                                              >  2:1461279/1‑150 (MQ=255)
                                                                                    ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTg                                                              >  1:573590/1‑150 (MQ=255)
                                                                                    ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTg                                                              <  1:4105400/150‑1 (MQ=255)
                                                                                        ccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAg                                                          <  2:474348/150‑1 (MQ=255)
                                                                                         cGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGa                                                         <  2:2602401/150‑1 (MQ=255)
                                                                                         cGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGa                                                         >  2:3130795/1‑150 (MQ=255)
                                                                                          gAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGAc                                                        <  2:3681753/150‑1 (MQ=255)
                                                                                           aaCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACa                                                       <  2:4382216/150‑1 (MQ=255)
                                                                                            aCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAg                                                      <  1:4562113/150‑1 (MQ=255)
                                                                                               cgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCtt                                                   >  2:2197563/1‑150 (MQ=255)
                                                                                                 cTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTg                                                 >  1:3629577/1‑150 (MQ=255)
                                                                                                 cTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTg                                                 <  1:4059476/150‑1 (MQ=255)
                                                                                                 cTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTg                                                 >  2:4926936/1‑150 (MQ=255)
                                                                                                    tGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACg                                              <  1:191241/150‑1 (MQ=255)
                                                                                                        ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCAGAAAGACAGCTTTGACGATTc                                          <  1:417765/150‑1 (MQ=255)
                                                                                                        ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCAGAAAGACAGCTTTGACGATTc                                          <  1:3157056/150‑1 (MQ=255)
                                                                                                            gTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt                                      >  1:2936485/1‑150 (MQ=255)
                                                                                                            gTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGt                                      <  1:1779743/150‑1 (MQ=255)
                                                                                                                      cacaGCATATTTCGCTGAACAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt                            >  2:281606/1‑150 (MQ=255)
                                                                                                                       acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                           >  2:5201723/1‑150 (MQ=255)
                                                                                                                       acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                           >  1:433363/1‑150 (MQ=255)
                                                                                                                        caGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTc                          >  1:3807510/1‑150 (MQ=255)
                                                                                                                        caGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCAGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTc                          <  1:2119769/150‑1 (MQ=255)
                                                                                                                         aGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCa                         >  1:1097099/1‑150 (MQ=255)
                                                                                                                           cATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAg                       <  1:758833/150‑1 (MQ=255)
                                                                                                                              atTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTa                    <  2:1665886/150‑1 (MQ=255)
                                                                                                                               tttCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAg                   >  1:3464739/1‑150 (MQ=255)
                                                                                                                                ttCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCAGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAgg                  <  2:4237093/150‑1 (MQ=255)
                                                                                                                                   gCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgcgc               >  1:2061781/1‑150 (MQ=255)
                                                                                                                                     tGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCgg             >  1:2200582/1‑150 (MQ=255)
                                                                                                                                       aCCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGgag           <  1:4964478/150‑1 (MQ=255)
                                                                                                                                        ccAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGgaga          <  2:3825458/150‑1 (MQ=255)
                                                                                                                                        ccAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGgaga          >  2:4877038/1‑150 (MQ=255)
                                                                                                                                          aaCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAAc        >  2:3434321/1‑150 (MQ=255)
                                                                                                                                          aaCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAAc        >  2:3199006/1‑150 (MQ=255)
                                                                                                                                             ttCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATc     <  1:1054950/150‑1 (MQ=255)
                                                                                                                                             ttCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATc     >  2:4617037/1‑150 (MQ=255)
                                                                                                                                               ccTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCtt   <  2:1034582/150‑1 (MQ=255)
                                                                                                                                               ccTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCtt   <  2:965832/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTC  >  CP009273/3766853‑3767146

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: