Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,467,027 C→T R161R (CGC→CGT intS → CPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,467,0270CT98.7% 571.1 / NA 150R161R (CGC→CGTintSCPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (83/65);  total (85/65)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

AACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAA  >  NC_000913/2466883‑2467173
                                                                                                                                                |                                                                                                                                                  
aaCTGAACTTGCCAAAATGTTTGCCGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTagaaga                                                                                                                                              <  1:1085779/151‑1 (MQ=255)
aaCTGAACTTGCAAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTagaaga                                                                                                                                              <  2:1369805/151‑1 (MQ=255)
   tGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATg                                                                                                                                            >  4:557009/1‑150 (MQ=255)
     aaCTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcgg                                                                                                                                         <  4:575165/151‑1 (MQ=255)
     aaCTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcgg                                                                                                                                         >  2:452696/1‑151 (MQ=255)
      aCTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcggc                                                                                                                                        >  3:466915/1‑151 (MQ=255)
       cTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcggcg                                                                                                                                       >  2:1556433/1‑151 (MQ=255)
        ttGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcggcg                                                                                                                                       >  1:1447045/1‑150 (MQ=255)
        ttGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcggcg                                                                                                                                       >  2:564825/1‑150 (MQ=255)
        ttGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcggcg                                                                                                                                       >  2:851086/1‑150 (MQ=255)
         tGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATgcggcg                                                                                                                                       <  1:182305/149‑1 (MQ=255)
           ccAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAgg                                                                                                                                    <  1:1233408/150‑1 (MQ=255)
           ccAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGt                                                                                                                                   <  2:1550470/151‑1 (MQ=255)
           ccAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGt                                                                                                                                   <  3:312755/151‑1 (MQ=255)
           ccAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGt                                                                                                                                   >  4:298376/1‑151 (MQ=255)
            cAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGt                                                                                                                                   <  1:1144667/150‑1 (MQ=255)
             aaaaTGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGttt                                                                                                                                 >  3:556863/1‑151 (MQ=255)
              aaaTGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGtttt                                                                                                                                >  1:243819/1‑151 (MQ=255)
               aaTGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTc                                                                                                                               >  2:1419956/1‑151 (MQ=255)
                aTGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTcc                                                                                                                              >  2:1346065/1‑151 (MQ=255)
                aTGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTcc                                                                                                                              >  4:735/1‑151 (MQ=255)
                  gTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCg                                                                                                                             <  4:166102/150‑1 (MQ=255)
                        cgacgacATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCt                                                                                                                      <  2:413665/151‑1 (MQ=255)
                          acgacATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTAt                                                                                                                    <  2:1617377/151‑1 (MQ=255)
                            gacATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTg                                                                                                                  <  1:1070291/151‑1 (MQ=255)
                             acATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGt                                                                                                                 <  3:8929/151‑1 (MQ=255)
                              cATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTc                                                                                                                <  1:1389297/151‑1 (MQ=255)
                              cATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTc                                                                                                                >  2:541169/1‑151 (MQ=255)
                              cATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTc                                                                                                                >  2:240795/1‑151 (MQ=255)
                               atgttACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTc                                                                                                                >  4:336356/4‑150 (MQ=255)
                               aTTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTc                                                                                                                >  4:515243/1‑150 (MQ=255)
                                  ttACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCAcc                                                                                                             <  2:501490/150‑1 (MQ=255)
                                   tACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCgg                                                                                                           >  2:729759/1‑151 (MQ=255)
                                    aCCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCgg                                                                                                           <  1:1458107/150‑1 (MQ=255)
                                    aCCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCgg                                                                                                           >  1:707761/1‑150 (MQ=255)
                                     ccTATCATTAGCGGCCTTGAAATTGAGGATAATGAGCCGATGCAACTGCTGGAAGTAATCAGCAGGTTTGAAGATCGCGGTGCTATGGAACGAGCCAACAAATCACGTAGAAGATGAGGCGAGGTTTTCCGTTACGGTATTGTAAcc                                                                                                             >  1:1612138/1‑147 (MQ=255)
                                       tATCATTGGAGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGCAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGacg                                                                                                        >  1:168831/1‑148 (MQ=255)
                                          cATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCt                                                                                                    >  1:1364770/1‑151 (MQ=255)
                                             tGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTaa                                                                                                  >  1:725237/1‑150 (MQ=255)
                                             tGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTaa                                                                                                  >  1:670931/1‑150 (MQ=255)
                                              ggcggcCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTaaa                                                                                                 <  1:613003/150‑1 (MQ=255)
                                                 tgcCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAatat                                                                                              >  1:814363/2‑150 (MQ=255)
                                                 ggcATTGAAATTCAGGATATTGAGCCGATGCAACAGCTGGAAGGAATCCGCAGGTTTGAAGATCGCGGTGCAATGAAACGAGCCAACAAAGCAAGTAGAAGATGCGGCGAGGTTTTCCGGTAAGCTATTGACACCGGAAGGGCTAAat                                                                                                >  1:1262359/1‑148 (MQ=255)
                                                  gcCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATaa                                                                                            <  1:175248/151‑1 (MQ=255)
                                                  gcCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATaa                                                                                            <  1:420475/151‑1 (MQ=255)
                                                     ttGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAAccc                                                                                         >  3:476738/1‑151 (MQ=255)
                                                      tGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAAcc                                                                                          <  2:433701/149‑1 (MQ=255)
                                                        aaaTTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGc                                                                                      >  2:294301/1‑151 (MQ=255)
                                                        aaaTTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGc                                                                                      <  2:1541243/151‑1 (MQ=255)
                                                           ttCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCAcc                                                                                   >  2:406960/1‑151 (MQ=255)
                                                              aGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTg                                                                                 >  1:419416/1‑150 (MQ=255)
                                                                   atTGAGCCGCTGCAACTGCTGGAAGTAATCCACAGGGTTGAAGATCGCGGTGCCATGGAACGAGCCACCCAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTGTTGTCAACGGAAGGGATAAATATAACCCGGCACCTGACCTTg                                                                           >  1:1698914/1‑151 (MQ=255)
                                                                   atTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCtt                                                                            <  2:293712/150‑1 (MQ=255)
                                                                      gAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTg                                                                        >  4:252613/1‑151 (MQ=255)
                                                                       aGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTg                                                                        >  1:762247/1‑150 (MQ=255)
                                                                       aGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTg                                                                        >  2:1214212/1‑150 (MQ=255)
                                                                             tGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCAt                                                                 >  1:533547/1‑151 (MQ=255)
                                                                                aaCTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGaa                                                              <  2:125270/151‑1 (MQ=255)
                                                                                aaCTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGa                                                               >  2:547149/1‑150 (MQ=255)
                                                                                aaCTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGAGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGa                                                               >  1:1311728/1‑150 (MQ=255)
                                                                                    gctgGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAgg                                                            >  3:91759/1‑149 (MQ=255)
                                                                                    gctgGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGa                                                          <  2:423463/151‑1 (MQ=255)
                                                                                       ggAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATAc                                                       >  3:170524/1‑151 (MQ=255)
                                                                                        gAAGTAATCCGCAGTTTTGAAGATCGCGGTGCAATGGAGCGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACAGGAAGGGCTAAATATAACCCGGCACCTGCCCTTGCTGACGCCATGAAGGGATGcc                                                      <  1:212358/151‑1 (MQ=255)
                                                                                        gAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATAcc                                                      <  1:1548207/151‑1 (MQ=255)
                                                                                        gAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATAcc                                                      >  3:415585/1‑151 (MQ=255)
                                                                                         aaGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCg                                                     >  2:172204/1‑151 (MQ=255)
                                                                                           gTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGc                                                    >  3:344706/1‑150 (MQ=255)
                                                                                             aaTCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCa                                                   >  1:1626460/1‑149 (MQ=255)
                                                                                              aTCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCa                                                   >  2:492394/1‑148 (MQ=255)
                                                                                                    aGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACt                                          >  2:577189/1‑151 (MQ=255)
                                                                                                     ggTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGGGGCGAGGGTTTTCGTTACGCTAATGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGAAACCGCAAGAAGAACtt                                         >  1:286234/1‑151 (MQ=255)
                                                                                                     ggTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACtt                                         <  1:857482/151‑1 (MQ=255)
                                                                                                     ggTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACtt                                         <  2:707864/151‑1 (MQ=255)
                                                                                                        ttGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTcc                                       <  1:544111/150‑1 (MQ=255)
                                                                                                          gAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGt                                    <  2:765251/151‑1 (MQ=255)
                                                                                                           aaGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGtt                                   <  2:1445072/151‑1 (MQ=255)
                                                                                                            aGATCGCGGGGCAATGCAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCCATTGTCACCGGACGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGCAGAACTTCCCGttt                                  >  1:337264/1‑151 (MQ=255)
                                                                                                             gATCGCGGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTttc                                 >  3:563079/1‑151 (MQ=255)
                                                                                                               tCGCGGTGCAATGGAACGAGCCAACAAATCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTTTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGTGATACCGCAAGGAGAACATCCCGTttctt                               >  2:471985/1‑151 (MQ=255)
                                                                                                                 gcgGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTcc                             <  2:230488/151‑1 (MQ=255)
                                                                                                                  cgGTGCAATGGAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCt                            <  2:1497400/151‑1 (MQ=255)
                                                                                                                           ggAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCa                    <  4:490990/150‑1 (MQ=255)
                                                                                                                           ggAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCa                    <  1:1214068/150‑1 (MQ=255)
                                                                                                                            gAACGAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGa                  <  2:1531982/151‑1 (MQ=255)
                                                                                                                                gAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATccc              <  1:1556274/151‑1 (MQ=255)
                                                                                                                                gAGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATccc              >  1:116311/1‑151 (MQ=255)
                                                                                                                                 aGCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATTAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCg             >  2:363384/1‑151 (MQ=255)
                                                                                                                                  gCCAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCgg            >  2:414871/1‑151 (MQ=255)
                                                                                                                                   ccAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc           <  2:1718983/151‑1 (MQ=255)
                                                                                                                                   ccAACAAAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc           >  2:1132097/1‑151 (MQ=255)
                                                                                                                                    caacaaAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc           >  3:5187/1‑150 (MQ=255)
                                                                                                                                    caacaaAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc           >  1:89065/1‑150 (MQ=255)
                                                                                                                                    caacaaAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc           <  1:760090/150‑1 (MQ=255)
                                                                                                                                     aacaaAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCAt         <  2:616484/151‑1 (MQ=255)
                                                                                                                                      acaaAGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCAtt        <  2:922449/151‑1 (MQ=255)
                                                                                                                                        aaaGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTca      <  4:556843/151‑1 (MQ=255)
                                                                                                                                          aGCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCATGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTcaac    <  1:794335/151‑1 (MQ=255)
                                                                                                                                           gCACGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTcaaca   <  3:508560/151‑1 (MQ=255)
                                                                                                                                             aCGTAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTcaacaa  >  2:1144622/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AACTGAACTTGCCAAAATGTTTGACGACGACATTTTACCTATCATTGGCGGCCTTGAAATTCAGGATATTGAGCCGATGCAACTGCTGGAAGTAATCCGCAGGTTTGAAGATCGCGGTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAA  >  NC_000913/2466883‑2467173

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