Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 28,729 C→T A119V (GCC→GTC)  dapB → dihydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091328,7290CT86.8% 719.9 / 64.7 250A119V (GCC→GTC) dapBdihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (103/114);  total (121/129)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

Only 100 of 250 total aligned reads displayed.

TCGCTTTTTGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGA  >  NC_000913/28639‑28819
                                                                                          |                                                                                          
tcGCTTTTTGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCg                                                                                          >  1:5615164/1‑93 (MQ=255)
    tttttGTCGCCAGCATGGCATCGGGAGGATGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGatct                                                                                      >  1:383197/1‑91 (MQ=255)
     ttttGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATAtt                                                                                     <  1:1603252/93‑1 (MQ=255)
     ttttGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATAtt                                                                                     >  2:6899210/1‑93 (MQ=255)
       ttGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGCTATTGc                                                                                   <  2:7097763/93‑1 (MQ=255)
        tGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCg                                                                                  >  1:3570283/1‑93 (MQ=255)
         gTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGa                                                                                 <  1:1145053/93‑1 (MQ=255)
         gTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGa                                                                                 <  1:3418689/93‑1 (MQ=255)
         gTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGa                                                                                 <  1:3359810/93‑1 (MQ=255)
         gTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGa                                                                                 >  2:3879720/1‑93 (MQ=255)
          tCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGAt                                                                                <  1:2574422/93‑1 (MQ=255)
          tCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGAt                                                                                >  2:1079029/1‑93 (MQ=255)
           cGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGAtt                                                                               >  2:4120589/1‑93 (MQ=255)
             ccAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGt                                                                             >  1:1005304/1‑93 (MQ=255)
              cAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTc                                                                            <  1:4676992/93‑1 (MQ=255)
              cAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTc                                                                            >  1:1821477/1‑93 (MQ=255)
                gCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCtt                                                                          <  1:6548044/93‑1 (MQ=255)
                  aTGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTtg                                                                        >  2:6022917/1‑93 (MQ=255)
                   tGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTtgc                                                                       >  1:7001800/1‑93 (MQ=255)
                        aaGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCa                                                                  >  1:2676256/1‑93 (MQ=255)
                        aaGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCa                                                                  >  1:91896/1‑93 (MQ=255)
                        aaGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCa                                                                  >  1:3100902/1‑93 (MQ=255)
                         aGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCaa                                                                 <  2:4634369/93‑1 (MQ=255)
                         aGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCaa                                                                 <  1:5022832/93‑1 (MQ=255)
                          gggATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAAt                                                                <  2:1647297/93‑1 (MQ=255)
                          gggATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAAt                                                                <  2:5712354/93‑1 (MQ=255)
                           ggATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAAtt                                                               <  2:4394528/93‑1 (MQ=255)
                            gATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAAttt                                                              >  2:5414199/1‑93 (MQ=255)
                              tGGCGATCGGCACTACGGGGCTTGACGAAGCCGGTAAACAAGCAATTTGTGACGCCGCTGTCGATAGTGAGATTGTCTCTGCTGCAAATTTTa                                                            >  1:5729331/1‑93 (MQ=255)
                               ggTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAg                                                           <  1:6440686/93‑1 (MQ=255)
                                  gATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGt                                                        <  2:4989493/93‑1 (MQ=255)
                                   aTCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGtt                                                       <  1:5169539/93‑1 (MQ=255)
                                    tCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCGATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTg                                                      <  2:7220960/93‑1 (MQ=255)
                                    tCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTg                                                      <  1:4231313/93‑1 (MQ=255)
                                    tCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTg                                                      >  2:577823/1‑93 (MQ=255)
                                    tCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTg                                                      <  2:1212564/93‑1 (MQ=255)
                                        cACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGt                                                  <  1:2700425/93‑1 (MQ=255)
                                        cACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTGTGCTGCCTATTTTAGCGTTGGCGt                                                  <  2:7362796/93‑1 (MQ=255)
                                         aCTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGtt                                                 >  2:6713350/1‑93 (MQ=255)
                                          cTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTa                                                <  2:4202794/93‑1 (MQ=255)
                                           gaCGGGTTTTGACGAAGCCGATAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTaa                                               <  2:2181342/92‑1 (MQ=255)
                                            aCGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAAc                                              >  2:381114/1‑93 (MQ=255)
                                            aCGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAAc                                              <  1:836592/93‑1 (MQ=255)
                                            aCGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTATTTGCTGCCAATTTTAGCGTTGGCGTTAAc                                              >  1:6657741/1‑93 (MQ=255)
                                               gggTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTc                                           <  1:5809572/93‑1 (MQ=255)
                                                ggTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCa                                          >  1:6223717/1‑93 (MQ=255)
                                                 gTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCAt                                         <  2:3422391/93‑1 (MQ=255)
                                                  tttGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATg                                        <  2:138999/93‑1 (MQ=255)
                                                    taaCGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCt                                      <  1:6376415/91‑1 (MQ=255)
                                                    tGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCt                                      >  2:1342650/1‑93 (MQ=255)
                                                     gACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCtt                                     >  1:7250792/1‑93 (MQ=255)
                                                       cGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTaa                                   >  2:6367775/1‑93 (MQ=255)
                                                       cGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTaa                                   <  2:4281614/93‑1 (MQ=255)
                                                        gAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAg                                  <  1:8601/93‑1 (MQ=255)
                                                         aaGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAgc                                 >  2:4761265/1‑93 (MQ=255)
                                                          aGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAgct                                >  2:194262/1‑93 (MQ=255)
                                                          aGCCGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAgct                                >  1:4653827/1‑93 (MQ=255)
                                                            ccGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAgctgc                              <  1:1790643/93‑1 (MQ=255)
                                                             cGGTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAgctgct                             >  2:2917458/1‑93 (MQ=255)
                                                             cGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAgctgct                             <  2:3971826/93‑1 (MQ=255)
                                                              ggTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGctgctg                            <  1:5370573/93‑1 (MQ=255)
                                                              ggTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGctgctg                            >  1:2579606/1‑93 (MQ=255)
                                                              ggTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGctgctg                            >  2:1884844/1‑93 (MQ=255)
                                                               gTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTgg                           >  2:6123591/1‑93 (MQ=255)
                                                               gTAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTgg                           >  2:7424285/1‑93 (MQ=255)
                                                               gTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTgg                           >  1:136643/1‑93 (MQ=255)
                                                               gTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTgg                           >  2:4613678/1‑93 (MQ=255)
                                                                tAAACAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGga                          >  1:1295122/1‑93 (MQ=255)
                                                                  aaCAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGgaga                        >  2:3302158/1‑93 (MQ=255)
                                                                  aaCAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGgaga                        >  2:5137349/1‑93 (MQ=255)
                                                                  aaCAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGgaga                        <  1:360886/93‑1 (MQ=255)
                                                                   aCAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGaa                       >  1:440476/1‑93 (MQ=255)
                                                                   aCAAGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGaa                       >  1:1970525/1‑93 (MQ=255)
                                                                     aaGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAag                     >  2:7933/1‑93 (MQ=255)
                                                                      aGCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAagc                    <  1:6084454/93‑1 (MQ=255)
                                                                       gCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAagca                   <  1:4218212/93‑1 (MQ=255)
                                                                       gCAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAagca                   <  2:5477070/93‑1 (MQ=255)
                                                                        cAATTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAagcag                  <  2:3198242/93‑1 (MQ=255)
                                                                         aaTTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAagcagc                 >  1:4763859/1‑93 (MQ=255)
                                                                         aaTTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAtgcagc                 >  1:4844365/1‑93 (MQ=255)
                                                                          aTTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGcc                <  1:7421083/93‑1 (MQ=255)
                                                                          aTTCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGcc                >  1:2866207/1‑93 (MQ=255)
                                                                           ttCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCa               >  1:594846/1‑93 (MQ=255)
                                                                           ttCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCa               <  1:5853734/93‑1 (MQ=255)
                                                                           ttCGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCa               >  1:5295010/1‑93 (MQ=255)
                                                                             cGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCaaa             <  2:1408041/93‑1 (MQ=255)
                                                                             cGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCaaa             <  2:7038908/93‑1 (MQ=255)
                                                                             cGTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCaaa             <  1:3912836/93‑1 (MQ=255)
                                                                              gTGACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAg            >  1:158112/1‑93 (MQ=255)
                                                                                gACGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTg          >  1:1059288/1‑93 (MQ=255)
                                                                                 aCGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGa         <  1:964088/93‑1 (MQ=255)
                                                                                 aCGCCGCTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGa         <  1:2074003/93‑1 (MQ=255)
                                                                                  cgccgcTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGAt        <  1:2247984/93‑1 (MQ=255)
                                                                                   gccgcTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATg       >  2:2455251/1‑93 (MQ=255)
                                                                                   gccgcTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATg       >  1:7318490/1‑93 (MQ=255)
                                                                                     cgcTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATggg     <  1:1472947/93‑1 (MQ=255)
                                                                                      gcTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGt    >  2:7523928/1‑93 (MQ=255)
                                                                                       cTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTg   >  2:3142807/1‑93 (MQ=255)
                                                                                       cTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTg   >  2:2663876/1‑93 (MQ=255)
                                                                                       cTGTCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTg   >  1:2211210/1‑93 (MQ=255)
                                                                                          |                                                                                          
TCGCTTTTTGTCGCCAGCATGGCAAAGGGATGGTGATCGGCACTACGGGGTTTGACGAAGCCGGTAAACAAGCAATTCGTGACGCCGCTGCCGATATTGCGATTGTCTTTGCTGCCAATTTTAGCGTTGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGA  >  NC_000913/28639‑28819

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