Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,227,795 A→G intergenic (‑74/‑138) pepE ← / → yjbC peptidase E/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,7950AG99.4% 1591.4 / ‑4.5 503intergenic (‑74/‑138)pepE/yjbCpeptidase E/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (262/236);  total (264/239)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 566 total aligned reads displayed.

ATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACTATTTTT  >  NC_000913/4227706‑4227886
                                                                                         |                                                                                           
aTAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGa                                                                                          >  1:8707239/1‑93 (MQ=255)
 tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGAt                                                                                         <  2:8179552/93‑1 (MQ=255)
 tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGAt                                                                                         >  1:2589818/1‑93 (MQ=255)
 tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCCCGAt                                                                                         >  1:2499962/1‑93 (MQ=255)
   aaaGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGt                                                                                       <  2:3908593/93‑1 (MQ=255)
   aaaGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGt                                                                                       >  2:492844/1‑93 (MQ=255)
       cAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGca                                                                                   <  1:8567398/93‑1 (MQ=255)
        aGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGgac                                                                                  >  1:6173700/1‑93 (MQ=255)
        aGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGcac                                                                                  <  1:8048879/93‑1 (MQ=255)
         gTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGcaca                                                                                 <  1:6039597/93‑1 (MQ=255)
         gTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGcaca                                                                                 <  1:4862219/93‑1 (MQ=255)
         gTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCGcaca                                                                                 >  1:9820714/1‑93 (MQ=255)
          ttCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACAt                                                                                <  2:4700552/93‑1 (MQ=255)
          ttCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACAt                                                                                >  1:1999824/1‑93 (MQ=255)
          ttCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACAt                                                                                >  1:6287393/1‑93 (MQ=255)
            ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCCCATTg                                                                              >  1:5027406/1‑93 (MQ=255)
            ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTg                                                                              <  2:4161050/93‑1 (MQ=255)
             cATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGGTGTCGCACATTGc                                                                             >  2:377116/1‑93 (MQ=255)
             cATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGc                                                                             >  2:5802422/1‑93 (MQ=255)
             cATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGc                                                                             <  2:8939614/93‑1 (MQ=255)
              aTTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCt                                                                            >  2:6204747/1‑93 (MQ=255)
              aTTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCt                                                                            <  2:4424223/93‑1 (MQ=255)
               ttACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCtt                                                                           <  1:9554398/93‑1 (MQ=255)
               ttACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATCTCGCACATTGCtt                                                                           >  1:4398981/1‑93 (MQ=255)
                tACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTc                                                                          <  2:1169190/93‑1 (MQ=255)
                tACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTc                                                                          <  2:1343479/93‑1 (MQ=255)
                tACTTATCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTc                                                                          <  1:9223654/93‑1 (MQ=255)
                 aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCt                                                                         >  1:3633920/1‑93 (MQ=255)
                  cTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCCCATTGCTTCTg                                                                        >  1:5108946/1‑93 (MQ=255)
                   ttCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTgg                                                                       <  2:7133839/93‑1 (MQ=255)
                   ttCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTgg                                                                       <  2:8781815/93‑1 (MQ=255)
                    tctcCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGa                                                                      <  2:6737053/93‑1 (MQ=255)
                      tcCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGaaa                                                                    <  2:3199297/93‑1 (MQ=255)
                      tcCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGaaa                                                                    >  1:5506206/1‑93 (MQ=255)
                       ccGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAg                                                                   >  2:2043461/1‑93 (MQ=255)
                         ggTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTg                                                                 >  1:1202991/1‑93 (MQ=255)
                         ggTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTg                                                                 >  2:7780716/1‑93 (MQ=255)
                          gTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGa                                                                >  2:5465830/1‑93 (MQ=255)
                          gTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGa                                                                >  2:594467/1‑93 (MQ=255)
                           ttAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGAt                                                               <  2:2459038/93‑1 (MQ=255)
                            tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGGAAGAGAtt                                                              >  1:6857950/1‑93 (MQ=255)
                            tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGAtt                                                              <  2:9499625/93‑1 (MQ=255)
                              gTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTg                                                            >  1:1003373/1‑93 (MQ=255)
                               tAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGa                                                           <  1:270299/93‑1 (MQ=255)
                               tAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGa                                                           <  1:1798036/93‑1 (MQ=255)
                                 aGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGc                                                         <  2:4314265/93‑1 (MQ=255)
                                 aGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGc                                                         <  2:5024556/93‑1 (MQ=255)
                                  gAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCa                                                        <  2:1773586/93‑1 (MQ=255)
                                  gAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCa                                                        >  2:7987324/1‑93 (MQ=255)
                                  gAAATCACGGGAAGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGAAGTCGCACATTGCTTCTGGAAAGTGATTTGAGcg                                                        >  2:2195673/1‑92 (MQ=255)
                                     aTCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGt                                                     <  1:9348299/93‑1 (MQ=255)
                                     aTCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGt                                                     <  1:4525444/93‑1 (MQ=255)
                                      tCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTa                                                    <  1:7504272/93‑1 (MQ=255)
                                      acACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATCGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTa                                                    <  1:5233709/92‑1 (MQ=255)
                                       cACGGGAGGCATGCTACCACTGCGAAGTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTaa                                                   <  1:8446038/93‑1 (MQ=255)
                                         cGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACg                                                 >  1:1665304/1‑93 (MQ=255)
                                          gggAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGa                                                <  1:2374374/93‑1 (MQ=255)
                                           ggAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGaa                                               >  2:2339686/1‑93 (MQ=255)
                                               gCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCa                                           <  2:8491051/93‑1 (MQ=255)
                                                cATGCTCCCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCaa                                          <  1:8859223/93‑1 (MQ=255)
                                                  tGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGc                                        >  2:1007101/1‑93 (MQ=255)
                                                  tGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGc                                        >  2:2999067/1‑93 (MQ=255)
                                                    cTACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCa                                      <  1:1518905/93‑1 (MQ=255)
                                                     tACCACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCaa                                     >  1:32720/1‑93 (MQ=255)
                                                       ccACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATg                                   <  2:1534569/93‑1 (MQ=255)
                                                        cACTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATgg                                  <  1:2788055/93‑1 (MQ=255)
                                                          cagCGAATTTAACCACTTACTGACATGGCTCTCGATGCCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaa                                <  1:6617794/91‑1 (MQ=255)
                                                          cTGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaa                                <  2:6331604/93‑1 (MQ=255)
                                                           tGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCGCATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaaa                               >  2:6085448/1‑93 (MQ=255)
                                                           tGCGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGaaa                               >  2:7569113/1‑93 (MQ=255)
                                                             cGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATa                             <  2:9830840/93‑1 (MQ=255)
                                                             cGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATa                             <  2:9143404/93‑1 (MQ=255)
                                                             cGAATTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATa                             <  1:4246327/93‑1 (MQ=255)
                                                               aaTTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAAc                           <  1:4070337/93‑1 (MQ=255)
                                                                aTTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGTTAACGAATTCAAGCCAATGGAAATAACa                          <  2:6439932/93‑1 (MQ=255)
                                                                aTTTAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACa                          >  1:3157870/1‑93 (MQ=255)
                                                                 tttAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACaa                         <  2:2770975/93‑1 (MQ=255)
                                                                  ttAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACaaa                        <  1:5746155/93‑1 (MQ=255)
                                                                   tAACCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAAt                       <  2:7229920/93‑1 (MQ=255)
                                                                    aaCCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATa                      >  1:3420441/1‑93 (MQ=255)
                                                                    aaCCACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATa                      <  1:3841287/93‑1 (MQ=255)
                                                                       cACTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAAAAAATAGtt                   >  2:7474110/1‑93 (MQ=255)
                                                                        aCTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTc                  >  2:1682295/1‑93 (MQ=255)
                                                                         cTTACTGACATGGCTCGCTATGTCGCACATTGTTTCGGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCg                 <  2:6839720/93‑1 (MQ=255)
                                                                         cTTACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTcc                 >  1:2113156/1‑92 (MQ=255)
                                                                          ttACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGc                >  2:6505341/1‑93 (MQ=255)
                                                                          ttACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGc                >  1:7020496/1‑93 (MQ=255)
                                                                          ttACTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGc                <  2:9408431/93‑1 (MQ=255)
                                                                             cTGACATGGCTGGCGATGTCGCAGATTGCTCCTGGAAAGTCATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCa             <  1:4734658/93‑1 (MQ=255)
                                                                             cTGACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCa             >  1:8209689/1‑93 (MQ=255)
                                                                               gACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATg           <  2:8298707/93‑1 (MQ=255)
                                                                               gACATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATg           <  1:8901426/93‑1 (MQ=255)
                                                                                aCATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGa          >  1:3036820/1‑93 (MQ=255)
                                                                                 cATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGAc         <  2:9743438/93‑1 (MQ=255)
                                                                                 cATGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGAc         <  2:3379377/93‑1 (MQ=255)
                                                                                  aTGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACt        >  2:5728801/1‑93 (MQ=255)
                                                                                  aTGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACt        <  1:7152361/93‑1 (MQ=255)
                                                                                   tGGCTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACTa       >  1:2264500/1‑93 (MQ=255)
                                                                                      cTCGCGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACTAttt    >  2:5619290/1‑93 (MQ=255)
                                                                                        cACGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACTAttttt  <  1:7719682/93‑1 (MQ=255)
                                                                                         |                                                                                           
ATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCGCACATTGCTTCTGGAAAGTGATTTGAGCAGGTAACGAATTCAAGCCAATGGAAATAACAAATAGTTCGCTCATGACTATTTTT  >  NC_000913/4227706‑4227886

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