Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,047 C→T 80.5% intergenic (‑306/+352) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,0470CT80.5% 405.4 / 137.3 293intergenic (‑306/+352)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base C (15/42);  new base T (102/132);  total (117/176)
Fisher's exact test for biased strand distribution p-value = 2.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 294 total aligned reads displayed.

GGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA  >  NC_000913/4325989‑4326145
                                                          |                                                                                                  
ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                           <  2:2010293/100‑1 (MQ=11)
 gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          <  2:1453528/100‑1 (MQ=14)
 gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          >  2:707473/1‑100 (MQ=14)
   gCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGt                                                        <  2:2526359/100‑1 (MQ=17)
    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATgggg                                                       >  1:1485421/1‑98 (MQ=9)
    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATggag                                                       >  1:841407/1‑98 (MQ=9)
    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGCGAACGCCTTATCCGGCCTACCGAtcac                                                        >  1:52216/1‑96 (MQ=11)
     gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGc                                                      >  2:673275/1‑100 (MQ=17)
     gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          >  2:2270819/1‑96 (MQ=31)
      ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg                                                     >  2:2468281/1‑100 (MQ=14)
      cagtgccaactgcctgatgcgacgcttgcgcgtcttATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                     <  1:1651264‑M1/64‑1 (MQ=255)
       aagatcgTAGGCCTGATAAGACGCGCAAGCGTCGCATCAAGAAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTAATCCGGCCTACAGATGGTGCgt                                                    <  1:1171779/94‑2 (MQ=11)
        aaTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCggg                                                   <  1:2715451/100‑1 (MQ=33)
          gccgactgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGCATGGtgccg                                                    <  1:176152‑M1/65‑3 (MQ=255)
          ccgatctgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGaa                                                 >  2:581641‑M1/33‑100 (MQ=255)
           ccgactgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGCGCGGGAAt                                                <  1:2024096‑M1/69‑1 (MQ=255)
           ccgactgcctgatgcgacgctcgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                <  2:2207361‑M1/69‑1 (MQ=255)
            tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAtt                                               <  1:2294312/100‑1 (MQ=18)
            tGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCTGGCCTACGGATGGTGCGGGAAt                                                <  1:1685319/99‑1 (MQ=11)
            cgactgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAAt                                                <  1:511840‑M1/69‑1 (MQ=255)
             gactgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                              <  2:1873924‑M1/71‑1 (MQ=255)
             gactgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                              <  1:733893‑M1/71‑1 (MQ=255)
             gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                              <  2:2517288/100‑1 (MQ=18)
              tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                             <  2:1543539/100‑1 (MQ=17)
              tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                             <  2:1101752/100‑1 (MQ=18)
               ctgcctgatgcgacgcttgcgcgtcttATCAGGCCGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGAGGGTGCGGGAACTTGt                                            <  1:2530970‑M1/73‑1 (MQ=255)
               ctgcctgatgcgacgcttgcgcgtcttATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                            <  2:2666479‑M1/73‑1 (MQ=255)
               ctgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgggggcggggaattg                                             <  1:2414460‑M1/72‑7 (MQ=255)
               ctgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                            <  2:1282902‑M1/73‑1 (MQ=255)
               aGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                            >  2:1675283/1‑100 (MQ=17)
               aGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                            >  2:847899/1‑100 (MQ=18)
                tgcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                           <  2:234079‑M1/74‑1 (MQ=255)
                ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGGATTTGt                                            >  1:1565552/1‑99 (MQ=12)
                ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                           <  1:1983252/100‑1 (MQ=18)
                ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                           >  1:2825968/1‑100 (MQ=14)
                ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                           >  2:2769449/1‑100 (MQ=18)
                 gcctgatgcgacgcttgcgcgtcttATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                          <  2:319167‑M1/75‑1 (MQ=255)
                 gcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGGCTACGGATGGCGCGGGAatgtgttg                                          <  2:1901084‑M1/75‑3 (MQ=255)
                 gcctgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtttggag                                          <  1:566656‑M1/75‑4 (MQ=255)
                 gactgatgcgacgattgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                          <  2:1976742‑M1/75‑1 (MQ=255)
                  ccagATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  1:1429962/97‑1 (MQ=16)
                  ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTTTAgg                                         >  2:1789061/1‑100 (MQ=9)
                  ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  2:1233449/100‑1 (MQ=17)
                   cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGc                                        <  2:1013323/100‑1 (MQ=18)
                    tgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                       <  1:372300‑M1/78‑1 (MQ=255)
                    tgatgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGcc                                       <  2:1407616‑M1/78‑1 (MQ=255)
                    tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                       <  2:1144610/100‑1 (MQ=18)
                      atgcgacgcttgcgcgtcgcATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTg                                     <  2:158169‑M1/80‑1 (MQ=255)
                      aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTg                                     >  1:2914503/1‑100 (MQ=14)
                      aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTg                                     <  2:1619525/100‑1 (MQ=18)
                       tgcgacgcttgcgcgtcttATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                    <  2:287618‑M1/81‑1 (MQ=255)
                       tAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                    >  1:2512028/1‑100 (MQ=32)
                        aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                   <  2:837546/100‑1 (MQ=18)
                        aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                   >  2:798995/1‑100 (MQ=17)
                         aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGGATTTGTAGGCCTGATa                                  >  1:779653/1‑100 (MQ=12)
                         aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                  <  2:90540/100‑1 (MQ=17)
                         aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                  <  2:77863/100‑1 (MQ=18)
                         aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                  <  2:2598600/100‑1 (MQ=33)
                         aGAAGCGCAAGCGGCGCGTCAGGCGGTCGGCACTGTTGCCGGATGTGGCGTAAACGCTTTCTCCGGCCTACGGATGGGGCGGGAATTTGTAGGcc                                       <  1:908577/95‑1 (MQ=11)
                          gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 >  2:448836/1‑100 (MQ=17)
                          gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 >  1:1276860/1‑100 (MQ=17)
                          gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 <  1:2751497/100‑1 (MQ=18)
                          gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 >  2:613708/1‑100 (MQ=18)
                          gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 <  1:2669497/100‑1 (MQ=18)
                          gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                  <  1:1373984/99‑1 (MQ=18)
                            cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGTGCGGGAATTTGTCGGCCTGTTAAg                                >  1:1963756/1‑99 (MQ=11)
                            cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGTTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:212868/1‑100 (MQ=12)
                            cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:1180825/1‑100 (MQ=17)
                            cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:1856386/1‑100 (MQ=17)
                            cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTTGTCCTGATAAGa                               >  1:611068/1‑100 (MQ=11)
                             gcgcAAGCGTCGCATCAGTCAGTCCGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAc                              <  1:2061536/100‑1 (MQ=255)
                             gcgcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAc                              <  2:2726764/100‑1 (MQ=17)
                             gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:1321350/1‑99 (MQ=33)
                             gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAc                              >  2:2505144/1‑100 (MQ=33)
                              cttgcgcgtcttATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGAc                              <  1:2863425‑M1/87‑1 (MQ=255)
                              cgcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                             <  1:2084485/100‑1 (MQ=18)
                               gcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                            >  1:2215031/1‑100 (MQ=17)
                               gcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                            >  2:1153261/1‑100 (MQ=14)
                               gcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                             <  1:2332746/99‑1 (MQ=18)
                                cAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                           >  2:2218667/1‑100 (MQ=11)
                                 aaGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                           >  2:1377116/1‑99 (MQ=1)
                                 aaGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                          <  2:2350226/100‑1 (MQ=1)
                                        gCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                   <  1:1046943/100‑1 (MQ=32)
                                          ATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                  >  1:2138537‑M1/1‑99 (MQ=255)
                                          ATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                 >  2:2420542‑M1/1‑100 (MQ=255)
                                          ATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                 >  2:556423‑M1/1‑100 (MQ=255)
                                             AGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc              <  1:111487‑M1/100‑1 (MQ=255)
                                             AGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc              <  1:2315890‑M1/100‑1 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCGTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa             <  2:2053631‑M1/100‑1 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGGTAGAACGCGCAGGCGGCGCATc              <  1:1189037‑M1/99‑1 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa             >  2:2835434‑M1/1‑100 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa             >  2:2397123‑M1/1‑100 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa             <  2:1843475‑M1/100‑1 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa             <  2:1513981‑M1/100‑1 (MQ=255)
                                              GGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAACGTCGCATCa             <  1:1915355‑M1/100‑1 (MQ=255)
                                               GCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg            >  2:1275227‑M1/1‑100 (MQ=255)
                                                      gCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTg     >  2:958226/1‑100 (MQ=17)
                                                      gCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATTTg     >  2:482014/1‑100 (MQ=18)
                                                       cACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTgg    >  2:2033706/1‑100 (MQ=17)
                                                       cACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTgg    <  1:593668/100‑1 (MQ=17)
                                                          |                                                                                                  
GGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA  >  NC_000913/4325989‑4326145

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