Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,621,612 T→C intergenic (+199/‑397) yhhH → / → yhhI hypothetical protein/putative receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,621,6120TC79.2% 536.5 / 170.4 294intergenic (+199/‑397)yhhH/yhhIhypothetical protein/putative receptor
Reads supporting (aligned to +/- strand):  ref base T (36/25);  new base C (151/82);  total (187/107)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 294 total aligned reads displayed.

AAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATT  >  NC_000913/3621523‑3621702
                                                                                         |                                                                                          
aaaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtga                                                                                         <  1:4545745/93‑1 (MQ=34)
aaaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtga                                                                                         <  2:7270970/93‑1 (MQ=34)
aaaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtga                                                                                         >  2:4625129/1‑93 (MQ=21)
aaaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtga                                                                                         <  1:1886088/93‑1 (MQ=21)
 aaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtgat                                                                                        <  2:6334137/93‑1 (MQ=21)
 aaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtgat                                                                                        >  1:3993941/1‑93 (MQ=21)
 aaTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtgat                                                                                        >  2:5689827/1‑93 (MQ=21)
  atacatacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtgatg                                                                                       >  2:5881601/1‑93 (MQ=21)
   tacatacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCtgatga                                                                                      >  2:5207086/1‑93 (MQ=18)
    acatacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTgatgat                                                                                     <  2:7524761/93‑1 (MQ=34)
     catacCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATc                                                                                    >  1:5871029/1‑93 (MQ=21)
         ccAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCaa                                                                                <  2:5461090/93‑1 (MQ=255)
          cAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAAt                                                                               >  2:5779558/1‑93 (MQ=255)
          cAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAAt                                                                               >  1:5159034/1‑93 (MQ=255)
          cAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAAt                                                                               <  1:2939012/93‑1 (MQ=255)
           agaaaagaaaGGCGAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAAtt                                                                              <  2:2069302/93‑1 (MQ=255)
            gaaaagaaaGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTc                                                                             <  1:1323483/93‑1 (MQ=39)
            gaaaagaaaGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTc                                                                             >  2:4577541/1‑93 (MQ=255)
             aaaagaaaGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTcc                                                                            >  2:3343964/1‑93 (MQ=255)
              aaagaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCt                                                                           <  2:4778687/93‑1 (MQ=38)
              aaagaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCt                                                                           >  2:3431930/1‑93 (MQ=37)
               aagaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCtt                                                                          <  1:1450427/93‑1 (MQ=255)
               aagaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCtt                                                                          <  2:1634197/93‑1 (MQ=255)
                agaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTg                                                                         >  2:4389684/1‑93 (MQ=255)
                 gaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGa                                                                        >  1:2197985/1‑93 (MQ=38)
                 gaaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGa                                                                        >  2:489007/1‑93 (MQ=255)
                  aaagGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGaa                                                                       >  2:572797/1‑93 (MQ=255)
                   aagGCAAGTAAAGAGTATGGTTTTTATAATAAATGGAAAAAACTAAATGATTATGAATTATTTCGCTTTCCTTATTGTCGCAATTTCTTGaaa                                                                      >  1:3640549/1‑93 (MQ=11)
                    agGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGaaaa                                                                     <  2:2633405/93‑1 (MQ=37)
                      gCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAgg                                                                   >  2:2260194/1‑93 (MQ=255)
                       cAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGa                                                                  <  1:5391776/93‑1 (MQ=255)
                         aGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGAtt                                                                <  2:1237635/93‑1 (MQ=37)
                           tAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTc                                                              >  1:4428083/1‑93 (MQ=255)
                           tAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTc                                                              >  2:2913937/1‑93 (MQ=255)
                             aaGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCAt                                                            >  1:2516751/1‑93 (MQ=255)
                              agagTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATc                                                           <  1:1168855/93‑1 (MQ=255)
                              agagTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATc                                                           <  2:357275/93‑1 (MQ=255)
                              agagTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATc                                                           >  2:5504757/1‑93 (MQ=255)
                               gagTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCt                                                          >  2:1877915/1‑93 (MQ=255)
                               gagTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCt                                                          >  1:1841255/1‑93 (MQ=255)
                                 gTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGc                                                        <  1:6064481/93‑1 (MQ=255)
                                  tATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGcc                                                       >  2:6991425/1‑93 (MQ=255)
                                    tGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCag                                                     >  1:2640668/1‑93 (MQ=39)
                                    tGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCag                                                     >  2:2823303/1‑93 (MQ=39)
                                    tGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCag                                                     >  1:1261993/1‑93 (MQ=255)
                                    tGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCag                                                     >  2:2138739/1‑93 (MQ=255)
                                     ggTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCaga                                                    <  1:3859807/93‑1 (MQ=255)
                                         tatataATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTAt                                                <  1:657726/93‑1 (MQ=38)
                                         tatataATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTAt                                                >  1:1102856/1‑93 (MQ=255)
                                           tataATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTa                                              >  1:5838240/1‑93 (MQ=255)
                                           tataATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTa                                              >  1:819866/1‑93 (MQ=255)
                                            ataataAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTAc                                             <  1:375268/93‑1 (MQ=255)
                                            ataataAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTAc                                             >  1:7018341/1‑93 (MQ=255)
                                            ataataAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTAc                                             <  1:946055/93‑1 (MQ=37)
                                             taataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACa                                            <  2:3110828/93‑1 (MQ=255)
                                               ataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGt                                          >  1:4696044/1‑93 (MQ=255)
                                               ataaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGt                                          >  1:7157354/1‑93 (MQ=255)
                                                taaATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGtt                                         >  1:5659165/1‑93 (MQ=255)
                                                  aaTGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTaa                                       >  1:6348248/1‑93 (MQ=21)
                                                   aTGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAag                                      >  2:6054718/1‑93 (MQ=34)
                                                   aTGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAag                                      >  2:2417350/1‑93 (MQ=18)
                                                    tGTAAGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAaga                                     >  2:6481334/1‑93 (MQ=21)
                                                     gTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAagag                                    >  2:4224184/1‑93 (MQ=33)
                                                       aaGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAggg                                  <  1:6774091/93‑1 (MQ=25)
                                                        aGAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAgggg                                 <  1:6756569/93‑1 (MQ=255)
                                                         gAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAggggg                                >  1:7370102/1‑93 (MQ=255)
                                                         gAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAggggg                                >  1:7469961/1‑93 (MQ=255)
                                                         gAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAggggg                                >  2:3795635/1‑93 (MQ=255)
                                                         gAAACTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAggggg                                <  1:5502830/93‑1 (MQ=255)
                                                           aaCTAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGc                              <  1:2873702/93‑1 (MQ=255)
                                                              tAAATGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAg                           >  1:1205781/1‑93 (MQ=37)
                                                               aaaTGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGa                          >  1:3243975/1‑93 (MQ=37)
                                                                aaTGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGAc                         >  1:4957382/1‑93 (MQ=38)
                                                                aaTGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGAc                         <  2:7341617/93‑1 (MQ=255)
                                                                aaTGATGATGAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGAc                         >  2:5373286/1‑93 (MQ=255)
                                                                 atgatgatgAATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACg                        >  1:2692482/1‑93 (MQ=255)
                                                                  tgatgatgaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGc                       <  1:1011032/93‑1 (MQ=38)
                                                                     tgatgaATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGt                    >  1:1709067/1‑93 (MQ=38)
                                                                     tgatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGt                    >  2:387907/1‑93 (MQ=255)
                                                                     tgatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGt                    >  2:2519623/1‑93 (MQ=255)
                                                                      gatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGtt                   >  2:1287023/1‑93 (MQ=255)
                                                                      gatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGtt                   >  1:5992886/1‑93 (MQ=255)
                                                                      gatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGtt                   >  2:611377/1‑93 (MQ=255)
                                                                      gatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGtt                   >  2:5638392/1‑93 (MQ=255)
                                                                      gatgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGtt                   <  1:212995/93‑1 (MQ=255)
                                                                       atgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTa                  >  2:7179585/1‑93 (MQ=255)
                                                                        tgaATTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAg                 <  2:1215038/93‑1 (MQ=255)
                                                                          aaTTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGAt               <  1:436828/93‑1 (MQ=38)
                                                                          aaTTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGAt               >  2:6188923/1‑93 (MQ=38)
                                                                          aaTTATTTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGAt               >  1:4191637/1‑93 (MQ=255)
                                                                           attattTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGAtt              <  2:4309826/93‑1 (MQ=38)
                                                                           attattTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGAtt              <  1:6149047/93‑1 (MQ=37)
                                                                              attTCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGc           >  2:6716131/1‑93 (MQ=255)
                                                                               tttCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCa          >  2:6609502/1‑93 (MQ=38)
                                                                               tttCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCa          >  2:1700347/1‑93 (MQ=255)
                                                                               tttCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCa          >  2:1480126/1‑93 (MQ=255)
                                                                                 tCGCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCAAt        <  2:4165071/93‑1 (MQ=255)
                                                                                   gCTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCAATTg      >  2:315039/1‑93 (MQ=255)
                                                                                    cTTACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCAATTGa     >  1:200065/1‑93 (MQ=255)
                                                                                      tACCTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCAATTGAAt   >  1:6944642/1‑93 (MQ=255)
                                                                                         |                                                                                          
AAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCGCAATTCCTTGAAAAGGATTTCATCTGCCAGAGTATTACAGTTAAGAGGGGGGCAAGACGCTGTTAGATTGGCAATTGAATT  >  NC_000913/3621523‑3621702

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