Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,978,702 C→T S68F (TCC→TTC)  yifK → putative amino acid/amine transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,978,7020CT99.5% 1477.1 / NA 409S68F (TCC→TTC) yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (206/201);  minor base G (0/2);  total (206/203)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 429 total aligned reads displayed.

TGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGT  >  NC_000913/3978612‑3978791
                                                                                          |                                                                                         
tGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTtcc                                                                                         <  2:1657181/93‑1 (MQ=255)
 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCt                                                                                        <  2:3014470/93‑1 (MQ=255)
 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCt                                                                                        <  1:3577550/93‑1 (MQ=255)
 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCt                                                                                        >  2:4354264/1‑93 (MQ=255)
 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCt                                                                                        >  2:4562228/1‑93 (MQ=255)
     tGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGac                                                                                    >  1:2716190/1‑92 (MQ=255)
     tGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                                    <  1:5206304/93‑1 (MQ=255)
      gggccgggccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAac                                                                                   <  1:100232/93‑1 (MQ=255)
          cgggccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggt                                                                               <  2:6159732/93‑1 (MQ=255)
           gggccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggtt                                                                              >  2:6681396/1‑93 (MQ=255)
             gccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttac                                                                            >  1:5801662/1‑93 (MQ=255)
             gccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttac                                                                            <  1:4520155/93‑1 (MQ=255)
             gccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttac                                                                            <  2:5439054/93‑1 (MQ=255)
             gccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttac                                                                            <  2:3947172/93‑1 (MQ=255)
              ccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttacc                                                                           <  2:2356313/93‑1 (MQ=255)
              ccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttacc                                                                           <  1:6962140/93‑1 (MQ=255)
              ccATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttacc                                                                           >  2:5877870/1‑93 (MQ=255)
               cATCCGTATTGTTGGCCTCTATCATCGCCGGGCTGTTCGTCTTTTTCATCATCAGTTCAATGGCCGCACTGTTCTTCCTCGAactggttaccg                                                                          >  1:5518800/1‑93 (MQ=255)
               cATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccg                                                                          <  1:2552959/93‑1 (MQ=255)
                aTCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccgg                                                                         >  1:128346/1‑93 (MQ=255)
                aTCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccgg                                                                         <  2:3338086/93‑1 (MQ=255)
                aTCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccgg                                                                         >  2:6399734/1‑93 (MQ=255)
                 tCCGTATTGTTGGCCTATATCATCTCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggt                                                                        <  2:4302179/93‑1 (MQ=255)
                 tCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggt                                                                        >  2:6050416/1‑93 (MQ=255)
                  ccGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                                       <  2:6563772/93‑1 (MQ=255)
                   cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttc                                                                      >  1:677732/1‑93 (MQ=255)
                    gTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcg                                                                     >  1:3507728/1‑93 (MQ=255)
                    gTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTAAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcg                                                                     <  1:2833385/93‑1 (MQ=255)
                     tATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgt                                                                    >  1:1945333/1‑93 (MQ=255)
                     tATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgt                                                                    >  2:2533685/1‑93 (MQ=255)
                     tATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgt                                                                    <  1:705934/93‑1 (MQ=255)
                      aTTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgtt                                                                   <  2:7865197/93‑1 (MQ=255)
                       ttgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgttc                                                                  >  1:7545141/1‑93 (MQ=255)
                         gttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGc                                                                >  2:3132758/1‑93 (MQ=255)
                         gtgggCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGc                                                                <  2:2817671/90‑1 (MQ=255)
                            ggCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGt                                                             <  2:7427880/93‑1 (MQ=255)
                            ggCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGt                                                             >  1:5886740/1‑93 (MQ=255)
                             gCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                            >  1:2205017/1‑93 (MQ=255)
                             gCCTATATCACCGCCGGGCTGTGCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                            <  1:6668785/93‑1 (MQ=255)
                              ccTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGttt                                                           <  2:4912919/93‑1 (MQ=255)
                              ccTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGttt                                                           <  1:7650331/93‑1 (MQ=255)
                               cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGTACCGCTTACCGGATCGTTCTCCGTTTa                                                          >  2:6656252/1‑93 (MQ=255)
                               cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                                          <  1:987459/93‑1 (MQ=255)
                               cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                                          <  1:1785125/93‑1 (MQ=255)
                               cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                                          >  1:2083232/1‑93 (MQ=255)
                               cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                                          <  1:7245464/93‑1 (MQ=255)
                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                         >  1:641879/1‑93 (MQ=255)
                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                         >  2:383552/1‑93 (MQ=255)
                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                         >  2:6126050/1‑93 (MQ=255)
                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                         >  2:262786/1‑93 (MQ=255)
                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                         >  2:2428807/1‑93 (MQ=255)
                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                                         <  2:1792245/93‑1 (MQ=255)
                                  tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                                       <  2:1890879/93‑1 (MQ=255)
                                  tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                                       <  2:5812337/93‑1 (MQ=255)
                                  tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                                       <  2:4678550/93‑1 (MQ=255)
                                     catcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCa                                                    <  1:3844199/93‑1 (MQ=255)
                                     catcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCa                                                    <  1:411787/93‑1 (MQ=255)
                                        ccccGGGTTGTTCGTCTTTTTCATAATGCGTTCAATGGGCAAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCg                                                 <  2:7767719/91‑1 (MQ=255)
                                         gCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGt                                                >  2:1458702/1‑93 (MQ=255)
                                            gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtat                                             >  1:1405729/1‑93 (MQ=255)
                                             ggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                                            >  2:701117/1‑93 (MQ=255)
                                             ggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                                            <  1:5695567/93‑1 (MQ=255)
                                              gCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTatat                                           >  2:7306442/1‑93 (MQ=255)
                                              gCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTatat                                           <  2:5075738/93‑1 (MQ=255)
                                                tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGa                                         >  2:5142792/1‑93 (MQ=255)
                                                tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGa                                         >  2:4795696/1‑93 (MQ=255)
                                                tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGa                                         >  2:7666522/1‑93 (MQ=255)
                                                   tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGcc                                      <  1:4897487/93‑1 (MQ=255)
                                                     gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                                    >  2:7590636/1‑93 (MQ=255)
                                                     gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                                    <  2:7149114/93‑1 (MQ=255)
                                                             catcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                            >  1:845013/1‑93 (MQ=255)
                                                              atcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGc                           <  2:2941462/93‑1 (MQ=255)
                                                               tcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCt                          <  2:2097491/93‑1 (MQ=255)
                                                                catGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTa                         <  2:3500147/93‑1 (MQ=255)
                                                                catGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTa                         >  2:6834741/1‑93 (MQ=255)
                                                                 atGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt                        >  1:4347189/1‑93 (MQ=255)
                                                                 atGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt                        >  1:4408319/1‑93 (MQ=255)
                                                                   gCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtct                      >  1:7441832/1‑93 (MQ=255)
                                                                     gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa                    >  2:3748649/1‑93 (MQ=255)
                                                                     gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa                    >  2:1349783/1‑93 (MQ=255)
                                                                      ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCAc                   >  1:1074759/1‑93 (MQ=255)
                                                                      ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCAc                   <  1:3155871/93‑1 (MQ=255)
                                                                         aaTGGGCGAAACGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGc                <  1:2446437/93‑1 (MQ=255)
                                                                          aTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGcc               >  1:6210742/1‑93 (MQ=255)
                                                                           tGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCt              <  1:7788955/93‑1 (MQ=255)
                                                                            gggCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTg             <  2:34637/93‑1 (MQ=255)
                                                                            gggCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTg             <  2:2899345/93‑1 (MQ=255)
                                                                            gggCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTg             <  2:722794/93‑1 (MQ=255)
                                                                                gAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCt         <  2:6300406/93‑1 (MQ=255)
                                                                                gAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCt         >  1:5311819/1‑93 (MQ=255)
                                                                                 aaaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCtt        >  1:3347159/1‑93 (MQ=255)
                                                                                  aaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTa       <  1:1603658/93‑1 (MQ=255)
                                                                                  aaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTa       <  2:1099486/93‑1 (MQ=255)
                                                                                   aTGTTGTTCCTCGAACCCGTTACCGGTTCGTTTGCCGGTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCTCCTTGtttttc      >  2:4230542/1‑91 (MQ=255)
                                                                                    ggttgttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACt     <  2:3181382/92‑1 (MQ=255)
                                                                                     gttgttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTg    >  1:5435274/1‑93 (MQ=255)
                                                                                      ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTgg   <  2:1990142/93‑1 (MQ=255)
                                                                                       tgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt  <  2:3653966/93‑1 (MQ=255)
                                                                                       tgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt  >  1:2942437/1‑93 (MQ=255)
                                                                                       tgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt  <  2:5561632/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGT  >  NC_000913/3978612‑3978791

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: