Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 981,066 A→T D266V (GAC→GTC)  ycbB → putative amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913981,0660AT99.6% 821.9 / NA 236D266V (GAC→GTC) ycbBputative amidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base T (117/117);  minor base G (1/0);  total (118/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

AACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAA  >  NC_000913/980976‑981155
                                                                                          |                                                                                         
aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                         <  1:3380883/93‑1 (MQ=255)
 aCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCgg                                                                                        <  2:6791388/93‑1 (MQ=255)
 aCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCgg                                                                                        >  2:441893/1‑93 (MQ=255)
   tGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCgggg                                                                                      <  1:7683176/93‑1 (MQ=255)
   tGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCgggg                                                                                      <  2:2030199/93‑1 (MQ=255)
    gACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCggggg                                                                                     >  1:2354217/1‑93 (MQ=255)
     aCCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCgggggg                                                                                    >  2:6315384/1‑93 (MQ=255)
      ccGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGc                                                                                   <  2:827021/93‑1 (MQ=255)
         gCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTAGGGGGGCCGa                                                                                >  2:4941864/1‑93 (MQ=255)
          cAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGaa                                                                               <  1:4240467/93‑1 (MQ=255)
              gCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAAtt                                                                           <  2:3443128/93‑1 (MQ=255)
                aaCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTAc                                                                         >  1:6129141/1‑93 (MQ=255)
                 aCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTAcc                                                                        >  2:6392666/1‑92 (MQ=255)
                   gTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActct                                                                      >  2:1586872/1‑93 (MQ=255)
                   gTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActct                                                                      >  1:8083685/1‑93 (MQ=255)
                    ttGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTa                                                                     >  1:355432/1‑93 (MQ=255)
                    ttGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTa                                                                     >  2:3913549/1‑93 (MQ=255)
                     tGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTAc                                                                    >  2:5075784/1‑93 (MQ=255)
                      gcgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGGCGAAAATTACTCTAcc                                                                   >  2:7520243/1‑93 (MQ=255)
                       cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCt                                                                  <  1:5863153/93‑1 (MQ=255)
                       cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCt                                                                  <  2:4230006/93‑1 (MQ=255)
                       cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCt                                                                  <  2:6034765/93‑1 (MQ=255)
                        gcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTg                                                                 >  1:3901233/1‑93 (MQ=255)
                           cAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCg                                                              <  1:2570873/93‑1 (MQ=255)
                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                             <  1:7564529/93‑1 (MQ=255)
                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                             <  2:8106958/93‑1 (MQ=255)
                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                             <  2:573874/93‑1 (MQ=255)
                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                             <  1:1139239/93‑1 (MQ=255)
                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                             >  1:1194335/1‑93 (MQ=255)
                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAACTACTCTACCTCGCGa                                                             >  2:3825831/1‑93 (MQ=255)
                             gggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGAt                                                            >  2:3246849/1‑93 (MQ=255)
                             gggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGAt                                                            <  2:8687127/93‑1 (MQ=255)
                              ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATg                                                           <  1:7776107/93‑1 (MQ=255)
                              ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATg                                                           <  2:6893463/93‑1 (MQ=255)
                               gCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGa                                                          <  2:7190742/93‑1 (MQ=255)
                                cAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCgacgac                                                         <  2:257137/93‑1 (MQ=255)
                                cAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGac                                                         <  1:6856649/93‑1 (MQ=255)
                                cAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGac                                                         >  2:1431364/1‑93 (MQ=255)
                                cAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGac                                                         <  2:6542108/93‑1 (MQ=255)
                                 aGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGaca                                                        >  2:2899545/1‑93 (MQ=255)
                                 aGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACATGGCGATGaca                                                        >  1:2911276/1‑93 (MQ=255)
                                  gtagAGTAACGACGTACCGGCGTTGCGCGAAATATTACAACGCACAGACATGTTGATCGGGGGGCCGAAAAATACTCTACATGGCGATGacac                                                       <  2:7815137/90‑1 (MQ=255)
                                  gTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGacac                                                       >  1:5481917/1‑93 (MQ=255)
                                  gTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGacac                                                       <  2:1156580/93‑1 (MQ=255)
                                   tGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACg                                                      >  1:1827650/1‑93 (MQ=255)
                                   cggAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACg                                                      <  2:5110080/92‑1 (MQ=255)
                                     gAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGcc                                                    >  2:5859117/1‑93 (MQ=255)
                                        tAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAAc                                                 <  1:4883547/93‑1 (MQ=255)
                                           cgacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGa                                              <  1:7117833/93‑1 (MQ=255)
                                           cgacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGa                                              <  1:633806/93‑1 (MQ=255)
                                             acgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcg                                            <  2:6218310/93‑1 (MQ=255)
                                                  ccGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTa                                       >  1:70386/1‑93 (MQ=255)
                                                   cGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAg                                      <  2:4611559/93‑1 (MQ=255)
                                                   cGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAg                                      <  1:8641338/93‑1 (MQ=255)
                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    <  2:2175382/93‑1 (MQ=255)
                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    <  2:5338037/93‑1 (MQ=255)
                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    >  1:8339640/1‑93 (MQ=255)
                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    <  2:7673023/93‑1 (MQ=255)
                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    <  1:1550531/93‑1 (MQ=255)
                                                      cGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCa                                   <  2:3992760/93‑1 (MQ=255)
                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCac                                  >  2:2539868/1‑92 (MQ=255)
                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAg                                  >  2:8477406/1‑93 (MQ=255)
                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAg                                  >  2:5539605/1‑93 (MQ=255)
                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAg                                  >  2:1023466/1‑93 (MQ=255)
                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCACACTGACGCGGTAGTCAg                                  >  2:1062153/1‑93 (MQ=255)
                                                             cgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                            <  2:2797550/93‑1 (MQ=255)
                                                             cgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                            <  1:5191752/93‑1 (MQ=255)
                                                              gAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTAAGCCCATcc                           >  2:5669795/1‑93 (MQ=255)
                                                              gAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGCCAGCCCATcc                           >  2:1833394/1‑93 (MQ=255)
                                                               aaaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCg                          >  1:7278839/1‑93 (MQ=255)
                                                                aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                         <  1:5063502/93‑1 (MQ=255)
                                                                aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                         <  2:4691897/93‑1 (MQ=255)
                                                                aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                         <  1:2179737/93‑1 (MQ=255)
                                                                aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                         >  1:8106393/1‑93 (MQ=255)
                                                                aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                         <  1:2526387/93‑1 (MQ=255)
                                                                 atatTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCt                        <  2:191666/93‑1 (MQ=255)
                                                                   atTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGt                      <  2:7632986/93‑1 (MQ=255)
                                                                    ttGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGtt                     <  2:597171/93‑1 (MQ=255)
                                                                    ttGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGtt                     >  1:330053/1‑93 (MQ=255)
                                                                    ttGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGtt                     <  1:5808990/93‑1 (MQ=255)
                                                                      gCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGCCGCGGTAGTCAGCCCATCCGCTGTTag                   >  2:5275901/1‑92 (MQ=255)
                                                                      gCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTAc                   <  1:4042256/93‑1 (MQ=255)
                                                                      gCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTAc                   >  2:4924473/1‑93 (MQ=255)
                                                                       cAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACt                  <  1:6948526/93‑1 (MQ=255)
                                                                        aaCGCACAGGCATGTTGGTCGGGGGGCCGAAAATTCCTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTg                 <  2:5411690/93‑1 (MQ=255)
                                                                         aCGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGt                <  2:3592396/93‑1 (MQ=255)
                                                                          cGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGtt               <  2:314742/93‑1 (MQ=255)
                                                                          cGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGtt               <  1:493425/93‑1 (MQ=255)
                                                                           gCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTg              >  2:8745191/1‑93 (MQ=255)
                                                                            cacaGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGa             <  2:7073311/93‑1 (MQ=255)
                                                                             acaGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGaa            >  1:3067019/1‑93 (MQ=255)
                                                                               aGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAc          >  1:5296036/1‑93 (MQ=255)
                                                                               aGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAc          >  2:560311/1‑93 (MQ=255)
                                                                                ggCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAca         >  2:8172467/1‑93 (MQ=255)
                                                                                  cATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagc       <  2:4812897/93‑1 (MQ=255)
                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      >  1:411110/1‑93 (MQ=255)
                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      <  1:6484685/93‑1 (MQ=255)
                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      >  2:2321496/1‑93 (MQ=255)
                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGCAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      >  2:8062360/1‑93 (MQ=255)
                                                                                    tGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagcag     <  2:7256927/93‑1 (MQ=255)
                                                                                          |                                                                                         
AACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAA  >  NC_000913/980976‑981155

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