Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,339,624 A→G K348R (AAG→AGG)  yciM → putative heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,339,6240AG98.4% 401.7 / ‑4.9 126K348R (AAG→AGG) yciMputative heat shock protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (73/51);  total (75/51)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

TCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGT  >  NC_000913/1339534‑1339713
                                                                                          |                                                                                         
tCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAggg                                                                                         <  1:3384748/93‑1 (MQ=255)
tCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAggg                                                                                         <  1:805431/93‑1 (MQ=255)
  cATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTa                                                                                       >  2:4334852/1‑93 (MQ=255)
  cATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTa                                                                                       <  1:138238/93‑1 (MQ=255)
   aTAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTAc                                                                                      >  1:463587/1‑93 (MQ=255)
   aTAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTAc                                                                                      <  1:1446164/93‑1 (MQ=255)
     aaGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgt                                                                                    >  2:481453/1‑93 (MQ=255)
      aGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgta                                                                                   <  2:3753242/93‑1 (MQ=255)
       gTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtag                                                                                  >  1:1029646/1‑93 (MQ=255)
         tAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtagta                                                                                <  2:4354021/93‑1 (MQ=255)
          aaTGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTaa                                                                               >  2:4089989/1‑93 (MQ=255)
               aTTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTc                                                                          >  1:3740675/1‑93 (MQ=255)
               aTTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTc                                                                          <  2:2097203/93‑1 (MQ=255)
                 tACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGt                                                                        >  1:3186379/1‑93 (MQ=255)
                 tACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGt                                                                        >  1:2047454/1‑93 (MQ=255)
                 tACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGt                                                                        <  1:4725942/93‑1 (MQ=255)
                    cACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTAt                                                                     >  1:1030037/1‑93 (MQ=255)
                    cACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTAt                                                                     <  2:3439483/93‑1 (MQ=255)
                      cTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCg                                                                   >  1:2838054/1‑93 (MQ=255)
                      cTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCg                                                                   >  1:317562/1‑93 (MQ=255)
                       ttAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGt                                                                  >  2:4499839/1‑93 (MQ=255)
                        tAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGtt                                                                 <  1:1537712/93‑1 (MQ=255)
                         aaaTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTg                                                                >  2:4693938/1‑93 (MQ=255)
                         aaaTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTg                                                                >  2:1291187/1‑93 (MQ=255)
                          aaTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGc                                                               <  2:738617/93‑1 (MQ=255)
                          aaTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGc                                                               >  1:4110906/1‑93 (MQ=255)
                           aTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGcc                                                              <  1:4394024/93‑1 (MQ=255)
                            tGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCa                                                             >  1:3386341/1‑93 (MQ=255)
                              aaGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGa                                                           <  2:487818/93‑1 (MQ=255)
                               aGCGGAAGAAGGGCGTGCCAAATAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaa                                                          >  1:2433060/1‑93 (MQ=255)
                                gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         <  1:2643566/93‑1 (MQ=255)
                                 cGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAt                                                        >  1:3607990/1‑93 (MQ=255)
                                  ggAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtg                                                       <  2:136183/93‑1 (MQ=255)
                                  ggAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtg                                                       <  2:136184/93‑1 (MQ=255)
                                   gaagaaGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgt                                                      >  1:1139794/1‑93 (MQ=255)
                                   gaagaaGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgt                                                      >  1:1139795/1‑93 (MQ=255)
                                    aagaagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgtg                                                     <  1:4353144/93‑1 (MQ=255)
                                     agaagGGCGTACCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTgg                                                    >  2:820744/1‑93 (MQ=255)
                                         gggCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtttt                                                >  2:2059465/1‑93 (MQ=255)
                                         gggCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtttt                                                >  2:3671032/1‑93 (MQ=255)
                                         gggCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtttt                                                <  1:1388636/93‑1 (MQ=255)
                                             gTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCg                                            <  1:1629662/93‑1 (MQ=255)
                                              tGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGc                                           <  1:4513792/93‑1 (MQ=255)
                                               gCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCa                                          <  2:2080600/93‑1 (MQ=255)
                                               gCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCa                                          >  2:469416/1‑93 (MQ=255)
                                               gCCAAAGAGAGCCTGATGGGGCTCCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTCCCAGAAATGTGGTTTTACCGCa                                          <  2:2006740/93‑1 (MQ=255)
                                                ccAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCAt                                         >  2:3378697/1‑93 (MQ=255)
                                                 cAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATa                                        <  2:2501024/93‑1 (MQ=255)
                                                 cAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATa                                        >  2:2752659/1‑93 (MQ=255)
                                                 cAAAGAAAGCCTGATGATGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATa                                        >  2:3130824/1‑93 (MQ=255)
                                                  aaaGAGAGCCTGATGGTGCTGCGTGACATGTTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATac                                       >  2:486269/1‑93 (MQ=255)
                                                  aaaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATac                                       >  2:2374670/1‑93 (MQ=255)
                                                   aaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATaca                                      <  1:3481234/93‑1 (MQ=255)
                                                    agagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATacac                                     >  1:4627563/1‑93 (MQ=255)
                                                    agagagCCTGATGGTGCTGCGGGACATGGTTGGCGAGAGGGAACGTAGCAAGCCTCGTTATCGTGGCCAGAAATGTGGTTTTACCGCCTaacc                                     >  1:1230289/1‑90 (MQ=255)
                                                     gagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAcc                                    >  2:1425946/1‑93 (MQ=255)
                                                     gagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAcc                                    >  2:3009491/1‑93 (MQ=255)
                                                      agagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAccc                                   <  2:4351702/93‑1 (MQ=255)
                                                      agagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAccc                                   >  1:3720725/1‑93 (MQ=255)
                                                      agagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAccc                                   <  1:3821697/93‑1 (MQ=255)
                                                       gagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCct                                  >  1:48929/1‑93 (MQ=255)
                                                       gagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCct                                  >  1:409743/1‑93 (MQ=255)
                                                       gagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGGACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCct                                  >  2:2398851/1‑93 (MQ=255)
                                                          ccTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTa                               <  2:4595665/93‑1 (MQ=255)
                                                           cTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTAc                              >  2:2680397/1‑93 (MQ=255)
                                                            tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             >  1:3516004/1‑93 (MQ=255)
                                                            tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             >  1:394603/1‑93 (MQ=255)
                                                            tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             <  2:749236/93‑1 (MQ=255)
                                                              aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgg                           >  1:3003826/1‑93 (MQ=255)
                                                                ggTGCTGCGTGACATGGTGGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCa                         <  1:1311892/93‑1 (MQ=255)
                                                                 gTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAt                        >  1:2709654/1‑93 (MQ=255)
                                                                 gTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAt                        <  1:1300645/93‑1 (MQ=255)
                                                                  tgctgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAtt                       >  2:4085960/1‑93 (MQ=255)
                                                                  tgctgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAtt                       >  2:4085961/1‑93 (MQ=255)
                                                                     tgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTc                    <  1:4607308/93‑1 (MQ=255)
                                                                     tgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTc                    >  2:4524938/1‑93 (MQ=255)
                                                                      gcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTcc                   >  2:1000918/1‑93 (MQ=255)
                                                                      gcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTcc                   >  2:1000919/1‑93 (MQ=255)
                                                                        gTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGt                 >  2:2558414/1‑93 (MQ=255)
                                                                        gTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGt                 >  1:584519/1‑93 (MQ=255)
                                                                         tGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTc                >  1:1704864/1‑93 (MQ=255)
                                                                         tGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTc                >  2:1583169/1‑93 (MQ=255)
                                                                         tGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCAGTc                <  1:3340952/93‑1 (MQ=255)
                                                                         tGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGAGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGtt                >  2:1865319/1‑92 (MQ=255)
                                                                          gACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCt               >  2:234200/1‑93 (MQ=255)
                                                                           aCATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCtt              >  2:94136/1‑93 (MQ=255)
                                                                           aCATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCtt              <  2:1621230/93‑1 (MQ=255)
                                                                           aCATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCtt              <  2:189143/93‑1 (MQ=255)
                                                                             aTGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGt            <  2:1926686/93‑1 (MQ=255)
                                                                              tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           <  1:252214/93‑1 (MQ=255)
                                                                              tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           >  1:1991339/1‑93 (MQ=255)
                                                                                   ggCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGcc      <  2:2966919/93‑1 (MQ=255)
                                                                                    gCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCt     <  2:1124092/93‑1 (MQ=255)
                                                                                    gCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCt     <  1:998874/93‑1 (MQ=255)
                                                                                    gCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCt     <  1:201533/93‑1 (MQ=255)
                                                                                    gCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCt     <  1:201532/93‑1 (MQ=255)
                                                                                     cGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTg    <  1:738562/93‑1 (MQ=255)
                                                                                      gagaGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTgg   >  2:3337958/1‑93 (MQ=255)
                                                                                       agagGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGt  >  2:1150244/1‑93 (MQ=255)
                                                                                       agagGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGt  >  1:837539/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGT  >  NC_000913/1339534‑1339713

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