Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,388,406 G→T 100% R308L (CGC→CTC)  mppA → murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,388,4060GT95.1% 1085.2 / 39.5 325R308L (CGC→CTC) mppAmurein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (3/13);  new base T (168/141);  total (171/154)
Fisher's exact test for biased strand distribution p-value = 8.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 332 total aligned reads displayed.

CGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAATGAGTC  >  CP009273/1388259‑1388552
                                                                                                                                                   |                                                                                                                                                  
   gtgATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTc                                                                                                                                               <  1:4455409/148‑1 (MQ=255)
   gggATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTc                                                                                                                                               >  2:2834678/1‑150 (MQ=255)
     gATATTGATATCACCGAATCCTTCCCGAACAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCtcgttcgt                                                                                                                                             <  2:2156587/150‑1 (MQ=255)
     gATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCtcgttcgt                                                                                                                                             >  2:4448662/1‑150 (MQ=255)
       tatTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCt                                                                                                                                           <  2:4102233/150‑1 (MQ=255)
          tGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGc                                                                                                                                        >  1:1612389/1‑150 (MQ=255)
            atatCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCAt                                                                                                                                      >  1:173708/1‑150 (MQ=255)
            atatCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCAt                                                                                                                                      >  1:2170904/1‑150 (MQ=255)
             tatCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCAtt                                                                                                                                     <  1:2661744/150‑1 (MQ=255)
                cACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAg                                                                                                                                  <  1:1265821/150‑1 (MQ=255)
                 aCCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGt                                                                                                                                 >  1:1865015/1‑150 (MQ=255)
                 aCCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGt                                                                                                                                 <  2:3486038/150‑1 (MQ=255)
                 aCCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATAGTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGt                                                                                                                                 >  2:2066879/1‑150 (MQ=255)
                   cGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTAt                                                                                                                               >  1:2016964/1‑150 (MQ=255)
                   cGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTAt                                                                                                                               >  2:2785282/1‑150 (MQ=255)
                    gAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATg                                                                                                                              >  1:194516/1‑150 (MQ=255)
                    gAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCAGTAAGTATg                                                                                                                              >  2:893840/1‑150 (MQ=255)
                      aTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGAc                                                                                                                            <  2:3761038/150‑1 (MQ=255)
                      aTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGAc                                                                                                                            <  2:1351081/150‑1 (MQ=255)
                      aTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGAc                                                                                                                            >  1:1656514/1‑150 (MQ=255)
                      aTCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCGATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGAc                                                                                                                            <  2:2305572/150‑1 (MQ=255)
                       tCCGTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACg                                                                                                                           <  2:3856591/150‑1 (MQ=255)
                                aaaaaTATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATcgc                                                                                                                  >  1:2978692/1‑150 (MQ=255)
                                 aaaaTATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATcgcc                                                                                                                 <  1:1379420/150‑1 (MQ=255)
                                  aaaTATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATcgccg                                                                                                                >  2:3604556/1‑150 (MQ=255)
                                  aaaTATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATcgccg                                                                                                                >  2:3210875/1‑150 (MQ=255)
                                    atatGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCgccgcc                                                                                                              >  1:3909649/1‑150 (MQ=255)
                                       tGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtga                                                                                                           <  2:1718177/150‑1 (MQ=255)
                                        gTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgat                                                                                                          <  1:712369/150‑1 (MQ=255)
                                         tATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatg                                                                                                         >  1:2178602/1‑150 (MQ=255)
                                         tATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatg                                                                                                         <  1:585700/150‑1 (MQ=255)
                                         tATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatg                                                                                                         >  2:3227438/1‑150 (MQ=255)
                                         tATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatg                                                                                                         >  2:842631/1‑150 (MQ=255)
                                          aTCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCtgatga                                                                                                        <  1:3901473/150‑1 (MQ=255)
                                            cAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGAcc                                                                                                      >  2:4629733/1‑150 (MQ=255)
                                             aGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCg                                                                                                     >  1:3050272/1‑150 (MQ=255)
                                              gAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGa                                                                                                    >  2:3569181/1‑150 (MQ=255)
                                              gAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGa                                                                                                    >  2:684950/1‑150 (MQ=255)
                                               aaGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGaa                                                                                                   <  2:2136298/150‑1 (MQ=255)
                                                 gCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGaaaa                                                                                                 <  1:4663645/150‑1 (MQ=255)
                                                 gCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGaaaa                                                                                                 >  1:3347800/1‑150 (MQ=255)
                                                   tGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAg                                                                                               >  1:2554020/1‑150 (MQ=255)
                                                   tGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAg                                                                                               >  1:1032756/1‑150 (MQ=255)
                                                        aaGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTa                                                                                          >  1:514348/1‑150 (MQ=255)
                                                          ggATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAgg                                                                                        <  2:2468282/150‑1 (MQ=255)
                                                          ggATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAgg                                                                                        <  1:3389639/150‑1 (MQ=255)
                                                          ggATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAgg                                                                                        >  2:3451976/1‑150 (MQ=255)
                                                           gATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAggg                                                                                       >  1:1811609/1‑150 (MQ=255)
                                                             tatTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGAc                                                                                     <  1:3607305/150‑1 (MQ=255)
                                                                   ggggCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGa                                                                               >  2:4734158/1‑150 (MQ=255)
                                                                   ggggCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGa                                                                               >  1:3453174/1‑150 (MQ=255)
                                                                    gggCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGaa                                                                              >  2:1409519/1‑150 (MQ=255)
                                                                      gCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGaaaa                                                                            <  1:2112697/150‑1 (MQ=255)
                                                                       cAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAg                                                                           <  2:1606642/150‑1 (MQ=255)
                                                                        aGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGc                                                                          <  1:4240840/150‑1 (MQ=255)
                                                                          gTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCa                                                                        >  2:4102518/1‑150 (MQ=255)
                                                                             tataCGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCg                                                                     <  1:3962102/150‑1 (MQ=255)
                                                                               taCGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTg                                                                   <  1:842631/150‑1 (MQ=255)
                                                                               taCGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTg                                                                   <  1:3227438/150‑1 (MQ=255)
                                                                                aCGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTgg                                                                  <  2:168744/150‑1 (MQ=255)
                                                                                  gccgccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCa                                                                >  2:956465/1‑150 (MQ=255)
                                                                                  gccgccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCa                                                                >  2:3352441/1‑150 (MQ=255)
                                                                                   ccgccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCAt                                                               <  2:4291445/150‑1 (MQ=255)
                                                                                    cgccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCAtt                                                              >  2:2875996/1‑150 (MQ=255)
                                                                                     gccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCAttt                                                             <  2:2550933/150‑1 (MQ=255)
                                                                                     gccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCAttt                                                             <  2:4001000/150‑1 (MQ=255)
                                                                                     gccgcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGAAttt                                                             <  1:1366416/150‑1 (MQ=255)
                                                                                       cgcAGCTCGGGAGCTATTTTTTTGCGTTTAACACGCAAATAGGGCCGTCGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGTTCGCCGCCTGATGAGCGAAAAAGTTTTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTa                                                           <  1:3449936/150‑1 (MQ=255)
                                                                                        gcAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTac                                                          >  2:438556/1‑150 (MQ=255)
                                                                                            cTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCa                                                      >  1:4190491/1‑150 (MQ=255)
                                                                                                 gACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGt                                                 >  2:3694728/1‑150 (MQ=255)
                                                                                                    cTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCAtttt                                                            >  1:4547688/1‑136 (MQ=255)
                                                                                                    cTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCAtttt                                                            <  2:4547688/136‑1 (MQ=255)
                                                                                                    cTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTAc                                              <  2:529119/150‑1 (MQ=255)
                                                                                                    cTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTAc                                              <  2:777702/150‑1 (MQ=255)
                                                                                                      attattatGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACcg                                            <  1:2066621/150‑1 (MQ=255)
                                                                                                      attattatGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACcg                                            >  1:3484915/1‑150 (MQ=255)
                                                                                                          ttatGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCggg                                        >  1:2372125/1‑150 (MQ=255)
                                                                                                           tatGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGa                                       <  2:4662518/150‑1 (MQ=255)
                                                                                                             tGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGAtt                                     >  2:1007961/1‑150 (MQ=255)
                                                                                                             tGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGAtt                                     >  2:2975107/1‑150 (MQ=255)
                                                                                                                 tttAACACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACg                                 >  2:3097654/1‑150 (MQ=255)
                                                                                                                    aaCACGCAAAAAGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCg                              <  2:149099/150‑1 (MQ=255)
                                                                                                                           aaaaaGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCtt                       >  2:1644350/1‑150 (MQ=255)
                                                                                                                           aaaaaGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCtt                       <  1:4385309/150‑1 (MQ=255)
                                                                                                                            aaaaGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTc                      <  2:1038336/150‑1 (MQ=255)
                                                                                                                               aGGGCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGcc                   >  2:1099446/1‑150 (MQ=255)
                                                                                                                                tggCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCg                  >  1:1474274/2‑150 (MQ=255)
                                                                                                                                  gCCGACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGtt                <  2:3996714/150‑1 (MQ=255)
                                                                                                                                     gACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGa             >  1:276283/1‑150 (MQ=255)
                                                                                                                                     gACGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGa             <  1:4734158/150‑1 (MQ=255)
                                                                                                                                       cGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAAc           <  2:2435705/150‑1 (MQ=255)
                                                                                                                                       cGGCAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGGTTGAAc           >  2:4030971/1‑150 (MQ=255)
                                                                                                                                          cAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACaaa        >  2:3034268/1‑150 (MQ=255)
                                                                                                                                          cAGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACaaa        >  1:2827779/1‑150 (MQ=255)
                                                                                                                                           aGATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAAt       >  2:1167720/1‑150 (MQ=255)
                                                                                                                                            gATCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATTGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAATg      >  1:3106100/1‑150 (MQ=255)
                                                                                                                                             aTCAGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAATGa     >  2:35133/1‑150 (MQ=255)
                                                                                                                                                aGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAATGAGTc  <  1:2435012/150‑1 (MQ=255)
                                                                                                                                                aGCTCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAATGAGTc  >  2:1916998/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGGGGATATTGATATCACCGAATCCTTCCCGAAAAATATGTATCAGAAGCTGTTGAAGGATATTCCGGGGCAGGTTTATACGCCGCCGCAGCTCGGGACCTATTATTATGCGTTTAACACGCAAAAAGGGCCGACGGCAGATCAGCGCGTTCGTCTGGCATTAAGTATGACGATAGATCGCCGCCTGATGACCGAAAAAGTATTAGGGACGGGCGAAAAGCCAGCGTGGCATTTTACACCAGATGTTACCGCGGGATTTACGCCGGAACCTTCGCCGTTTGAACAAATGAGTC  >  CP009273/1388259‑1388552

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: