Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,036,102 C→A 12.0% S85* (TCG→TAG)  yodB → cytochrome b561 homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,036,1020CA12.0% 1065.6 / 76.8 324S85* (TCG→TAG) yodBcytochrome b561 homolog
Reads supporting (aligned to +/- strand):  ref base C (158/127);  new base A (16/23);  total (174/150)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 329 total aligned reads displayed.

TAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAA  >  CP009273/2035953‑2036251
                                                                                                                                                     |                                                                                                                                                     
tagtagtaCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTc                                                                                                                                                       <  2:2194371/150‑1 (MQ=255)
  gtagtaATTATCTGATGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATTTTTTCACGCCTGATTATAAAACACCGTTATTGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGc                                                                                                                                                     <  1:4250661/150‑1 (MQ=255)
    agtaCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTa                                                                                                                                                   <  1:2979177/150‑1 (MQ=255)
    agtaCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTa                                                                                                                                                   >  2:83239/1‑150 (MQ=255)
     gtaCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTaa                                                                                                                                                  >  1:1849557/1‑150 (MQ=255)
      taCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAAt                                                                                                                                                 <  2:1989360/150‑1 (MQ=255)
      taCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAAt                                                                                                                                                 <  1:1210501/150‑1 (MQ=255)
       aCTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATg                                                                                                                                                <  1:5146312/150‑1 (MQ=255)
          tATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGcac                                                                                                                                             >  1:420073/1‑150 (MQ=255)
          tATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGcac                                                                                                                                             <  2:3495085/150‑1 (MQ=255)
           aTCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGcaca                                                                                                                                            >  2:2072043/1‑150 (MQ=255)
             ctgctgATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCAcatc                                                                                                                                          >  2:4690644/1‑150 (MQ=255)
              tgctgATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATTTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCAcatca                                                                                                                                         <  1:1812801/150‑1 (MQ=255)
              tgctgATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCAcatca                                                                                                                                         <  2:444661/150‑1 (MQ=255)
               gctgATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCAcatcat                                                                                                                                        <  1:1679213/150‑1 (MQ=255)
                ctgATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATg                                                                                                                                       <  1:762247/150‑1 (MQ=255)
                    tGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCttt                                                                                                                                   <  2:4369063/150‑1 (MQ=255)
                    tGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCttt                                                                                                                                   <  1:3185270/150‑1 (MQ=255)
                     gCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTta                                                                                                                                  <  1:1784987/150‑1 (MQ=255)
                     gCGAGAAACACATTACAATGCGGGTATATTCGTATGGGGGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTTATCCTTATATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTta                                                                                                                                  >  1:1075811/1‑150 (MQ=255)
                        agaAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTtatat                                                                                                                               <  1:1186897/150‑1 (MQ=255)
                        agaAACACATTACAAGGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTtatat                                                                                                                               <  1:5024835/150‑1 (MQ=255)
                           aaCACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATNCTTTATATAAc                                                                                                                            <  1:2758064/150‑1 (MQ=255)
                           aaCACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAAc                                                                                                                            >  1:284086/1‑150 (MQ=255)
                           aaCACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAAc                                                                                                                            >  2:4251208/1‑150 (MQ=255)
                           aaCACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCTTTATATAAc                                                                                                                            <  2:1175463/150‑1 (MQ=255)
                              acaTTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAAcctt                                                                                                                         <  2:4649645/150‑1 (MQ=255)
                              acaTTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAAcctt                                                                                                                         >  2:5144542/1‑150 (MQ=255)
                               caTTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAAccttc                                                                                                                        <  1:1384238/150‑1 (MQ=255)
                                aTTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTGTATATAAcattcc                                                                                                                       >  1:3296714/1‑150 (MQ=255)
                                 ttACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAAccttcct                                                                                                                      <  2:2737675/150‑1 (MQ=255)
                                     aaTGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCa                                                                                                                  <  2:2715588/150‑1 (MQ=255)
                                      aTGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCTTTATATAACCTTCCTTGCAc                                                                                                                 <  1:4372131/150‑1 (MQ=255)
                                      aTGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCTTTATATAACCTTCCTTGCAc                                                                                                                 <  1:3931136/150‑1 (MQ=255)
                                       tGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACt                                                                                                                <  2:232260/150‑1 (MQ=255)
                                         cGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTc                                                                                                              >  1:562547/1‑150 (MQ=255)
                                            gTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCctc                                                                                                           >  2:3646072/1‑150 (MQ=255)
                                            gTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCctc                                                                                                           <  2:1707105/150‑1 (MQ=255)
                                            gTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCctc                                                                                                           >  1:619267/1‑150 (MQ=255)
                                             tataTTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCctct                                                                                                          <  2:926295/150‑1 (MQ=255)
                                             tataTTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCctct                                                                                                          <  2:5087209/150‑1 (MQ=255)
                                               tatTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTAGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTctgc                                                                                                        >  2:3066929/1‑150 (MQ=255)
                                                    gTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTgggg                                                                                                   >  2:4091369/1‑150 (MQ=255)
                                                    gTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCCGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTgggg                                                                                                   >  1:4543710/1‑150 (MQ=255)
                                                         tgtgtTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGc                                                                                              >  2:2450149/2‑150 (MQ=255)
                                                         ggTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGc                                                                                              >  2:4335867/1‑150 (MQ=255)
                                                           tgtTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCtt                                                                                            <  1:1910841/150‑1 (MQ=255)
                                                            gtTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCttt                                                                                           >  1:1035669/1‑150 (MQ=255)
                                                               aaTGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGAt                                                                                        >  2:2242246/1‑150 (MQ=255)
                                                                 tGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATgg                                                                                      <  2:3308658/150‑1 (MQ=255)
                                                                   ttttCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCt                                                                                    >  1:987296/1‑150 (MQ=255)
                                                                    tttCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCAAGCTTTATATAAAATTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCtt                                                                                   >  2:3590592/1‑150 (MQ=255)
                                                                       cACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCGTACa                                                                                >  1:3045994/1‑150 (MQ=255)
                                                                        aCGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAg                                                                               >  1:2697803/1‑150 (MQ=255)
                                                                         cGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGt                                                                              <  2:790671/150‑1 (MQ=255)
                                                                          gCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTg                                                                             <  1:44710/150‑1 (MQ=255)
                                                                           ccTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTgg                                                                            <  1:3418476/150‑1 (MQ=255)
                                                                            cTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGa                                                                           >  1:2428987/1‑150 (MQ=255)
                                                                               aTTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGaaaa                                                                        >  2:3686064/1‑150 (MQ=255)
                                                                               aTTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGaaaa                                                                        >  1:2164503/1‑150 (MQ=255)
                                                                                ttATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAAt                                                                       <  1:3708637/150‑1 (MQ=255)
                                                                                   taAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGt                                                                    >  2:5059871/1‑150 (MQ=255)
                                                                                        cacCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGt                                                               >  1:230796/1‑150 (MQ=255)
                                                                                               ataGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTg                                                        <  1:2577115/150‑1 (MQ=255)
                                                                                               ataGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTg                                                        >  2:1346579/1‑150 (MQ=255)
                                                                                                  gTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACATTGGAAAATCGTGGAGTTTCCTTGGtt                                                     <  1:1428091/150‑1 (MQ=255)
                                                                                                    gATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTc                                                   >  2:4466172/1‑150 (MQ=255)
                                                                                                     aTCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCa                                                  >  2:3965427/1‑150 (MQ=255)
                                                                                                     aTCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTAGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCa                                                  >  2:3892713/1‑150 (MQ=255)
                                                                                                        cTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAAtg                                               >  1:4834731/1‑150 (MQ=255)
                                                                                                        cTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAAtg                                               >  2:4185407/1‑150 (MQ=255)
                                                                                                        cTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAAtg                                               >  1:431329/1‑150 (MQ=255)
                                                                                                         ttCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAAtgt                                              >  2:152024/1‑150 (MQ=255)
                                                                                                         ttCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAAtgt                                              <  1:2488576/150‑1 (MQ=255)
                                                                                                         ttCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAAtgt                                              <  1:292020/150‑1 (MQ=255)
                                                                                                            tATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGtc                                           >  1:200318/1‑150 (MQ=255)
                                                                                                            tATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGtc                                           <  1:3522638/150‑1 (MQ=255)
                                                                                                             aTTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGtct                                          >  1:4064659/1‑150 (MQ=255)
                                                                                                              ttGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGtctc                                         >  2:698706/1‑150 (MQ=255)
                                                                                                                  gccaccgccaccTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCtt                                     <  2:5121672/150‑1 (MQ=255)
                                                                                                                   ccaccgccaccTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCttt                                    >  1:4871593/1‑150 (MQ=255)
                                                                                                                      ccgccaccTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGtt                                 <  2:5246339/150‑1 (MQ=255)
                                                                                                                            cctgcctgGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCa                           >  2:376019/1‑150 (MQ=255)
                                                                                                                             ctgcctgGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCaa                          >  1:1962080/1‑150 (MQ=255)
                                                                                                                                cctgGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACa                       >  1:4634436/1‑150 (MQ=255)
                                                                                                                                 ctgGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAg                      >  1:4793838/1‑150 (MQ=255)
                                                                                                                                 ctgGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAg                      >  1:2179683/1‑150 (MQ=255)
                                                                                                                                   ggCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCg                    <  1:378931/150‑1 (MQ=255)
                                                                                                                                        aTGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAAtt               >  2:4761253/1‑150 (MQ=255)
                                                                                                                                           aaaGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTaaa            <  1:3963669/150‑1 (MQ=255)
                                                                                                                                           aaaGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTaaa            >  1:4145198/1‑150 (MQ=255)
                                                                                                                                            aaGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAg           <  1:3685891/150‑1 (MQ=255)
                                                                                                                                               cGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCAc        >  1:1535352/1‑150 (MQ=255)
                                                                                                                                               cGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCAc        >  2:5329640/1‑150 (MQ=255)
                                                                                                                                                 gCTTAGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTg      <  2:110109/150‑1 (MQ=255)
                                                                                                                                                   ttCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGAt    >  2:4140105/1‑150 (MQ=255)
                                                                                                                                                   ttCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGAt    <  2:3436713/150‑1 (MQ=255)
                                                                                                                                                    tCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATa   <  2:1847908/150‑1 (MQ=255)
                                                                                                                                                     cGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATaa  <  2:655982/150‑1 (MQ=255)
                                                                                                                                                     cGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATaa  >  1:3384637/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACCTGCCTGGCAAATGAAAGCGGCTTCGCTAATGCACATCATGCTTTATATAACCTTCCTTGCACTTCCTCTGCTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAA  >  CP009273/2035953‑2036251

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: