Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,321 C→T 15.6% A22V (GCT→GTT)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3210CT15.6% 442.4 / 59.7 179A22V (GCT→GTT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (67/84);  new base T (16/12);  total (83/96)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 193 total aligned reads displayed.

AATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCA  >  NC_000913/2560225‑2560420
                                                                                                |                                                                                                   
aaTAAATTATTCAATATCAAGAATTTAAACCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACTCAAATAAAAGCCCCCGGAAtgctg                                                                                                   <  2:805020/99‑1 (MQ=255)
aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                                  <  2:3267091/100‑1 (MQ=255)
 aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  2:3316554/100‑1 (MQ=255)
 aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  2:3704315/100‑1 (MQ=255)
 aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACCAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                   <  2:2827023/98‑1 (MQ=255)
 aTAAATTATTAAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  2:2720759/100‑1 (MQ=255)
   aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAGTAAAAGCCCCCGGAATGCTGCAAc                                                                                               >  1:3554743/1‑100 (MQ=255)
    aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGTAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                                >  2:2500773/1‑98 (MQ=255)
       tattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                           <  1:1148066/100‑1 (MQ=255)
        attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGcc                                                                                          <  1:3852614/100‑1 (MQ=255)
         ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                         <  1:30767/100‑1 (MQ=255)
          tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGTATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgg                                                                                        <  2:15904/100‑1 (MQ=255)
           cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                       <  1:3137197/100‑1 (MQ=38)
           cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                       <  1:160905/100‑1 (MQ=38)
           cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                       <  2:1032608/100‑1 (MQ=38)
           cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                       <  1:2438152/100‑1 (MQ=255)
           cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                       >  1:1690507/1‑100 (MQ=255)
           cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTAGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                       <  2:616555/100‑1 (MQ=38)
            aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                      <  2:1751561/100‑1 (MQ=38)
            aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                      >  2:1154605/1‑100 (MQ=255)
             atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                     <  2:3016425/100‑1 (MQ=255)
              tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAATATCGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                      <  2:1443692/98‑1 (MQ=37)
              tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                    <  1:3828839/100‑1 (MQ=255)
              tatCAAGAACTTAAATCCATGCTTGCAACGAAAGATCCATAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                    <  2:1378956/100‑1 (MQ=255)
               atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGct                                                                                   <  2:2582786/100‑2 (MQ=38)
               atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                   >  2:1661559/1‑100 (MQ=255)
               atCAAAAACTTAAGTCCGCGTTTGCAACGAAAGATCCAGGGTTGGAAAATAGGATCAACACAAATAAAAGCCCCCTGCATGTTGCAGCATTCCggggg                                                                                     <  2:2393602/98‑1 (MQ=25)
                tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  >  1:2161839/1‑100 (MQ=255)
                tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  >  1:1889903/1‑100 (MQ=255)
                 caacaaCTTCAATCCCTCTTTGCATCTTATGATCCTACGTCCACAAACAGGATCACCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  <  2:4080503/99‑1 (MQ=255)
                     aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                             <  1:1516337/100‑1 (MQ=255)
                       cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                           >  1:804197/1‑100 (MQ=255)
                       cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTGCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                           <  2:1879526/100‑1 (MQ=255)
                        ttcaaTCCATGCTTCCAACAAAAGATCCACAGTTAGCAAATAGGATCACCACAACAACACGACCCAGGAAAGCTGCAACATGCCGGGGGAGGTATgagac                                                                          <  2:3293000/97‑2 (MQ=255)
                        ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                          >  2:3839917/1‑100 (MQ=255)
                          aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                        <  1:1742438/100‑1 (MQ=255)
                           aaTCCATGTTTGCAACGAGAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                       <  2:2079152/100‑1 (MQ=255)
                              ccATGTTTGCAACGTAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                    <  1:2500787/100‑1 (MQ=255)
                              ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                    <  2:3071805/100‑1 (MQ=255)
                              ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                    <  2:3396432/100‑1 (MQ=255)
                               cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACGAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                    >  1:2380008/1‑99 (MQ=255)
                               cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCGACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                   >  1:210624/1‑100 (MQ=255)
                                aTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTa                                                                  <  2:3000525/100‑1 (MQ=255)
                                  gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaca                                                                >  2:2514571/1‑98 (MQ=255)
                                  gTTTGCAACGAAAGATCCAGAGCTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCAGGGGCGGTATGAGAAACAAGCTaa                                                                 >  2:849528/1‑99 (MQ=255)
                                   tttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                               >  1:708555/1‑100 (MQ=255)
                                    ttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                               <  2:90818/99‑1 (MQ=255)
                                     tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAACAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGga                                                             <  2:3137638/100‑1 (MQ=255)
                                        aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                          <  1:2419027/100‑1 (MQ=255)
                                           gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                       >  1:2457343/1‑100 (MQ=255)
                                             aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                     <  2:1933166/100‑1 (MQ=255)
                                             aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                     >  2:499818/1‑100 (MQ=255)
                                              cttccgatctttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  1:1129350/90‑1 (MQ=255)
                                              aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  1:1927834/100‑1 (MQ=255)
                                              aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  2:216789/100‑1 (MQ=255)
                                              aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  1:3031671/100‑1 (MQ=255)
                                               tccgatctgTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                   <  1:2565936/92‑1 (MQ=255)
                                               gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                   <  2:1763998/100‑1 (MQ=255)
                                               gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                   <  1:2529546/100‑1 (MQ=255)
                                                aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                   >  2:3087978/1‑99 (MQ=255)
                                                   cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                               >  1:1380201/1‑100 (MQ=255)
                                                   cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                               <  1:561455/100‑1 (MQ=255)
                                                    agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                              >  1:2416671/1‑100 (MQ=255)
                                                       gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTTCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                           <  2:722613/100‑1 (MQ=255)
                                                       gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                           >  1:3929824/1‑100 (MQ=255)
                                                         tCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                         <  2:769476/100‑1 (MQ=255)
                                                           gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                       >  1:2683777/1‑100 (MQ=255)
                                                             aaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                     >  1:3048785/1‑100 (MQ=255)
                                                                   gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                               >  2:3508798/1‑100 (MQ=255)
                                                                       aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                           >  2:3790298/1‑100 (MQ=255)
                                                                        acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                          >  1:2883643/1‑100 (MQ=255)
                                                                        acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                          <  2:2416657/100‑1 (MQ=255)
                                                                         cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                         <  2:2369845/100‑1 (MQ=255)
                                                                           caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                       >  1:4028872/1‑100 (MQ=255)
                                                                             tcataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                     >  1:249072/8‑100 (MQ=25)
                                                                             aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAAAAGCTAAAGGAGAATTTCTCATGTCTGATATGATAGCAGTAGAACAAGaaa                     >  1:2108570/1‑100 (MQ=255)
                                                                              cataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                    >  1:1587780/7‑100 (MQ=25)
                                                                              aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                    <  2:3769489/100‑1 (MQ=255)
                                                                                taccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTAGTATGCTAGCAGTAGGACAAGaaaaaa                  >  2:2733983/5‑100 (MQ=25)
                                                                                aaaaGCCCCCAGAATGCGGCAACATGACGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAAAACAAGaaa                     <  2:155272/97‑1 (MQ=255)
                                                                                 aaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa                 >  2:1239967/1‑100 (MQ=255)
                                                                                  ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                >  1:1493733/3‑100 (MQ=38)
                                                                                  aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                <  2:3619615/100‑1 (MQ=255)
                                                                                  aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                <  1:1393148/100‑1 (MQ=255)
                                                                                   aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca               >  1:2323754/1‑100 (MQ=255)
                                                                                    gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa              >  1:399479/1‑100 (MQ=38)
                                                                                    gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTACCAGTAGAACAAGAaataaaca               >  2:3982801/1‑99 (MQ=37)
                                                                                     cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAaaaaacaa              >  1:1464475/1‑99 (MQ=37)
                                                                                     cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat             >  1:1048252/1‑100 (MQ=38)
                                                                                     cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat             >  1:1273934/1‑100 (MQ=255)
                                                                                      ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg            >  1:327971/1‑100 (MQ=255)
                                                                                       cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga           >  2:562040/1‑100 (MQ=38)
                                                                                       cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga           <  1:3198629/100‑1 (MQ=255)
                                                                                         cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg         <  2:749749/100‑1 (MQ=255)
                                                                                          ggAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGt        >  1:1843769/1‑100 (MQ=255)
                                                                                               gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc   >  1:2051005/1‑100 (MQ=255)
                                                                                               gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc   <  2:708549/100‑1 (MQ=255)
                                                                                               gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc   >  1:4103715/1‑100 (MQ=255)
                                                                                                ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCg     <  2:1594883/97‑1 (MQ=255)
                                                                                                ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa  >  2:2841089/1‑100 (MQ=255)
                                                                                                |                                                                                                   
AATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCA  >  NC_000913/2560225‑2560420

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