Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,432,031 G→A 85.9% intergenic (‑38/‑168) msrA ← / → tamA methionine sulfoxide reductase A/translocation and assembly module for autotransporter export, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,432,0310GA85.9% 579.0 / 80.2 213intergenic (‑38/‑168)msrA/tamAmethionine sulfoxide reductase A/translocation and assembly module for autotransporter export, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (7/23);  new base A (71/112);  total (78/135)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

CATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGA  >  CP009273/4431886‑4432175
                                                                                                                                                 |                                                                                                                                                
cATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGAtt                                                                                                                                                >  2:2149781/1‑148 (MQ=255)
 aTTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTg                                                                                                                                               <  2:860711/148‑1 (MQ=255)
 aTTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTg                                                                                                                                               >  2:892815/1‑148 (MQ=255)
 aTTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTg                                                                                                                                               >  1:864193/1‑148 (MQ=255)
  ttGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGt                                                                                                                                              >  1:1822248/1‑148 (MQ=255)
    gAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTaa                                                                                                                                            <  1:2554860/148‑1 (MQ=255)
        gACCGTTGACCGCATGGAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGTCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAag                                                                                                                                        <  1:2549049/148‑1 (MQ=255)
        gACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAag                                                                                                                                        <  1:102311/148‑1 (MQ=255)
         aCCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAaga                                                                                                                                       >  1:2526952/1‑148 (MQ=255)
         aCCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAaga                                                                                                                                       <  1:3137313/148‑1 (MQ=255)
          ccGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAagag                                                                                                                                      <  1:671061/148‑1 (MQ=255)
              tGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGt                                                                                                                                  >  2:101971/1‑148 (MQ=255)
              tGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGt                                                                                                                                  >  2:767930/1‑148 (MQ=255)
                aCCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTAt                                                                                                                                >  1:267707/1‑148 (MQ=255)
                aCCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTAt                                                                                                                                >  1:267434/1‑148 (MQ=255)
                      tGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGcac                                                                                                                          >  2:746431/1‑148 (MQ=255)
                            gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATa                                                                                                                    <  1:2921790/148‑1 (MQ=255)
                             tGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATaa                                                                                                                   >  1:977901/1‑148 (MQ=255)
                               gCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAacc                                                                                                                 <  1:2905092/148‑1 (MQ=255)
                               gCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAacc                                                                                                                 <  2:231660/148‑1 (MQ=255)
                                cTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAacca                                                                                                                <  1:419967/148‑1 (MQ=255)
                                cTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAacca                                                                                                                >  2:415858/1‑148 (MQ=255)
                                       cATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGaa                                                                                                         <  2:2579594/148‑1 (MQ=255)
                                       cATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGaa                                                                                                         <  2:2884481/148‑1 (MQ=255)
                                          cGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAg                                                                                                      >  1:221415/1‑148 (MQ=255)
                                          cGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAg                                                                                                      <  2:2517274/148‑1 (MQ=255)
                                          cGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAg                                                                                                      >  1:2082532/1‑148 (MQ=255)
                                           ggggTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAgg                                                                                                     <  1:1661430/148‑1 (MQ=255)
                                            gggTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACGAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGt                                                                                                    >  1:1110610/1‑148 (MQ=255)
                                             ggTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTc                                                                                                   <  2:940174/148‑1 (MQ=255)
                                              gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCt                                                                                                  <  2:2611440/148‑1 (MQ=255)
                                                gtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTaa                                                                                                <  1:1973171/148‑1 (MQ=255)
                                                 ttACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTaaa                                                                                               <  2:1448396/148‑1 (MQ=255)
                                                      tCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTaaa                                                                                          <  2:125075/148‑1 (MQ=255)
                                                       ccAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTaaaa                                                                                         >  2:749069/1‑148 (MQ=255)
                                                       ccAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTaaaa                                                                                         >  2:1427309/1‑148 (MQ=255)
                                                       ccAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTaaaa                                                                                         <  1:2684195/148‑1 (MQ=255)
                                                        caggcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAAc                                                                                        <  2:2541280/148‑1 (MQ=255)
                                                        caggcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTaaaaa                                                                                        <  2:1802925/148‑2 (MQ=255)
                                                          ggcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACAt                                                                                      <  2:474909/148‑1 (MQ=255)
                                                             aggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAg                                                                                   <  2:669024/148‑1 (MQ=255)
                                                                gCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGaaa                                                                                <  2:399242/148‑1 (MQ=255)
                                                                 cATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTACAGTAAAACATAAGAAat                                                                               <  1:1340421/148‑1 (MQ=255)
                                                                 cATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAat                                                                               <  2:1914972/148‑1 (MQ=255)
                                                                     ggcggGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATAttt                                                                           <  1:1037860/148‑1 (MQ=255)
                                                                     ggcggGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATAttt                                                                           <  2:221255/148‑1 (MQ=255)
                                                                      gcggGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTg                                                                          <  2:313735/148‑1 (MQ=255)
                                                                          gggAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGt                                                                      <  2:282702/148‑1 (MQ=255)
                                                                          gggAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGt                                                                      >  1:1257966/1‑148 (MQ=255)
                                                                           ggAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGtt                                                                     <  2:2985751/148‑1 (MQ=255)
                                                                            gAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGttt                                                                    <  1:1427148/148‑1 (MQ=255)
                                                                            gAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGttt                                                                    <  1:2882930/148‑1 (MQ=255)
                                                                              aaCCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAg                                                                  >  2:1197635/1‑148 (MQ=255)
                                                                                ccAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTc                                                                <  1:1810401/148‑1 (MQ=255)
                                                                                   gATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGc                                                             >  2:534248/1‑148 (MQ=255)
                                                                                   gATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAAGATTTGGGGTTTAGTCTGc                                                             >  2:2009194/1‑148 (MQ=255)
                                                                                    aTGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCt                                                            >  2:1748765/1‑148 (MQ=255)
                                                                                     tGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCtt                                                           <  2:359028/148‑1 (MQ=255)
                                                                                     tGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCtt                                                           <  1:2298517/148‑1 (MQ=255)
                                                                                      gCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCttt                                                          >  1:2396233/1‑148 (MQ=255)
                                                                                       cTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCtttt                                                         <  1:3074944/148‑1 (MQ=255)
                                                                                       cTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCtttt                                                         <  1:1708900/148‑1 (MQ=255)
                                                                                          tttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTaa                                                      >  1:2155726/1‑148 (MQ=255)
                                                                                           ttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAAt                                                     <  1:1884174/148‑1 (MQ=255)
                                                                                            tATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATc                                                    <  2:113741/148‑1 (MQ=255)
                                                                                               cAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCat                                                 >  1:2762670/1‑148 (MQ=255)
                                                                                               cAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCat                                                 <  1:358304/148‑1 (MQ=255)
                                                                                                  aTAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATAtt                                              <  2:1191647/148‑1 (MQ=255)
                                                                                                  aTAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATAtt                                              <  2:469769/148‑1 (MQ=255)
                                                                                                  aTAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATAtt                                              <  2:1191786/148‑1 (MQ=255)
                                                                                                   tAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTa                                             >  2:2925989/1‑148 (MQ=255)
                                                                                                   tAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTa                                             >  1:1552800/1‑148 (MQ=255)
                                                                                                        tCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGa                                        <  1:1076609/148‑1 (MQ=255)
                                                                                                        tCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGa                                        <  1:892815/148‑1 (MQ=255)
                                                                                                         cATGGTGTCGCTCTCCCGAAAATCGGTCATTTGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATTACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGAt                                       <  1:759873/148‑1 (MQ=255)
                                                                                                           tGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGAttt                                     <  2:1331680/148‑1 (MQ=255)
                                                                                                           tGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGAttt                                     <  1:1580990/148‑1 (MQ=255)
                                                                                                           tGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGAttt                                     <  1:1014455/148‑1 (MQ=255)
                                                                                                               gtCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGt                                 <  1:2128457/148‑1 (MQ=255)
                                                                                                                tCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGtt                                <  2:1730272/148‑1 (MQ=255)
                                                                                                                 cGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTa                               >  2:3048776/1‑148 (MQ=255)
                                                                                                                 cGCTCTCCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTa                               >  1:2025712/1‑148 (MQ=255)
                                                                                                                    tctcCCGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGc                            >  1:570250/1‑148 (MQ=255)
                                                                                                                       cccGAAAATCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGt                         <  2:597972/148‑1 (MQ=255)
                                                                                                                             aaTCGGTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACg                   <  1:746431/148‑1 (MQ=255)
                                                                                                                                  gTCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGtt              >  2:572490/1‑148 (MQ=255)
                                                                                                                                   tCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTc             <  1:3205986/148‑1 (MQ=255)
                                                                                                                                   tCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTc             <  1:1185841/148‑1 (MQ=255)
                                                                                                                                   tCATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACAAAATAACCACCGTGAATAGGACAAAAGTAAAACATAAGAAATATTTGGGGTTAAGTCTGCTTTTTAATCCATATTACTGTATTTTTGTTAAGCCGTTTAACGTCGTTc             >  1:3249719/1‑148 (MQ=255)
                                                                                                                                    cATTGGGGTTAAGATTTTAACAAGAGTCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTcc            <  1:39614/148‑1 (MQ=255)
                                                                                                                                    cATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTcc            <  1:3102461/148‑1 (MQ=255)
                                                                                                                                    cATTGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACAGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCCTATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTcc            <  1:2127191/148‑1 (MQ=255)
                                                                                                                                     aTTGGTGTTAAGATTGTAACAAGATTCGTATTTGCACTCTATAACCACCGTGAATAGGTCTAAAGTAAAACATATGAAATATTTTGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCa           <  2:188999/148‑1 (MQ=255)
                                                                                                                                      ttGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAg          <  2:32349/148‑1 (MQ=255)
                                                                                                                                      ttGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAg          <  2:3271044/148‑1 (MQ=255)
                                                                                                                                      ttGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGACGTTTAACGGCGTTCCAg          >  2:2392922/1‑148 (MQ=255)
                                                                                                                                       tGGGGTTAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAgg         <  1:918708/148‑1 (MQ=255)
                                                                                                                                             tAAGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAgg   >  2:1218956/1‑148 (MQ=255)
                                                                                                                                              aaGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGa  >  2:310166/1‑148 (MQ=255)
                                                                                                                                              aaGATTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGa  >  2:2578601/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGA  >  CP009273/4431886‑4432175

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: