Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,756,017 A→C 6.3% intergenic (+639/‑142) smpB → / → intA tmRNA‑binding trans‑translation protein/CP4‑57 prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,756,0170AC 6.3% 337.5 / 7.2 115intergenic (+639/‑142)smpB/intAtmRNA‑binding trans‑translation protein/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base A (47/60);  new base C (3/5);  total (50/65)
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 118 total aligned reads displayed.

GCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAG  >  NC_000913/2755870‑2756166
                                                                                                                                                   |                                                                                                                                                     
 cAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCATAAGAGCCCGCACGGCGCAAAcc                                                                                                                                                    >  2:1621004/1‑150 (MQ=255)
 cAAGCGAATGTAAACACTGACTAACCATGTAGTAACGAGTATGAAGCAATTTCGGAAGCGGTTTAAACTCCAGCAAGCAACACCAAAATACTCCATCGGTAATTACAAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGcc                                                                                                                                                    >  2:214873/1‑150 (MQ=255)
   aGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCt                                                                                                                                                  >  1:1734296/1‑150 (MQ=255)
     cGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGc                                                                                                                                                >  2:110781/1‑150 (MQ=255)
     cGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGc                                                                                                                                                >  1:1622159/1‑150 (MQ=255)
     cGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGc                                                                                                                                                >  2:1551054/1‑150 (MQ=255)
       aaTGTAAAGACTGACTAAGCTTGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCgg                                                                                                                                              <  1:233795/150‑1 (MQ=255)
        aTGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCggg                                                                                                                                             <  1:634788/150‑1 (MQ=255)
         tGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCAGGGCTTGCGCCGTGCGGGc                                                                                                                                            <  2:464967/150‑1 (MQ=255)
         tGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCAGGGCTTGCGCCGTGCGGGc                                                                                                                                            <  2:1667059/150‑1 (MQ=255)
          gTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCATCTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGcg                                                                                                                                           >  1:1298737/1‑149 (MQ=255)
           tAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCACCTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCtt                                                                                                                                          >  1:456675/1‑150 (MQ=255)
           tAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCAGGGCTTGCGCCGTGCGGGCtc                                                                                                                                          <  2:558579/150‑2 (MQ=255)
           tAAAGACTGACTAAGCATGTAGTACCGAGGAAGTAGGAATTTCGGACGCGGGTTCAAATCCCGCAAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCtt                                                                                                                                          <  2:1539965/150‑1 (MQ=255)
            aaaGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCttt                                                                                                                                         <  1:846002/150‑1 (MQ=255)
            aaaGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCttt                                                                                                                                         <  1:724182/150‑1 (MQ=255)
             aaGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCtttt                                                                                                                                        <  2:296799/150‑1 (MQ=255)
              aGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCttttt                                                                                                                                       <  2:1403356/150‑1 (MQ=255)
               gactgactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCATCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCtttttt                                                                                                                                      <  2:1620366/150‑1 (MQ=255)
                actgcctAAGCAAGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAAAGCCCGCAGGGCTTGCGCCGTGCGGGCtcttagg                                                                                                                                     <  1:1826693/150‑4 (MQ=38)
                actgactAAGCATGTTGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCCTCCGATGAAGTCCTAAGCGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTtg                                                                                                                                     <  1:911907/150‑1 (MQ=255)
                actgactAAGCATGTAGAACCGAGGATGTAGGAATTTCGGACGCTGGATCAACTCCCGCAAGCTCCACCAAAATTCTCCATCGGTGATAACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCACGCCCAGCGGGCTTTTTtg                                                                                                                                     >  1:285645/1‑150 (MQ=255)
                  tgactACGCCTGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTtgtg                                                                                                                                   <  2:356499/150‑1 (MQ=255)
                   gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCCGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATACGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGc                                                                                                                                  >  1:647031/1‑150 (MQ=255)
                   gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGc                                                                                                                                  >  1:15233/1‑150 (MQ=255)
                   gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGc                                                                                                                                  <  2:1439539/150‑1 (MQ=255)
                   gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGc                                                                                                                                  >  2:128037/1‑150 (MQ=255)
                   gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGc                                                                                                                                  <  1:230384/150‑1 (MQ=255)
                   gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCACCATCGGTGATTACCAGAGGCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGc                                                                                                                                  <  1:680060/150‑1 (MQ=255)
                       aaGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTc                                                                                                                              <  2:1566147/150‑1 (MQ=255)
                       aaGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCGTGTTTGTGCCCTc                                                                                                                              >  1:239398/1‑150 (MQ=255)
                        aGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGAGTACCAGAGTCATCCGATGAAGTCCTAAGAGCCAGCACGGCGCAAGCCCTGCGGGCTTTTTTGGGCCCGCa                                                                                                                             >  1:244601/1‑150 (MQ=255)
                         gCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCaa                                                                                                                            >  1:733184/1‑150 (MQ=255)
                           aTGTAGTACCGAGGATGTAGGAATTTCGGACGCGTGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAAtt                                                                                                                          <  2:1630611/150‑1 (MQ=255)
                            tGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAAttt                                                                                                                         >  2:406100/1‑150 (MQ=255)
                             gtagtaCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTg                                                                                                                        >  1:1112282/1‑150 (MQ=255)
                             gtagtaCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTg                                                                                                                        <  2:1255736/150‑1 (MQ=255)
                              tagtaCCGAGGATGTAGGAATTTCGGACGCAGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGTGGGCTTTTTTGTGCCCTCAATTTGt                                                                                                                       <  1:569192/150‑1 (MQ=255)
                               agtaCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTc                                                                                                                      <  1:874056/150‑1 (MQ=255)
                               agtaCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTc                                                                                                                      <  1:859355/150‑1 (MQ=255)
                               agtaCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTc                                                                                                                      <  2:1611556/150‑1 (MQ=255)
                                gtaCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTcc                                                                                                                     >  1:1683132/1‑150 (MQ=255)
                                 taCCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTccc                                                                                                                    >  1:1199933/1‑150 (MQ=255)
                                            tAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGCGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAAGTTGTCCCGCGAAGGCCgc                                                                                                         >  2:1609425/1‑149 (MQ=255)
                                            tAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCAACCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGACCTCAATTTGTCCCGCGAAGTCCGa                                                                                                         >  2:1199550/1‑150 (MQ=255)
                                              ggAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAag                                                                                                       >  2:542073/1‑150 (MQ=255)
                                              ggAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAag                                                                                                       >  1:1225337/1‑150 (MQ=255)
                                                aaTTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAagag                                                                                                     <  1:1188249/150‑1 (MQ=255)
                                                  tttCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAagagca                                                                                                   >  2:123007/1‑148 (MQ=255)
                                                       gACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAAt                                                                                              >  1:992272/1‑150 (MQ=255)
                                                         cgcgGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTa                                                                                            >  1:1738114/1‑150 (MQ=255)
                                                          gcgGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGATAACTAATTaa                                                                                           <  1:840501/150‑1 (MQ=255)
                                                                tCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGa                                                                                     <  2:1616935/150‑1 (MQ=255)
                                                                 cAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGaa                                                                                    >  2:288221/1‑150 (MQ=255)
                                                                  aaCTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAAc                                                                                   <  1:1303445/150‑1 (MQ=255)
                                                                  aaCTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAAc                                                                                   <  1:110781/150‑1 (MQ=255)
                                                                    cTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAACTCCGACCCt                                                                                 >  1:1447327/1‑150 (MQ=255)
                                                                     tCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCtt                                                                                >  1:1043944/1‑150 (MQ=255)
                                                                      cccGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAGCCttt                                                                               >  2:1429575/1‑150 (MQ=255)
                                                                       ccGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCtttt                                                                              <  1:1611533/150‑1 (MQ=255)
                                                                              cTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCATAAGAGCCCGCACGGCGCAAGCCCTGCGCGCTTTTTTGTGCCCTCAATTTGTCACGCGAAGTCCTAAGAGAACAAATTTAATCCGAACCTTTTAGGCCCa                                                                       >  2:397198/1‑150 (MQ=255)
                                                                                 caccaAAATTCTCCATCGGTGATTACCAGAGGCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTg                                                                    <  1:1005465/150‑1 (MQ=255)
                                                                                  accaAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGa                                                                   >  2:1455669/1‑150 (MQ=255)
                                                                                    caAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCGCGGCGCAAGCCCTGCGGGCATTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACGATTTATGCCCATTGAtt                                                                 >  1:1066108/1‑149 (MQ=255)
                                                                                     aaaaTTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAg                                                                <  1:1551054/150‑1 (MQ=255)
                                                                                      aaaTTCTCCATCGGTGATTACCAGAGGCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAgg                                                               <  1:1188659/150‑1 (MQ=255)
                                                                                       aaTTCTCCATCGGTGATTACCAGAGGCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGc                                                              <  2:536213/150‑1 (MQ=255)
                                                                                           ctcCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCaa                                                          >  2:574971/1‑150 (MQ=255)
                                                                                             ccATCGGTGATTACCAGAGTCATCCTATGAAGTCCTAAGAGACCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACg                                                        <  1:871880/150‑1 (MQ=255)
                                                                                               aTCGGTGATTACCAGAGTCATCCGAAGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCACTCAATTTGTCCAGCGAAGTCCGAAGAGAACTAATTAAATCCGCACCTTTTAGGCCCATTGATAGGCACAACGaa                                                      >  1:1216885/1‑150 (MQ=255)
                                                                                                     gATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGctc                                                <  1:934919/150‑1 (MQ=255)
                                                                                                      aTTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGctct                                               <  2:1465495/150‑1 (MQ=255)
                                                                                                       ttACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTa                                              <  2:1734296/150‑1 (MQ=255)
                                                                                                          ccAGAGTCATCCGATGAAGTCCTAAGAGCCCGCCCGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTg                                           <  1:159062/150‑1 (MQ=255)
                                                                                                          ccAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTAttt                                           >  2:1527199/1‑149 (MQ=255)
                                                                                                            agagTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGtt                                         <  1:542073/150‑1 (MQ=255)
                                                                                                              agTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGCCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTa                                       <  2:1622159/150‑1 (MQ=255)
                                                                                                                tCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACg                                     <  1:967674/150‑1 (MQ=255)
                                                                                                                tCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACg                                     <  2:247333/150‑1 (MQ=255)
                                                                                                                  aTCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGtt                                   >  1:1255904/1‑150 (MQ=255)
                                                                                                                  aTCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGtt                                   >  2:450644/1‑150 (MQ=255)
                                                                                                                    ccGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTgg                                 <  2:48642/150‑1 (MQ=255)
                                                                                                                     cGATGAATTCCTAAGATCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCGTCAATTTGTCCCGCGAAGTCCGAAGATAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTggg                                <  2:865087/150‑1 (MQ=255)
                                                                                                                     cGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTggg                                <  1:214873/150‑1 (MQ=255)
                                                                                                                     cGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTggg                                >  2:689576/1‑150 (MQ=255)
                                                                                                                       aTGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGcc                              <  2:1179352/150‑1 (MQ=255)
                                                                                                                             tCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACg                        >  1:1677725/1‑150 (MQ=255)
                                                                                                                                tAAGAGCCCGCACGGCGCAAGCCATGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAg                     <  2:1692626/150‑1 (MQ=255)
                                                                                                                                 aaGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAgg                    >  1:514963/1‑150 (MQ=255)
                                                                                                                                   aagCCCGCAGGGCTTGCGCCGTGCGGGGTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTGGGGCCTTAACGCAGGga                  >  2:932491/2‑150 (MQ=255)
                                                                                                                                      cccGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGCACTAAGTAACTCCGAACCTTTTAGGCCCATTGATAGGCCCACCGAAAAGCTCGTTTGTTTACGTTGGGCCTAAACGCAGGgagcc               >  1:631794/1‑148 (MQ=255)
                                                                                                                                       ccGCACGGCGCACGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACt              >  2:954794/1‑150 (MQ=255)
                                                                                                                                       ccGCACCGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTACCGCGAAGTCCGAAGAGAACTAATTTAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAACGCTCTATTGTTTACTTTGGGCCTAAACGCAGGGAGACt              >  1:1294628/1‑150 (MQ=255)
                                                                                                                                           aCGGCGCACCCCCGGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTTTTTACGTTGGGCCTAAACGCAGGGAGACTcccc          <  1:906409/150‑1 (MQ=255)
                                                                                                                                             ggCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCAt        >  1:615791/1‑150 (MQ=255)
                                                                                                                                              gcgcAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGACCCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATg       <  1:1429575/150‑1 (MQ=255)
                                                                                                                                              gcgcAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATg       <  1:1005800/150‑1 (MQ=255)
                                                                                                                                                   aGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAg  <  2:1738114/150‑1 (MQ=255)
                                                                                                                                                   aGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAg  <  2:1243079/150‑1 (MQ=255)
                                                                                                                                                   aGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAg  <  1:1612887/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAG  >  NC_000913/2755870‑2756166

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: