Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,916,745 T→C V163A (GTA→GCA)  yebS → intermembrane transport protein YebS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,916,7450TC100.0% 1286.4 / NA 310V163A (GTA→GCA) yebSintermembrane transport protein YebS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (145/165);  total (145/165)

Only 100 of 310 total aligned reads displayed.

CATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTTCTCTTTTGTG  >  NC_000913/1916647‑1916842
                                                                                                  |                                                                                                 
 aTTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAAt                                                                                                >  1:3210426/1‑101 (MQ=255)
   tCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGc                                                                                              <  2:3209936/101‑1 (MQ=255)
    cTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCt                                                                                             >  2:3210432/1‑101 (MQ=255)
    cTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCt                                                                                             >  1:3210430/1‑101 (MQ=255)
        tGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGat                                                                                         <  1:3210158/101‑1 (MQ=255)
            cTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCt                                                                                     >  2:3210436/1‑101 (MQ=255)
             tCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTa                                                                                    >  1:3210439/1‑101 (MQ=255)
             tCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTa                                                                                    <  1:3210008/101‑1 (MQ=255)
              ccATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTAc                                                                                   <  2:3210119/101‑1 (MQ=255)
              ccATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTAc                                                                                   >  2:3210440/1‑101 (MQ=255)
                aTAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCt                                                                                 >  2:3210444/1‑101 (MQ=255)
                 tAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTg                                                                                >  1:3210448/1‑101 (MQ=255)
                 tAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTg                                                                                >  2:3210449/1‑101 (MQ=255)
                     ttatttatGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCg                                                                            >  2:3210454/1‑101 (MQ=255)
                      tatttatGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCgg                                                                           >  1:3210456/1‑101 (MQ=255)
                       atttatGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGc                                                                          <  2:3210039/101‑1 (MQ=255)
                          tatGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCAtt                                                                       <  2:3209985/101‑1 (MQ=255)
                          tatGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCAtt                                                                       >  2:3210460/1‑101 (MQ=255)
                          tatGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCAtt                                                                       <  2:3209871/101‑1 (MQ=255)
                           atGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTg                                                                      <  1:3209574/101‑1 (MQ=255)
                            tGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTgg                                                                     <  1:3210018/101‑1 (MQ=255)
                             ggTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGc                                                                    <  1:3209845/101‑1 (MQ=255)
                             ggTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGc                                                                    <  1:3209809/101‑1 (MQ=255)
                             ggTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGc                                                                    >  1:3210464/1‑101 (MQ=255)
                             ggTGTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGc                                                                    <  1:3209937/101‑1 (MQ=255)
                               tttGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGt                                                                  >  2:3210465/1‑101 (MQ=255)
                                ttGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGtt                                                                 <  1:3210101/101‑1 (MQ=255)
                                ttGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGtt                                                                 <  2:3209618/101‑1 (MQ=255)
                                    tAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCtt                                                             >  1:3210471/1‑101 (MQ=255)
                                    tAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCtt                                                             >  2:3210469/1‑101 (MQ=255)
                                      aCCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCt                                                           >  1:3210476/1‑101 (MQ=255)
                                      aCCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCt                                                           <  2:3209906/101‑1 (MQ=255)
                                       ccGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCta                                                          <  2:3210043/101‑1 (MQ=255)
                                        cGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCtat                                                         >  2:3210478/1‑101 (MQ=255)
                                        cGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCtat                                                         <  1:3210284/101‑1 (MQ=255)
                                        cGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCtat                                                         <  2:3210277/101‑1 (MQ=255)
                                        cGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCtat                                                         <  1:3209918/101‑1 (MQ=255)
                                         gACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCtata                                                        >  2:3210481/1‑101 (MQ=255)
                                          aCTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATaa                                                       >  1:3210485/1‑101 (MQ=255)
                                          aCTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATaa                                                       <  2:3209843/101‑1 (MQ=255)
                                           cTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATaaa                                                      <  2:3210010/101‑1 (MQ=255)
                                           cTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATaaa                                                      <  1:3210085/101‑1 (MQ=255)
                                            tGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAg                                                     >  2:3210489/1‑101 (MQ=255)
                                             gggAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAgg                                                    <  2:3210064/101‑1 (MQ=255)
                                             gggAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAgg                                                    >  1:3210492/1‑101 (MQ=255)
                                             gggAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAgg                                                    >  2:3210493/1‑101 (MQ=255)
                                              ggAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCa                                                                         <  1:3209951/79‑1 (MQ=255)
                                              ggAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGt                                                   <  1:3209801/101‑1 (MQ=255)
                                                aatgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTAc                                                 <  1:3209879/101‑1 (MQ=255)
                                                aatgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTAc                                                 >  1:3210498/1‑101 (MQ=255)
                                                 atgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACa                                                <  2:3209884/101‑1 (MQ=255)
                                                  tgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAg                                               <  1:3209925/101‑1 (MQ=255)
                                                  tgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAg                                               <  2:3210206/101‑1 (MQ=255)
                                                  tgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAg                                               <  2:3210186/101‑1 (MQ=255)
                                                  tgaatTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAg                                               >  1:3210501/1‑101 (MQ=255)
                                                     atTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGAt                                            >  2:3210503/1‑101 (MQ=255)
                                                       ttACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTa                                          <  1:3210051/101‑1 (MQ=255)
                                                       ttACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTa                                          <  2:3210102/101‑1 (MQ=255)
                                                        tACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTAt                                         <  1:3210222/101‑1 (MQ=255)
                                                         aCGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATg                                        <  2:3210086/101‑1 (MQ=255)
                                                           gTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGcc                                      >  2:3210507/1‑101 (MQ=255)
                                                             ccGGTACTGCTGATGCTTGAGCGACTCAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCa                                    >  1:3210511/1‑101 (MQ=255)
                                                              cGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCat                                   <  1:3210190/101‑1 (MQ=255)
                                                               ggTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCata                                  >  2:3210513/1‑101 (MQ=255)
                                                                gTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCatat                                 <  2:3209971/101‑1 (MQ=255)
                                                                gTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCatat                                 >  2:3210518/1‑101 (MQ=255)
                                                                gTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCatat                                 <  2:3210352/101‑1 (MQ=255)
                                                                 tACTGCTGATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATc                                >  1:3210519/1‑101 (MQ=255)
                                                                   ctgctgATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCa                              >  2:3210524/1‑101 (MQ=255)
                                                                    tgctgATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAg                             <  1:3210151/101‑1 (MQ=255)
                                                                    tgctgATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAg                             >  2:3210527/1‑101 (MQ=255)
                                                                    tgctgATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAg                             <  1:3210164/101‑1 (MQ=255)
                                                                     gctgATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAgg                            >  1:3210529/1‑101 (MQ=255)
                                                                      ctgATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGc                           <  2:3210162/101‑1 (MQ=255)
                                                                        gATGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTAACTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTg                         >  1:3210533/1‑101 (MQ=255)
                                                                         aTGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTgg                        <  2:3210133/101‑1 (MQ=255)
                                                                         aTGCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTgg                        <  2:3210078/101‑1 (MQ=255)
                                                                           gCTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGgtg                      <  1:3210070/101‑1 (MQ=255)
                                                                            cTTGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGgtgt                     >  2:3210535/1‑101 (MQ=255)
                                                                             ttGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCGGGCTGGTGTc                    <  1:3209992/101‑1 (MQ=255)
                                                                             ttGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTc                    >  1:3210539/1‑101 (MQ=255)
                                                                              tGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCg                   <  1:3210326/101‑1 (MQ=255)
                                                                              tGAGCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCg                   >  2:3210540/1‑101 (MQ=255)
                                                                                 gCGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCt                <  1:3210136/101‑1 (MQ=255)
                                                                                  cGACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCtt               >  2:3210544/1‑101 (MQ=255)
                                                                                   gACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTg              <  2:3209978/101‑1 (MQ=255)
                                                                                   gACTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTg              >  1:3210546/1‑101 (MQ=255)
                                                                                    aCTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGt             <  2:3210217/101‑1 (MQ=255)
                                                                                    aCTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGt             <  2:3209949/101‑1 (MQ=255)
                                                                                     cTTAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGtt            >  2:3210549/1‑101 (MQ=255)
                                                                                      ttAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTtc           >  1:3210552/1‑101 (MQ=255)
                                                                                      ttAAAGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTtc           <  1:3209966/101‑1 (MQ=255)
                                                                                        aaaGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTtctc         >  2:3210558/1‑101 (MQ=255)
                                                                                        aaaGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTtctc         >  1:3210557/1‑101 (MQ=255)
                                                                                        aaaGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTtctc         <  1:3210037/101‑1 (MQ=255)
                                                                                         aaGAGTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTTctct        <  1:3210114/101‑1 (MQ=255)
                                                                                          agagTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTtct          >  2:3210561/1‑98 (MQ=255)
                                                                                          agagTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTTCTCtt       <  2:3209812/101‑1 (MQ=255)
                                                                                            agTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTTCTCtttt     >  1:3210565/1‑101 (MQ=255)
                                                                                            agTGGGCAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTTCTCtttt     <  1:3210415/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
CATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTCGATATCTACCTGGTCGGCATTGGCGTTGCTTCTATAAAGGTACAGGATTATGCCCATATCCAGGCTGGTGTCGGCTTGTTCTCTTTTGTG  >  NC_000913/1916647‑1916842

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