Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,264,143 G→A 100% intergenic (+119/‑62) mepS → / → rtn murein DD‑endopeptidase, space‑maker hydrolase, mutational suppressor of prc thermosensitivity, outer membrane lipoprotein/resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,264,1430GA99.1% 802.9 / ‑2.3 226intergenic (+119/‑62)mepS/rtnmurein DD‑endopeptidase, space‑maker hydrolase, mutational suppressor of prc thermosensitivity, outer membrane lipoprotein/resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (117/107);  total (119/107)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 227 total aligned reads displayed.

AAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCC  >  CP009273/2263996‑2264292
                                                                                                                                                   |                                                                                                                                                     
aaGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACt                                                                                                                                                     <  1:166925/150‑1 (MQ=255)
aaGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACt                                                                                                                                                     <  2:2964633/150‑1 (MQ=255)
 aGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACta                                                                                                                                                    >  2:3424714/1‑150 (MQ=255)
   cACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACtata                                                                                                                                                  >  1:2795036/1‑150 (MQ=255)
    aCGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAg                                                                                                                                                 >  1:1368379/1‑150 (MQ=255)
      gCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTc                                                                                                                                               >  2:306087/1‑150 (MQ=255)
       ccGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTcc                                                                                                                                              >  1:2457787/1‑150 (MQ=255)
       ccGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTcc                                                                                                                                              <  2:3277446/150‑1 (MQ=255)
       ccGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTcc                                                                                                                                              <  1:2629218/150‑1 (MQ=255)
        cGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCg                                                                                                                                             <  2:1416470/150‑1 (MQ=255)
           gTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTc                                                                                                                                          >  1:885907/1‑150 (MQ=255)
              ctcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCAtt                                                                                                                                       <  2:2940795/150‑1 (MQ=255)
                   ccGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTAt                                                                                                                                  >  1:3746486/1‑150 (MQ=255)
                       aGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCttt                                                                                                                              <  1:276721/150‑1 (MQ=255)
                        gCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCtttt                                                                                                                             >  2:2920950/1‑150 (MQ=255)
                         cTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTa                                                                                                                            >  1:3207994/1‑150 (MQ=255)
                          taataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTaa                                                                                                                           >  2:1331997/1‑150 (MQ=255)
                            ataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAAt                                                                                                                         >  1:1836028/1‑150 (MQ=255)
                             taaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAAtg                                                                                                                        >  1:240108/1‑150 (MQ=255)
                              aaaCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAAtgt                                                                                                                       <  2:3613672/150‑1 (MQ=255)
                                aCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAAtgttg                                                                                                                     <  1:2383490/150‑1 (MQ=255)
                                 ccGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAAtgttgt                                                                                                                    <  1:3766366/150‑1 (MQ=255)
                                  cGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACGGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                                                                                   <  2:3006305/150‑1 (MQ=255)
                                   gTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCg                                                                                                                  >  1:3490947/1‑150 (MQ=255)
                                    tttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGt                                                                                                                 <  2:3393944/150‑1 (MQ=255)
                                      tGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTaa                                                                                                               <  1:1273409/150‑1 (MQ=255)
                                       ggATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCCGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAAt                                                                                                              >  2:2898758/1‑150 (MQ=255)
                                       ggATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAAt                                                                                                              <  1:3505616/150‑1 (MQ=255)
                                         aTGCAATCACTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAAttt                                                                                                            <  1:223295/150‑1 (MQ=255)
                                            cAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAg                                                                                                         >  1:2950404/1‑150 (MQ=255)
                                              aTCCCTTGGCTATCCTTACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGa                                                                                                       <  1:3531518/150‑1 (MQ=255)
                                              aTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGa                                                                                                       >  1:241668/1‑150 (MQ=255)
                                               tCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGaa                                                                                                      <  1:2021225/150‑1 (MQ=255)
                                                cccTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGaaa                                                                                                     >  1:801384/1‑150 (MQ=255)
                                                    tGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTaa                                                                                                 >  1:3305594/1‑150 (MQ=255)
                                                      gCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACg                                                                                               >  1:2888239/1‑150 (MQ=255)
                                                       cTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACgg                                                                                              >  1:467232/1‑150 (MQ=255)
                                                        tATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGa                                                                                             >  2:1115252/1‑150 (MQ=255)
                                                         aTCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATTTTTCCGCGTCTTGATACAGACTATAGTCCTGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAAAGGaa                                                                                            >  2:2142097/1‑150 (MQ=255)
                                                         aTCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTATTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCATGATCCAGACTATAGTCCGGTCATTGTTATCTATTAAATGTTGTCGTAATTTCAGGAAATTAACGTaa                                                                                            >  1:36229/1‑150 (MQ=255)
                                                             tGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCAt                                                                                        >  2:3427411/1‑150 (MQ=255)
                                                              gACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCTTGAAATTAACGGAATCATg                                                                                       >  2:2243022/1‑150 (MQ=255)
                                                               aCGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGt                                                                                      <  2:2023423/150‑1 (MQ=255)
                                                                 gAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTc                                                                                    <  2:1019687/150‑1 (MQ=255)
                                                                  aGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCa                                                                                   >  1:3840175/1‑150 (MQ=255)
                                                                  aGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCa                                                                                   >  2:425801/1‑150 (MQ=255)
                                                                     tAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAc                                                                                >  1:1075713/1‑150 (MQ=255)
                                                                     tAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAc                                                                                >  1:48747/1‑150 (MQ=255)
                                                                         tGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGTAATTAACTGAATCATGTTCATACgcgc                                                                            >  2:3742255/1‑150 (MQ=255)
                                                                         tGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACgcgc                                                                            >  2:2521826/1‑150 (MQ=255)
                                                                         tGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACgcgc                                                                            >  2:1267920/1‑150 (MQ=255)
                                                                          gAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCt                                                                           <  2:2978549/150‑1 (MQ=255)
                                                                                    cTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTg                                                                 <  2:1483381/150‑1 (MQ=255)
                                                                                       aGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGAc                                                              <  2:953653/150‑1 (MQ=255)
                                                                                          cAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTa                                                           <  1:2272466/150‑1 (MQ=255)
                                                                                          cAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTa                                                           >  1:317258/1‑150 (MQ=255)
                                                                                          cAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTa                                                           >  2:1808471/1‑150 (MQ=255)
                                                                                           aGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTaa                                                          >  1:3109503/1‑150 (MQ=255)
                                                                                           aGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTaa                                                          >  1:4007336/1‑150 (MQ=255)
                                                                                             tGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGc                                                        <  2:2132098/150‑1 (MQ=255)
                                                                                              gCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCt                                                       <  1:1090637/150‑1 (MQ=255)
                                                                                               cTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTc                                                      >  1:3897892/1‑150 (MQ=255)
                                                                                               cTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTc                                                      >  1:699679/1‑150 (MQ=255)
                                                                                                ttttttGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTcc                                                     <  2:2363480/150‑1 (MQ=255)
                                                                                                    ttGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCt                                                 >  2:1296222/1‑150 (MQ=255)
                                                                                                    ttGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCt                                                 >  1:436684/1‑150 (MQ=255)
                                                                                                    ttGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCt                                                 >  1:656063/1‑150 (MQ=255)
                                                                                                           cattcatCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGc                                          >  1:3560868/1‑150 (MQ=255)
                                                                                                             ttcatCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCAt                                        >  2:2737942/1‑150 (MQ=255)
                                                                                                              tcatcaGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCAtt                                       >  2:477971/1‑150 (MQ=255)
                                                                                                              tcatcaGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCAtt                                       <  2:2204426/150‑1 (MQ=255)
                                                                                                                  caGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCAttgttg                                   >  2:2646717/1‑150 (MQ=255)
                                                                                                                  caGAGAAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCAttgttg                                   >  1:2060170/1‑150 (MQ=255)
                                                                                                                    gagaAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCa                                 <  2:1683190/150‑1 (MQ=255)
                                                                                                                     agaAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAg                                <  1:2291083/150‑1 (MQ=255)
                                                                                                                      gaAAATGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAgg                               >  2:2597770/1‑150 (MQ=255)
                                                                                                                       aaaaTGATGTTTCCGCGTCTTGATACAGACTATAGTCCGGTCATTGTTATCTTTTATATGTTGTCGTAATTTCAGGAAATTAACGGAATCAGGTTCATAAGCTCTCCCAATGTTGGACTTAAGCTCGTGCTTACCTGCATTGTTGCAGGc                              >  2:3644051/1‑150 (MQ=255)
                                                                                                                        aaaTGATGTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCg                             <  1:1480090/150‑1 (MQ=255)
                                                                                                                              tgTTTCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAACTTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGa                       >  2:2746143/1‑150 (MQ=255)
                                                                                                                                tttCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGAtt                     <  2:161507/150‑1 (MQ=255)
                                                                                                                                tttCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGAtt                     <  1:2977023/150‑1 (MQ=255)
                                                                                                                                 ttCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTg                    >  2:911254/1‑150 (MQ=255)
                                                                                                                                 ttCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTg                    <  2:876319/150‑1 (MQ=255)
                                                                                                                                  tCCGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGc                   <  2:4010738/150‑1 (MQ=255)
                                                                                                                                   ccGCGTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCg                  >  1:1313695/1‑150 (MQ=255)
                                                                                                                                     gcgTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGAt                <  2:117349/150‑1 (MQ=255)
                                                                                                                                      cgTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATa               <  1:1122917/150‑1 (MQ=255)
                                                                                                                                      cgTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATa               >  1:3408170/1‑150 (MQ=255)
                                                                                                                                      cgTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTGGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATa               <  1:849705/150‑1 (MQ=255)
                                                                                                                                       gTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATAc              <  2:3495531/150‑1 (MQ=255)
                                                                                                                                       gTCTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATAc              >  2:3747251/1‑150 (MQ=255)
                                                                                                                                         cTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTg            <  2:218258/150‑1 (MQ=255)
                                                                                                                                         cTTGATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTg            >  2:1279026/1‑150 (MQ=255)
                                                                                                                                            gATCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGAATTGTTGCAGGCGTAAGGATTGCGATACTGGTg         >  1:2567092/1‑150 (MQ=255)
                                                                                                                                              tCCAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAg       >  2:1212289/1‑150 (MQ=255)
                                                                                                                                               ccAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGt      >  1:1958499/1‑150 (MQ=255)
                                                                                                                                               ccAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGt      >  2:2506652/1‑150 (MQ=255)
                                                                                                                                                cAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGtt     >  2:3593465/1‑150 (MQ=255)
                                                                                                                                                cAGACTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGtt     <  2:543306/150‑1 (MQ=255)
                                                                                                                                                   gCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGcc  <  2:354578/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCC  >  CP009273/2263996‑2264292

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: