Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,713,571 G→T 6.3% H97N (CAC→AAC)  anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,713,5710GT 6.3% 755.2 / 17.3 200H97N (CAC→AAC) anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (91/96);  new base T (5/8);  total (96/104)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

GTGCGCCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTAC  >  CP009273/1713424‑1713718
                                                                                                                                                   |                                                                                                                                                   
gTGCGCCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTg                                                                                                                                                     <  1:682054/148‑1 (MQ=255)
 tGCGCCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGa                                                                                                                                                    >  1:352986/1‑148 (MQ=255)
   cgcCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACa                                                                                                                                                  >  1:2348085/1‑148 (MQ=255)
    gcCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACaa                                                                                                                                                 <  1:775167/148‑1 (MQ=255)
    gcCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGAGTGACACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACaa                                                                                                                                                 <  2:1697694/148‑1 (MQ=255)
     ccTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAAc                                                                                                                                                >  2:696322/1‑148 (MQ=255)
     ccTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAAc                                                                                                                                                >  2:695503/1‑148 (MQ=255)
     ccTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAAc                                                                                                                                                >  1:211917/1‑148 (MQ=255)
        tGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGa                                                                                                                                             <  1:1316328/148‑1 (MQ=255)
        tGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGa                                                                                                                                             <  1:1612032/148‑1 (MQ=255)
        tGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGa                                                                                                                                             <  1:1875619/148‑1 (MQ=255)
         gTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGAt                                                                                                                                            <  1:2227472/148‑1 (MQ=255)
         gTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGAt                                                                                                                                            <  2:1658824/148‑1 (MQ=255)
             tccCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCa                                                                                                                                        >  1:797832/1‑148 (MQ=255)
                  aGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTTACAACCGATCGCAACtat                                                                                                                                   >  2:466229/1‑148 (MQ=255)
                  aGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACtat                                                                                                                                   <  1:2653931/148‑1 (MQ=255)
                   ggCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACtata                                                                                                                                  >  1:1121665/1‑148 (MQ=255)
                   ggCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACtata                                                                                                                                  <  1:716567/148‑1 (MQ=255)
                   ggCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACtata                                                                                                                                  <  2:1595890/148‑1 (MQ=255)
                    gCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTatat                                                                                                                                 <  2:778384/148‑1 (MQ=255)
                    gCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTatat                                                                                                                                 >  2:3313371/1‑148 (MQ=255)
                      aaTATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATAtct                                                                                                                               <  2:3059602/148‑1 (MQ=255)
                       atatCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTTACAACCGATCGCAACTATAtctc                                                                                                                              <  1:104503/148‑1 (MQ=255)
                       atatCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATAtctc                                                                                                                              <  2:2766915/148‑1 (MQ=255)
                        tatCCCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATctct                                                                                                                             >  2:411339/1‑148 (MQ=255)
                          tCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTTACAACCGATCGCAACTATATCTCTTg                                                                                                                           <  1:3156173/148‑1 (MQ=255)
                          tCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTg                                                                                                                           >  1:457904/1‑148 (MQ=255)
                          tCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTg                                                                                                                           >  1:329696/1‑148 (MQ=255)
                          tCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTg                                                                                                                           >  1:1662224/1‑148 (MQ=255)
                          tCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTg                                                                                                                           >  1:1546027/1‑148 (MQ=255)
                           cgcgACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGc                                                                                                                          <  2:13994/148‑1 (MQ=255)
                           cgcgACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGc                                                                                                                          >  1:822811/1‑148 (MQ=255)
                           cgcgACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGc                                                                                                                          <  1:1120463/148‑1 (MQ=255)
                           cgcgACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGc                                                                                                                          >  2:1132633/1‑148 (MQ=255)
                            gcgACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGcc                                                                                                                         <  1:2753759/148‑1 (MQ=255)
                               aCGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGc                                                                                                                      >  2:2617507/1‑148 (MQ=255)
                                  gCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGa                                                                                                                   <  1:2572858/148‑1 (MQ=255)
                                  gCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGa                                                                                                                   <  1:2634528/148‑1 (MQ=255)
                                     aaaaTCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGAttt                                                                                                                <  2:2000155/148‑1 (MQ=255)
                                       aaTCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTg                                                                                                              >  2:2687067/1‑148 (MQ=255)
                                           aCCGACAACCGTAAGTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTcc                                                                                                          <  2:1763422/148‑1 (MQ=255)
                                            ccGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCt                                                                                                         <  1:2068907/148‑1 (MQ=255)
                                                   ccGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCaa                                                                                                  >  1:1888360/1‑148 (MQ=255)
                                                    cGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAg                                                                                                 >  2:2195326/1‑148 (MQ=255)
                                                    cGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAg                                                                                                 >  1:2593117/1‑148 (MQ=255)
                                                     gTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAgg                                                                                                <  1:324384/148‑1 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               >  2:2281947/1‑148 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               >  2:2359097/1‑148 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               >  2:2656931/1‑148 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               >  2:1020850/1‑148 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               >  2:2674967/1‑148 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               <  2:2793271/148‑1 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               <  2:2800358/148‑1 (MQ=255)
                                                      tAATTCCGGTGCGTGCCACAATTTGATGGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGc                                                                                               <  2:2345720/148‑1 (MQ=255)
                                                       aaTTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCa                                                                                              >  2:1837236/1‑148 (MQ=255)
                                                        aTTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTTACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCAt                                                                                             >  2:2253149/1‑148 (MQ=255)
                                                         ttCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCAtt                                                                                            >  1:2095914/1‑148 (MQ=255)
                                                         ttCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCAtt                                                                                            >  1:1166152/1‑148 (MQ=255)
                                                          tCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTg                                                                                           <  1:2707058/148‑1 (MQ=255)
                                                             ggTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGAcc                                                                                        >  1:645796/1‑148 (MQ=255)
                                                                 cgtgcCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCAt                                                                                    <  2:2142122/148‑1 (MQ=255)
                                                                 cgagcCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCAt                                                                                    >  1:791112/4‑148 (MQ=255)
                                                                   tgcCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCa                                                                                  >  1:2424620/1‑148 (MQ=255)
                                                                         aaTTTGATTGTTATCTCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCaaaa                                                                            <  2:2936361/148‑1 (MQ=255)
                                                                         aaTTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCaaaa                                                                            <  1:1066775/148‑1 (MQ=255)
                                                                         aaTTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCaaaa                                                                            <  1:610887/148‑1 (MQ=255)
                                                                           tttGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAg                                                                          <  2:3023258/148‑1 (MQ=255)
                                                                           tttGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAg                                                                          <  2:1128818/148‑1 (MQ=255)
                                                                              gATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTg                                                                       <  1:144258/148‑1 (MQ=255)
                                                                              gATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTg                                                                       <  1:3190867/148‑1 (MQ=255)
                                                                               aTTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGc                                                                      >  2:2019843/1‑148 (MQ=255)
                                                                                ttGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGcc                                                                     >  1:87571/1‑148 (MQ=255)
                                                                                 tGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGccc                                                                    <  2:1688706/148‑1 (MQ=255)
                                                                                 tGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGccc                                                                    >  2:2844157/1‑148 (MQ=255)
                                                                                    tATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAg                                                                 <  1:2687067/148‑1 (MQ=255)
                                                                                       cGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTg                                                              <  1:3031942/148‑1 (MQ=255)
                                                                                          cAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGt                                                           >  1:2785828/1‑148 (MQ=255)
                                                                                           aaTCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTa                                                          >  1:2349576/1‑148 (MQ=255)
                                                                                             tCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTTACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATc                                                        <  1:466229/148‑1 (MQ=255)
                                                                                              cTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCa                                                       <  2:822811/148‑1 (MQ=255)
                                                                                               tGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCaa                                                      <  1:133630/148‑1 (MQ=255)
                                                                                                 tAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGc                                                    <  2:337931/148‑1 (MQ=255)
                                                                                                     gtgtgtGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTc                                                <  1:375506/148‑1 (MQ=255)
                                                                                                     gtgtgtGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACAGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTc                                                >  2:3156139/1‑148 (MQ=255)
                                                                                                     gtgtgtGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACAGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTc                                                >  2:3156052/1‑148 (MQ=255)
                                                                                                                  cGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAg                                   >  2:1674977/1‑148 (MQ=255)
                                                                                                                  cGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAg                                   >  1:392486/1‑148 (MQ=255)
                                                                                                                        tCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAg                             <  1:695602/148‑1 (MQ=255)
                                                                                                                         cGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAgc                            >  1:2774904/1‑148 (MQ=255)
                                                                                                                           gTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAgctg                          <  1:1914478/148‑1 (MQ=255)
                                                                                                                            ttCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAgctgc                         >  2:3100017/1‑148 (MQ=255)
                                                                                                                            ttCATGCCAGACGGTTTGACAGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAgctgc                         <  1:3156139/148‑1 (MQ=255)
                                                                                                                              cATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGctgctg                       >  1:2820114/1‑148 (MQ=255)
                                                                                                                              cATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGctgctg                       <  1:1193628/148‑1 (MQ=255)
                                                                                                                               aTGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTgg                      >  1:3271019/1‑148 (MQ=255)
                                                                                                                                       cGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCa              >  2:1458930/1‑148 (MQ=255)
                                                                                                                                          tttGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAat           >  1:2788207/1‑148 (MQ=255)
                                                                                                                                               ccGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCa      >  1:809262/1‑148 (MQ=255)
                                                                                                                                                cGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAg     <  2:352986/148‑1 (MQ=255)
                                                                                                                                                cGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAg     <  2:1888360/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGCGCCTTGTCCTCCCAAGGCAATATCGCGACGGCGAAAATCACCGACAACCGTAATTCCGGTGCGTGCCACAATTTGATTGTTATCGCCAATCTGTAAAGTGTGTGGTGCCACGCCCGTCGGTTCATGCCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTAC  >  CP009273/1713424‑1713718

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: