Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,842,228 G→A G213D (GGT→GAT)  yicP → putative adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,842,2280GA85.1% 973.6 / 115.5 355G213D (GGT→GAT) yicPputative adenine deaminase
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (139/163);  total (163/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 355 total aligned reads displayed.

CTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAG  >  NC_000913/3842137‑3842317
                                                                                           |                                                                                         
cTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGt                                                                                          >  1:3023146/1‑93 (MQ=255)
 ttGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATc                                                                                         >  2:1682209/1‑93 (MQ=255)
 ttGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATc                                                                                         >  2:6949089/1‑93 (MQ=255)
   gCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCAc                                                                                       <  2:2486124/93‑1 (MQ=255)
     aGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTg                                                                                     >  1:7052423/1‑93 (MQ=255)
     aGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTg                                                                                     >  2:5734950/1‑93 (MQ=255)
     aGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTg                                                                                     <  1:4286414/93‑1 (MQ=255)
      gAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGc                                                                                    >  1:794231/1‑93 (MQ=255)
        aaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGccc                                                                                  <  2:4986447/93‑1 (MQ=255)
          tgatgGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCgg                                                                                >  2:5796639/1‑93 (MQ=255)
             tgGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGtt                                                                             <  1:4968266/93‑1 (MQ=255)
             tgGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGtt                                                                             <  2:6213527/93‑1 (MQ=255)
              ggACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGttt                                                                            <  2:3169876/93‑1 (MQ=255)
               gACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTg                                                                           <  1:3368015/93‑1 (MQ=255)
                aCTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTgg                                                                          <  2:489925/93‑1 (MQ=255)
                 cTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTggg                                                                         >  2:1195841/1‑93 (MQ=255)
                 cTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTggg                                                                         >  1:7313408/1‑93 (MQ=255)
                 cTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTggg                                                                         <  2:7059983/93‑1 (MQ=255)
                  tACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGggt                                                                        >  1:4531956/1‑93 (MQ=255)
                   aCCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGggtg                                                                       <  2:7371148/93‑1 (MQ=255)
                     ccTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGggtggt                                                                     >  2:4072355/1‑93 (MQ=255)
                     ccTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGggtggt                                                                     <  1:3561774/93‑1 (MQ=255)
                       tGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCGCGGGTTTGGGTGGTaa                                                                   >  2:3506200/1‑93 (MQ=255)
                       tGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTaa                                                                   <  2:3957564/93‑1 (MQ=255)
                        ggCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTaaa                                                                  <  1:4054394/93‑1 (MQ=255)
                        ggCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTaaa                                                                  <  2:3639947/93‑1 (MQ=255)
                         gCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAg                                                                 >  1:1701241/1‑93 (MQ=255)
                         gCGNAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAg                                                                 <  1:4386486/93‑1 (MQ=255)
                           gTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGaa                                                               <  1:1372444/93‑1 (MQ=255)
                            tAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAAc                                                              <  1:102477/93‑1 (MQ=255)
                            tAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAAc                                                              >  2:4458896/1‑93 (MQ=255)
                              aTTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACtt                                                            >  1:2327917/1‑93 (MQ=255)
                              aTTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACtt                                                            <  1:2110673/93‑1 (MQ=255)
                               ttAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTa                                                           >  2:5290687/1‑93 (MQ=255)
                                  gCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACg                                                        <  1:3204527/93‑1 (MQ=255)
                                       cAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCtat                                                   <  2:2024226/93‑1 (MQ=255)
                                        aGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCtata                                                  >  1:3381340/1‑93 (MQ=255)
                                        aGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCtata                                                  <  2:1296011/93‑1 (MQ=255)
                                         gAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTatat                                                 >  2:4147189/1‑93 (MQ=255)
                                         gAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTatat                                                 <  2:6437967/93‑1 (MQ=255)
                                         gAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTatat                                                 >  2:6514392/1‑93 (MQ=255)
                                           aTGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTa                                               <  2:3278855/93‑1 (MQ=255)
                                           aTGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTa                                               <  2:5461828/93‑1 (MQ=255)
                                            tGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTAc                                              >  2:6677814/1‑93 (MQ=255)
                                            tGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTAc                                              >  2:2654517/1‑93 (MQ=255)
                                            tGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTAc                                              <  2:5571333/93‑1 (MQ=255)
                                             gcgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACt                                             <  1:6697276/93‑1 (MQ=255)
                                             gcgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACt                                             >  1:2955181/1‑93 (MQ=255)
                                              cgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTg                                            <  1:4664551/93‑1 (MQ=255)
                                               gctgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGc                                           <  1:3835327/93‑1 (MQ=255)
                                                 tgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCgg                                         <  1:148223/93‑1 (MQ=255)
                                                 tgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCgg                                         >  1:2655670/1‑93 (MQ=255)
                                                  gctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCggg                                        >  2:2377019/1‑93 (MQ=255)
                                                  gctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCggg                                        >  1:5323497/1‑93 (MQ=255)
                                                   ctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGt                                       <  2:7402388/93‑1 (MQ=255)
                                                   ctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGt                                       >  2:2440134/1‑93 (MQ=255)
                                                    tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTa                                      >  2:2550854/1‑93 (MQ=255)
                                                    tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTa                                      <  2:4425544/93‑1 (MQ=255)
                                                    tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTa                                      <  2:5868181/93‑1 (MQ=255)
                                                     cGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAt                                     <  1:2508681/93‑1 (MQ=255)
                                                      gATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                    >  1:3741855/1‑93 (MQ=255)
                                                      gATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                    <  1:5677614/93‑1 (MQ=255)
                                                      gATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                    >  2:6261834/1‑93 (MQ=255)
                                                      gATAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                    <  1:1093407/93‑1 (MQ=255)
                                                       aTAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTg                                   >  1:2148192/1‑93 (MQ=255)
                                                        tAAACTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGa                                  >  1:6427901/1‑93 (MQ=255)
                                                         aaaCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaa                                 >  1:5565235/1‑93 (MQ=255)
                                                         aaaCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaa                                 <  2:6760837/93‑1 (MQ=255)
                                                          aaCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaaa                                <  1:4797786/93‑1 (MQ=255)
                                                          aaCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaaa                                >  2:6122510/1‑93 (MQ=255)
                                                          aaCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaaa                                >  2:1817671/1‑93 (MQ=255)
                                                           aCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaaaa                               <  1:1943938/93‑1 (MQ=255)
                                                           aCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaaaa                               <  2:4770035/93‑1 (MQ=255)
                                                           aCTGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaaaa                               <  1:7419006/93‑1 (MQ=255)
                                                             tGGATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACt                             >  1:5871442/1‑93 (MQ=255)
                                                              ggATGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTg                            <  2:1647422/93‑1 (MQ=255)
                                                                 tGCATTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCa                         <  1:6770075/93‑1 (MQ=255)
                                                                    aTTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGa                      >  1:4831191/1‑93 (MQ=255)
                                                                    aTTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGa                      >  2:1779035/1‑93 (MQ=255)
                                                                    aTTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGa                      >  2:80982/1‑93 (MQ=255)
                                                                    aTTTCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGa                      <  2:1549338/93‑1 (MQ=255)
                                                                     tttCGCCACCTGACGCTGGACGATCACTGCCCGAGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGaa                     <  2:2165508/93‑1 (MQ=255)
                                                                      ttCGCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGaaa                    <  1:4234996/93‑1 (MQ=255)
                                                                       tCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAg                   <  1:5909897/93‑1 (MQ=255)
                                                                         gCCACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                 <  2:3311923/93‑1 (MQ=255)
                                                                          ccACCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTa                <  1:1258519/93‑1 (MQ=255)
                                                                            aCCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATc              >  2:429724/1‑93 (MQ=255)
                                                                            aCCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATc              <  1:7248326/93‑1 (MQ=255)
                                                                            aCCTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATc              >  2:2884981/1‑93 (MQ=255)
                                                                             ccTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCa             <  1:1795110/93‑1 (MQ=255)
                                                                              cTGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAg            <  2:3201163/93‑1 (MQ=255)
                                                                               tGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGc           >  2:5991035/1‑93 (MQ=255)
                                                                               tGACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGc           >  2:5224042/1‑93 (MQ=255)
                                                                                gACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCt          >  1:6117507/1‑93 (MQ=255)
                                                                                gACGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCt          >  1:5349638/1‑93 (MQ=255)
                                                                                 aCGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTg         <  1:7530220/93‑1 (MQ=255)
                                                                                  cGCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTgg        <  2:2781587/93‑1 (MQ=255)
                                                                                   gCTGGACGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGga       <  1:2839522/93‑1 (MQ=255)
                                                                                       gACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaagaa   <  2:3813061/93‑1 (MQ=255)
                                                                                        aCGATCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGaagaag  <  1:6175029/93‑1 (MQ=255)
                                                                                           |                                                                                         
CTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAG  >  NC_000913/3842137‑3842317

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