Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,165,425 C→A 10.5% D136Y (GAC→TAC)  ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,165,4250CA10.5% 899.4 / 54.4 256D136Y (GAC→TAC) ldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (104/125);  new base A (18/9);  total (122/134)
Fisher's exact test for biased strand distribution p-value = 4.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 259 total aligned reads displayed.

CCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACG  >  CP009273/1165277‑1165573
                                                                                                                                                    |                                                                                                                                                    
cccATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTca                                                                                                                                                     >  1:1750946/1‑150 (MQ=255)
     tGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCaccacc                                                                                                                                                <  1:708816/150‑1 (MQ=255)
     tGGGATATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGGGTCAGCGTGTCaccacc                                                                                                                                                <  1:2724845/150‑1 (MQ=255)
       ggTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCTCGAAAGACAGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTa                                                                                                                                              >  2:2504181/1‑150 (MQ=255)
       ggTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTa                                                                                                                                              >  2:2503968/1‑150 (MQ=255)
         ttATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAAc                                                                                                                                            <  2:422400/150‑1 (MQ=255)
           aTCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTg                                                                                                                                          >  2:3757068/1‑150 (MQ=255)
           aTCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTg                                                                                                                                          >  1:110460/1‑150 (MQ=255)
              cagtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGAcc                                                                                                                                       >  1:2843676/1‑150 (MQ=255)
              cagtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCAGTGTCGGTGTCAGCGTGTCACCACCTAACTGAcc                                                                                                                                       >  2:4451874/1‑150 (MQ=255)
                gtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCaa                                                                                                                                     <  1:3315647/150‑1 (MQ=255)
                 tccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAAt                                                                                                                                    <  1:1108541/150‑1 (MQ=255)
                   cagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATAc                                                                                                                                  >  1:4011693/1‑150 (MQ=255)
                   cagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATAc                                                                                                                                  <  1:4197704/150‑1 (MQ=255)
                    agCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATAcc                                                                                                                                 <  1:4353023/150‑1 (MQ=255)
                      ccGGCACTACCGCAGGCAACTCAATTCCCTGTGCTGTATATCGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTAGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTa                                                                                                                               <  2:4408457/150‑1 (MQ=255)
                            cTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGAt                                                                                                                         >  2:179051/1‑150 (MQ=255)
                               ccGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATaca                                                                                                                      <  2:1513229/150‑1 (MQ=255)
                               ccGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATaca                                                                                                                      <  2:1729394/150‑1 (MQ=255)
                                 gCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACg                                                                                                                    <  1:1115155/150‑1 (MQ=255)
                                   aGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGt                                                                                                                  >  2:816887/1‑150 (MQ=255)
                                    ggCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGtt                                                                                                                 >  2:2960028/1‑150 (MQ=255)
                                      cAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTAc                                                                                                               >  2:2152164/1‑150 (MQ=255)
                                       aaCTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTAcc                                                                                                              >  1:960512/1‑150 (MQ=255)
                                       aaCTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTAcc                                                                                                              >  1:2679172/1‑150 (MQ=255)
                                         cTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGa                                                                                                            >  1:2045305/1‑150 (MQ=255)
                                         cTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGa                                                                                                            <  2:3141631/150‑1 (MQ=255)
                                           cAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAt                                                                                                          <  2:345652/150‑1 (MQ=255)
                                               tCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAttttt                                                                                                      <  2:4034769/150‑1 (MQ=255)
                                               tCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAttttt                                                                                                      <  2:4034825/150‑1 (MQ=255)
                                                 ccTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTAc                                                                                                    <  2:4129024/150‑1 (MQ=255)
                                                  cTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTAcc                                                                                                   >  1:2633369/1‑150 (MQ=255)
                                                      gCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggc                                                                                               <  1:1475441/150‑1 (MQ=255)
                                                      gCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggc                                                                                               <  1:1474892/150‑1 (MQ=255)
                                                        tttATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggcgg                                                                                             >  1:3894553/1‑150 (MQ=255)
                                                         ttATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCggg                                                                                            >  2:15814/1‑150 (MQ=255)
                                                           ataACGGGCGGGGATGTTTGCCGTTGGCGTGCTGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTTTCAGCTTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgta                                                                                          <  1:3163179/150‑1 (MQ=255)
                                                           ataACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgta                                                                                          >  2:3173154/1‑150 (MQ=255)
                                                            taACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgtag                                                                                         >  2:3516739/1‑150 (MQ=255)
                                                            taACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgtag                                                                                         >  1:752315/1‑150 (MQ=255)
                                                             aaCGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgtagt                                                                                        >  1:2243582/1‑150 (MQ=255)
                                                                      ggATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGAc                                                                               <  2:364706/150‑1 (MQ=255)
                                                                      ggATGTTGGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGAc                                                                               <  2:2093946/150‑1 (MQ=255)
                                                                       gATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACg                                                                              <  2:3259386/150‑1 (MQ=255)
                                                                        aTGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTGTTACCCGGCGGGTAGTAATAGAGACGc                                                                             >  2:3674446/1‑150 (MQ=255)
                                                                         tGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCa                                                                            >  1:3040975/1‑150 (MQ=255)
                                                                           tttGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGc                                                                          >  2:137748/1‑150 (MQ=255)
                                                                              gCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTcc                                                                       <  1:429071/150‑1 (MQ=255)
                                                                                cGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGc                                                                     >  1:762435/1‑150 (MQ=255)
                                                                                 gTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCa                                                                    <  1:1573233/150‑1 (MQ=255)
                                                                                   tGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAAt                                                                  >  1:1202577/1‑150 (MQ=255)
                                                                                   tGGCGACCAGGTTGGGTTTGCACGTTTGTCTGAAACTGTGGTCACCATTGTCGGAGACAGCGTGAAACCACCTAACTGACCAATACCAAGTGGATACACGGTTACCGTATTTTTACCCGGCGGGTAGGAATAGAGACGCAGCTACGCagt                                                                  >  2:3382837/1‑148 (MQ=255)
                                                                                    ggCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATg                                                                 >  1:583836/1‑150 (MQ=255)
                                                                                    ggCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACAGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATg                                                                 >  1:584863/1‑150 (MQ=255)
                                                                                              ttGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACaa                                                       >  2:1961786/1‑150 (MQ=255)
                                                                                               tGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAAt                                                      >  1:4126784/1‑150 (MQ=255)
                                                                                               tGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATAGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAAt                                                      >  1:1877375/1‑150 (MQ=255)
                                                                                                gggTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATg                                                     >  1:385688/1‑150 (MQ=255)
                                                                                                gggTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATg                                                     >  1:385656/1‑150 (MQ=255)
                                                                                                gggTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATg                                                     >  1:384589/1‑150 (MQ=255)
                                                                                                 ggTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGc                                                    >  2:4517899/1‑150 (MQ=255)
                                                                                                    ttGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCtt                                                 <  2:3216960/150‑1 (MQ=255)
                                                                                                    ttGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCtt                                                 <  1:360731/150‑1 (MQ=255)
                                                                                                    ttGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCtt                                                 <  1:155413/150‑1 (MQ=255)
                                                                                                     tGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTc                                                <  2:948396/150‑1 (MQ=255)
                                                                                                     tGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTc                                                <  1:2056410/150‑1 (MQ=255)
                                                                                                       cACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTcgc                                              >  1:2550778/1‑150 (MQ=255)
                                                                                                             tGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCgcggcg                                        >  1:715773/1‑150 (MQ=255)
                                                                                                              gTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGgcgc                                       >  2:608380/1‑150 (MQ=255)
                                                                                                              gTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGgcgc                                       >  1:2162079/1‑150 (MQ=255)
                                                                                                              gTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGgcgc                                       >  1:2774844/1‑150 (MQ=255)
                                                                                                               tCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCa                                      <  2:1232501/150‑1 (MQ=255)
                                                                                                                   aaaCGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTg                                  <  2:3894540/150‑1 (MQ=255)
                                                                                                                     aCGGTGGTCACCATTGTCGGTGGCAGCGTGTCACCACCTAACTGACCAATACCTATTGGAGACACGGGTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGa                                >  2:4361969/1‑150 (MQ=255)
                                                                                                                      cGGTGGTCACCATTGTAGGAGTCAGCGTGTCACCACCTAACTGACCAATACCGATTGGTTACACGGTTACGGAATTTTTACCCGGCGGGTAGTAGTAGTGTCGCAGGTCCGCGCTGTTGTTGAGAAGGCCTTCGGGGGGCGCagcgggga                               >  2:3006137/1‑143 (MQ=255)
                                                                                                                        gtggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGt                             <  2:232125/150‑1 (MQ=255)
                                                                                                                        gtggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAACTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGt                             <  1:190245/150‑1 (MQ=255)
                                                                                                                        gtggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAACTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGt                             <  1:137618/150‑1 (MQ=255)
                                                                                                                             cACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGt                        <  1:4386793/150‑1 (MQ=255)
                                                                                                                               ccATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGttt                      >  1:2666590/1‑150 (MQ=255)
                                                                                                                               ccATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGttt                      >  1:3392264/1‑150 (MQ=255)
                                                                                                                               ccATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGttt                      >  1:3391792/1‑150 (MQ=255)
                                                                                                                                cATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACCCGGTTACCGAATTTGTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTg                     <  2:4462037/150‑1 (MQ=255)
                                                                                                                                cATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTg                     <  2:595277/150‑1 (MQ=255)
                                                                                                                                 aTTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGACGTAGGGTTTGc                    >  2:3948525/1‑150 (MQ=255)
                                                                                                                                  ttGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCa                   <  1:1425781/150‑1 (MQ=255)
                                                                                                                                    gTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAAc                 >  1:1101172/1‑150 (MQ=255)
                                                                                                                                    gTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAAc                 >  2:1073422/1‑150 (MQ=255)
                                                                                                                                     tCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACg                <  2:821559/150‑1 (MQ=255)
                                                                                                                                     tCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACg                <  1:232140/150‑1 (MQ=255)
                                                                                                                                         tgtCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGAt            <  2:3700523/150‑1 (MQ=255)
                                                                                                                                         tgtCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGAt            <  2:1679850/150‑1 (MQ=255)
                                                                                                                                          gtCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATc           <  2:4004565/150‑1 (MQ=255)
                                                                                                                                             aGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGtt        >  2:3313607/1‑150 (MQ=255)
                                                                                                                                             aGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGtt        <  2:1108714/150‑1 (MQ=255)
                                                                                                                                                  gtCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAacac   >  2:4155714/1‑150 (MQ=255)
                                                                                                                                                   tCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACg  >  2:2820091/1‑150 (MQ=255)
                                                                                                                                                   tCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACg  <  1:3674446/150‑1 (MQ=255)
                                                                                                                                                   tCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACg  <  1:556291/150‑1 (MQ=255)
                                                                                                                                                   tCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTCATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACg  >  2:2334339/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACG  >  CP009273/1165277‑1165573

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: