Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,330 G→T 15.1% C25F (TGC→TTC)  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,3300GT15.1% 511.1 / 69.7 186C25F (TGC→TTC) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (74/84);  new base T (16/12);  total (90/96)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 189 total aligned reads displayed.

TATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAA  >  NC_000913/2560232‑2560429
                                                                                                  |                                                                                                   
tattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                                    <  1:1148066/100‑1 (MQ=255)
  ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                  <  1:30767/100‑1 (MQ=255)
    cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                                <  1:3137197/100‑1 (MQ=38)
    cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                                >  2:2075335/1‑100 (MQ=255)
    cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                                >  1:1690507/1‑100 (MQ=255)
    cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                                <  1:2438152/100‑1 (MQ=255)
     aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                               >  2:1154605/1‑100 (MQ=255)
      atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                              <  2:1927339/100‑1 (MQ=38)
      atatCAAGAACTTAAAGCCATGTTTGCAATTATAGATGCTGCGGCGGACACCAAGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                              <  2:264867/100‑1 (MQ=9)
       tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAATATCGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                               <  2:1443692/98‑1 (MQ=37)
       tatCAAGAACTTAAATCCATGCTTGCAACGAAAGATCCATAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                             <  2:1378956/100‑1 (MQ=255)
        atCAAAAACTTAAGTCCGCGTTTGCAACGAAAGATCCAGGGTTGGAAAATAGGATCAACACAAATAAAAGCCCCCTGCATGTTGCAGCATTCCggggg                                                                                              <  2:2393602/98‑1 (MQ=25)
         tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           <  2:302222/100‑1 (MQ=255)
         tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           >  1:1889903/1‑100 (MQ=255)
          caacaaCTTCAATCCCTCTTTGCATCTTATGATCCTACGTCCACAAACAGGATCACCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           <  2:4080503/99‑1 (MQ=255)
          cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGACATGTTGCAGCATTCCGGGGGcttt                                                                                          <  1:2107368/100‑4 (MQ=37)
              aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                                      <  1:1516337/100‑1 (MQ=255)
               aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                     <  2:3458884/100‑1 (MQ=255)
                tttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                    <  1:3839942/99‑1 (MQ=255)
                ctcagaTCCATGTTTGCAACGAAACATCCATAGCTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTCTga                                                                                      <  2:3078636/93‑1 (MQ=255)
                cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                    <  1:1614272/100‑1 (MQ=255)
                cTTAAATCCATGTTTGCAACGAAAGATCAAGAGTTAGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAAATGCCGGGGGCGGTATgaga                                                                                    >  2:553095/1‑100 (MQ=255)
                 ttcaaTCCATGCTTCCAACAAAAGATCCACAGTTAGCAAATAGGATCACCACAACAACACGACCCAGGAAAGCTGCAACATGCCGGGGGAGGTATgagac                                                                                   <  2:3293000/97‑2 (MQ=255)
                 ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                   <  2:366410/100‑1 (MQ=255)
                 ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                   >  2:3839917/1‑100 (MQ=255)
                  tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                                  >  2:722641/1‑100 (MQ=255)
                  naaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                                  >  2:721706/2‑100 (MQ=255)
                   aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                                 <  1:1742438/100‑1 (MQ=255)
                        cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACGAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                             >  1:2380008/1‑99 (MQ=255)
                          tGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaa                                                                          <  1:3964255/100‑1 (MQ=255)
                           gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                         >  1:15971/1‑100 (MQ=255)
                           gTTTGCAACGAAAGATCCAGAGCTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCAGGGGCGGTATGAGAAACAAGCTaa                                                                          >  2:849528/1‑99 (MQ=255)
                             ttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                        <  2:90818/99‑1 (MQ=255)
                                cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                                    >  2:3378621/1‑100 (MQ=255)
                                 aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                                   <  1:2419027/100‑1 (MQ=255)
                                 aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                                   >  1:3650712/1‑100 (MQ=255)
                                    gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                                >  1:2457343/1‑100 (MQ=255)
                                      aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                              >  2:499818/1‑100 (MQ=255)
                                       cttccgatctttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                             <  1:1129350/90‑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)
                                            cAGAGTTCGAAAATAGGATCAACACAGATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                        >  1:2512211/1‑100 (MQ=255)
                                            cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                        <  1:561455/100‑1 (MQ=255)
                                            cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                        <  2:752120/100‑1 (MQ=255)
                                            cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                        >  1:1380201/1‑100 (MQ=255)
                                             agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                                       >  1:2416671/1‑100 (MQ=255)
                                                gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCaagatcgga                                                    >  2:2565921/1‑92 (MQ=255)
                                                gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                    <  2:1577294/100‑1 (MQ=255)
                                                gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                    <  2:2151890/100‑1 (MQ=255)
                                                 ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                                   <  1:2487140/100‑1 (MQ=255)
                                                      aaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                              >  1:3048785/1‑100 (MQ=255)
                                                         tAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                                           <  1:836880/100‑1 (MQ=255)
                                                            gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                                        >  2:3508798/1‑100 (MQ=255)
                                                              tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGa                                      >  2:3195418/1‑100 (MQ=255)
                                                                aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                    >  2:313615/1‑100 (MQ=255)
                                                                aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                    >  2:3790298/1‑100 (MQ=255)
                                                                aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                    >  2:2932807/1‑100 (MQ=255)
                                                                aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                    <  1:1952052/100‑1 (MQ=255)
                                                                 acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                                   >  2:1593206/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)
                                                                     aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaa                               <  2:3015600/100‑1 (MQ=255)
                                                                      aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCGAAAGGAAAATTTCTCATGTCTGATATGCTAGCAGTCGAACAAGaaa                              <  2:562238/100‑1 (MQ=255)
                                                                      aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAAAAGCTAAAGGAGAATTTCTCATGTCTGATATGATAGCAGTAGAACAAGaaa                              >  1:2108570/1‑100 (MQ=255)
                                                                       cataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                             >  1:1587780/7‑100 (MQ=25)
                                                                       acaaaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                             <  2:2505383/98‑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)
                                                                           ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                         >  1:2495375/3‑100 (MQ=37)
                                                                           caacccccGGAATACTGCAACATGCCGCTGGCGGTATGAGAAGCAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGTCGTAGAACAAGaatcaacc                         <  2:3571888/97‑3 (MQ=255)
                                                                           aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                         <  1:1393148/100‑1 (MQ=255)
                                                                            aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                        >  2:623562/1‑100 (MQ=255)
                                                                            aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                        <  2:804191/100‑1 (MQ=255)
                                                                             gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                       >  1:458121/1‑100 (MQ=38)
                                                                             gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTACCAGTAGAACAAGAaataaaca                        >  2:3982801/1‑99 (MQ=37)
                                                                             gCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                       >  1:3883326/1‑100 (MQ=255)
                                                                              cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAaaaaacaa                       >  1:1464475/1‑99 (MQ=37)
                                                                              cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      >  1:1048252/1‑100 (MQ=38)
                                                                              cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      >  2:1703203/1‑100 (MQ=38)
                                                                              cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      >  1:742203/1‑100 (MQ=38)
                                                                              cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      >  1:1273934/1‑100 (MQ=255)
                                                                               ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg                     >  1:327971/1‑100 (MQ=255)
                                                                                cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                    >  1:2319239/1‑100 (MQ=38)
                                                                                cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                    <  1:3198629/100‑1 (MQ=255)
                                                                                   ggAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGt                 >  1:1843769/1‑100 (MQ=255)
                                                                                         ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa           >  2:2841089/1‑100 (MQ=255)
                                                                                          tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCaa          >  1:2355530/1‑100 (MQ=255)
                                                                                           gcAACATGCCGGGGGCGGTATGAGAAACACGCTAAAGGAGAATTTATCATGTCTGATATGCTAGCAGCAGAACACGAAACAAACAATGCTGTTCGTCAAt         >  2:1633/1‑100 (MQ=255)
                                                                                               cATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttttt     <  2:3761442/100‑1 (MQ=255)
                                                                                               cATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttttt     <  1:2120035/100‑1 (MQ=255)
                                                                                                aTGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa    >  1:2081226/1‑100 (MQ=255)
                                                                                                 tGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTaa   >  2:1111749/1‑100 (MQ=255)
                                                                                                  gCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTaaa  <  2:1872537/100‑1 (MQ=255)
                                                                                                  gCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTaaa  <  1:404940/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAA  >  NC_000913/2560232‑2560429

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: