Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,838,201 C→A 6.9% E310* (GAA→TAA)  ynjI ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,838,2010CA 6.9% 1140.6 / 37.0 295E310* (GAA→TAA) ynjIinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (133/141);  new base A (10/11);  total (143/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAA  >  CP009273/1838052‑1838350
                                                                                                                                                     |                                                                                                                                                     
gggCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTc                                                                                                                                                       >  2:2970974/1‑150 (MQ=255)
gggCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTc                                                                                                                                                       <  1:1890411/150‑1 (MQ=255)
gggCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTc                                                                                                                                                       >  1:2299808/1‑150 (MQ=255)
 ggCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCa                                                                                                                                                      <  2:3740200/150‑1 (MQ=255)
  gCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAg                                                                                                                                                     <  2:2770384/150‑1 (MQ=255)
    aaGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                                                                                                                                   >  2:2121597/1‑150 (MQ=255)
     agcgcagcgcCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                                                                                                                                  >  2:3969624/1‑150 (MQ=255)
           gcgcCCCCGGCACAATGTCAGGCGTTTATGCGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTGCAAGTTCAGCATCCGTTa                                                                                                                                            <  2:3156859/150‑1 (MQ=255)
           gcgcCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTa                                                                                                                                            <  2:298024/150‑1 (MQ=255)
             gcCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAAc                                                                                                                                          >  2:3100396/1‑150 (MQ=255)
               ccccGGCAAAATTTAAGGCGTTTATGAGTATTTAACGGAAGATTCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACtt                                                                                                                                        >  2:1762875/1‑150 (MQ=255)
                   ggCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTg                                                                                                                                    >  1:3792474/1‑150 (MQ=255)
                      aaaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGa                                                                                                                                 >  1:3614423/1‑150 (MQ=255)
                       aaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGAt                                                                                                                                <  2:650359/150‑1 (MQ=255)
                        aaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATa                                                                                                                               <  1:2078910/150‑1 (MQ=255)
                           ttCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagca                                                                                                                            <  1:3986168/150‑1 (MQ=255)
                            tCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcag                                                                                                                           <  1:2126280/150‑1 (MQ=255)
                            tCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcag                                                                                                                           <  2:2257224/150‑1 (MQ=255)
                              aGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAgcagca                                                                                                                         >  2:429892/1‑150 (MQ=255)
                               ggCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGcagcag                                                                                                                        <  2:3374390/150‑1 (MQ=255)
                                 cGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGat                                                                                                                      <  1:1787692/150‑1 (MQ=255)
                                  gTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                                                                                                                     >  1:65823/1‑150 (MQ=255)
                                  gTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                                                                                                                     >  2:2416811/1‑150 (MQ=255)
                                  gTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                                                                                                                     >  2:268756/1‑150 (MQ=255)
                                    ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg                                                                                                                   <  1:2978494/150‑1 (MQ=255)
                                     tATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGc                                                                                                                  <  2:3383115/150‑1 (MQ=255)
                                         aGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGt                                                                                                              >  2:1056744/1‑150 (MQ=255)
                                          gTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTa                                                                                                             >  1:2519307/1‑150 (MQ=255)
                                          gTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTa                                                                                                             >  1:495599/1‑150 (MQ=255)
                                           tATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaa                                                                                                            >  1:3097329/1‑150 (MQ=255)
                                            aTTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaaa                                                                                                           <  1:1155607/150‑1 (MQ=255)
                                             tttAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAg                                                                                                          <  1:3402803/150‑1 (MQ=255)
                                               tAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaa                                                                                                        <  2:2122291/150‑1 (MQ=255)
                                                aaCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaaa                                                                                                       <  1:2263869/150‑1 (MQ=255)
                                                 aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAt                                                                                                      <  2:3226200/150‑1 (MQ=255)
                                                 aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAt                                                                                                      <  2:3816694/150‑1 (MQ=255)
                                                 aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAt                                                                                                      >  1:2334218/1‑150 (MQ=255)
                                                  cGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATc                                                                                                     >  1:3739472/1‑150 (MQ=255)
                                                   ggATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCATCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATcc                                                                                                    >  2:2496172/1‑150 (MQ=255)
                                                    gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATccc                                                                                                   >  2:1896595/1‑150 (MQ=255)
                                                    gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATccc                                                                                                   >  2:3441056/1‑150 (MQ=255)
                                                     atgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCa                                                                                                  <  2:3650843/150‑1 (MQ=255)
                                                      tgatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCAt                                                                                                 <  2:1608117/150‑1 (MQ=255)
                                                       gatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCAtt                                                                                                >  2:1724815/1‑150 (MQ=255)
                                                          gCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGa                                                                                             <  2:1004130/150‑1 (MQ=255)
                                                           cTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGAc                                                                                            >  2:2568095/1‑150 (MQ=255)
                                                            tCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACt                                                                                           >  1:1926968/1‑150 (MQ=255)
                                                                cACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttt                                                                                       >  2:2357968/1‑150 (MQ=255)
                                                                 aCGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAtttt                                                                                      <  1:843040/150‑1 (MQ=255)
                                                                   ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAtttttt                                                                                    >  1:4095832/1‑150 (MQ=255)
                                                                   ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAtttttt                                                                                    <  1:1509798/150‑1 (MQ=255)
                                                                           tCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtc                                                                            <  1:977301/150‑1 (MQ=255)
                                                                           tCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtc                                                                            <  2:80774/150‑1 (MQ=255)
                                                                             ttATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtctt                                                                          >  1:3459548/1‑150 (MQ=255)
                                                                             ttATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtctt                                                                          >  1:1861652/1‑150 (MQ=255)
                                                                              tATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttc                                                                         >  1:199481/1‑150 (MQ=255)
                                                                                tGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATcttctt                                                                       >  1:884629/1‑150 (MQ=255)
                                                                                tGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTGGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATcttctt                                                                       >  1:3353711/1‑150 (MQ=255)
                                                                                 gggCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCttcttc                                                                      >  2:2064951/1‑150 (MQ=255)
                                                                                         ggCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGa                                                              >  2:1237176/1‑150 (MQ=255)
                                                                                         ggCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGa                                                              >  1:843680/1‑150 (MQ=255)
                                                                                           cATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaac                                                            >  2:1115968/1‑150 (MQ=255)
                                                                                            aTTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaaca                                                           >  2:4112449/1‑150 (MQ=255)
                                                                                            aTTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaaca                                                           >  1:452050/1‑150 (MQ=255)
                                                                                               tCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACt                                                        <  2:1523610/150‑1 (MQ=255)
                                                                                                 ttACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCg                                                      <  1:3893320/150‑1 (MQ=255)
                                                                                                   aCTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg                                                    >  2:143370/1‑150 (MQ=255)
                                                                                                    cTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTgc                                                   >  1:3217571/1‑150 (MQ=255)
                                                                                                    cTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTgc                                                   <  2:1172039/150‑1 (MQ=255)
                                                                                                       tAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGcc                                                >  2:2744008/1‑150 (MQ=255)
                                                                                                          tGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCttt                                             <  1:3612453/150‑1 (MQ=255)
                                                                                                           gCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTc                                            <  1:3265408/150‑1 (MQ=255)
                                                                                                           gCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTc                                            >  2:3110529/1‑150 (MQ=255)
                                                                                                              cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCgag                                         <  2:1923709/150‑1 (MQ=255)
                                                                                                               ccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCgaga                                        >  2:765862/1‑150 (MQ=255)
                                                                                                               ccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCgaga                                        >  1:3049048/1‑150 (MQ=255)
                                                                                                                 aaaaCTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGaag                                      >  2:567941/1‑150 (MQ=255)
                                                                                                                 aaaaCTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGaag                                      >  1:3065480/1‑150 (MQ=255)
                                                                                                                      tGCTTGTCGTAACGATAACACGCTTCAAGTTAAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGaagcaagc                                 >  2:3389867/1‑150 (MQ=255)
                                                                                                                        cTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCAt                               <  2:1607632/150‑1 (MQ=255)
                                                                                                                          tGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATta                             >  1:2681874/1‑150 (MQ=255)
                                                                                                                          tGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATta                             <  1:1301659/150‑1 (MQ=255)
                                                                                                                          tGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATta                             >  1:1012052/1‑150 (MQ=255)
                                                                                                                             cGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATtatat                          <  2:1860255/150‑1 (MQ=255)
                                                                                                                             cGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATtatat                          <  1:124654/150‑1 (MQ=255)
                                                                                                                              gTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATtatata                         <  1:4140655/150‑1 (MQ=255)
                                                                                                                              gTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATtatata                         >  1:3715687/1‑150 (MQ=255)
                                                                                                                              gTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATtatata                         <  1:930668/150‑1 (MQ=255)
                                                                                                                                    aTAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgcca                   <  2:917621/150‑1 (MQ=255)
                                                                                                                                    aTAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgcca                   <  1:224960/150‑1 (MQ=255)
                                                                                                                                     tAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgccag                  >  2:2847640/1‑150 (MQ=255)
                                                                                                                                     tAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgccag                  >  2:3055493/1‑150 (MQ=255)
                                                                                                                                        cacGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgccaggcc               >  2:3504276/1‑150 (MQ=255)
                                                                                                                                          cGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGAGCGCCTTTCGAGAAGCAAGCATTATATAATgacaggccag             >  2:3699991/1‑150 (MQ=255)
                                                                                                                                           gCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGt            >  2:1691132/1‑150 (MQ=255)
                                                                                                                                                 aGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGttcttct      <  1:2579087/150‑1 (MQ=255)
                                                                                                                                                 aGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGttcttct      >  1:1629267/1‑150 (MQ=255)
                                                                                                                                                   ttCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGttcttcttc    <  2:1378345/150‑1 (MQ=255)
                                                                                                                                                    tCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCa   >  1:1658955/1‑150 (MQ=255)
                                                                                                                                                     cAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCaa  >  1:2873190/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCAAGCGCAGCGCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAA  >  CP009273/1838052‑1838350

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: