Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,386,365 T→A F169Y (TTC→TAC)  yfbL → putative M28A family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,3650TA97.5% 401.2 / ‑0.7 120F169Y (TTC→TAC) yfbLputative M28A family peptidase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (58/59);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 121 total aligned reads displayed.

GTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACG  >  NC_000913/2386269‑2386460
                                                                                                |                                                                                               
  gtgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTc                                                                                            >  2:1105075/1‑100 (MQ=255)
    ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGt                                                                                          >  2:528062/1‑100 (MQ=255)
      cAGGATTACTCGAACTGGCGCGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGtcg                                                                                        >  1:3172716/1‑98 (MQ=255)
      cAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAg                                                                                        <  2:260874/100‑1 (MQ=255)
      cAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAg                                                                                        <  2:1049482/100‑1 (MQ=255)
       aGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGc                                                                                       <  1:497741/100‑1 (MQ=255)
        ggATTACTCGAACTGGCACGTTTGTTACAGCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCg                                                                                      <  2:2684860/100‑1 (MQ=255)
          aTTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGAt                                                                                    <  1:3046741/100‑1 (MQ=255)
          aTTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGAt                                                                                    <  2:1146896/100‑1 (MQ=255)
          aTTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGAt                                                                                    <  1:2007293/100‑1 (MQ=255)
             aCTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGaa                                                                                 >  1:2478891/1‑100 (MQ=255)
              cTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGaaa                                                                                >  2:1213700/1‑100 (MQ=255)
              cTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGaaa                                                                                >  2:474664/1‑100 (MQ=255)
              cTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGaaa                                                                                >  2:2872062/1‑100 (MQ=255)
                cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGACCAGCCCCACTTTCGTAGCCATGaaagg                                                                              >  1:3156489/1‑98 (MQ=255)
                cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATg                                                                              <  2:2433524/100‑1 (MQ=255)
                cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATg                                                                              <  2:2433568/100‑1 (MQ=255)
                cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATg                                                                              >  2:743508/1‑100 (MQ=255)
                  aaCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATggg                                                                            <  2:1819601/100‑1 (MQ=255)
                  aaCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATggg                                                                            <  2:155662/100‑1 (MQ=255)
                   aCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATgggg                                                                           <  1:1856173/100‑1 (MQ=255)
                   aCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATgggg                                                                           <  2:2953788/100‑1 (MQ=255)
                   aCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATgggg                                                                           <  2:502292/100‑1 (MQ=255)
                     tGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAg                                                                         >  1:653768/1‑100 (MQ=255)
                     tGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAg                                                                         >  1:3393900/1‑100 (MQ=255)
                     tGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAg                                                                         <  1:3246628/100‑1 (MQ=255)
                      ggCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAgc                                                                        >  1:1759103/1‑100 (MQ=255)
                      ggCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAgc                                                                        >  1:1190160/1‑100 (MQ=255)
                            tttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtg                                                                  <  1:1600659/100‑1 (MQ=255)
                             ttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgc                                                                 <  1:900645/100‑1 (MQ=255)
                               gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                               >  1:1653435/1‑100 (MQ=255)
                               gTTACATCAGCAAGTACAGAAAACAGGCGTGCAGCTGGGCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgca                                                                <  1:1389969/99‑1 (MQ=255)
                                ttACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatg                                                              <  2:136796/100‑1 (MQ=255)
                                ttACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatg                                                              <  2:2828928/100‑1 (MQ=255)
                                 tACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                             <  2:2138102/100‑1 (MQ=255)
                                  aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                            >  1:1023308/1‑100 (MQ=255)
                                  aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                            <  1:51602/100‑1 (MQ=255)
                                     tCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCt                                                         >  1:132718/1‑100 (MQ=255)
                                     tCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCt                                                         >  2:3450280/1‑100 (MQ=255)
                                      cAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCtt                                                        <  2:653776/100‑1 (MQ=255)
                                        gCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCg                                                      >  2:242512/1‑100 (MQ=255)
                                          aaGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCt                                                    <  1:1397158/100‑1 (MQ=255)
                                           aGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCtt                                                   <  2:2189847/100‑1 (MQ=255)
                                           aGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCtt                                                   >  1:3356045/1‑100 (MQ=255)
                                            gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                  >  1:3453405/1‑100 (MQ=255)
                                            gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                  <  2:1203017/100‑1 (MQ=255)
                                               ccGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGc                                               >  2:1593434/1‑100 (MQ=255)
                                               ccGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGc                                               >  1:2742404/1‑100 (MQ=255)
                                                 gAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGt                                             >  2:48391/1‑100 (MQ=255)
                                                 gAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGt                                             >  2:121486/1‑100 (MQ=255)
                                                  aaaaCGGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTc                                            >  1:338323/1‑100 (MQ=255)
                                                   aaaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGGACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTcc                                           <  2:2437944/100‑1 (MQ=255)
                                                    aaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCa                                          <  1:2749680/100‑1 (MQ=255)
                                                      cAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGt                                        >  1:881085/1‑100 (MQ=255)
                                                      cAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGt                                        >  2:1495226/1‑100 (MQ=255)
                                                         gCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaa                                     <  2:546030/100‑1 (MQ=255)
                                                         gCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaa                                     >  2:503197/1‑100 (MQ=255)
                                                          cGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaaa                                    <  1:824101/100‑1 (MQ=255)
                                                          cGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaaa                                    >  1:2544959/1‑100 (MQ=255)
                                                           gTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAAt                                   <  1:2607920/100‑1 (MQ=255)
                                                           gTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAAt                                   >  1:3115038/1‑100 (MQ=255)
                                                              cAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTa                                 >  1:3012532/1‑99 (MQ=255)
                                                               aGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAAt                               <  2:889886/100‑1 (MQ=255)
                                                                gCTGGTCGCCTATGTGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                              >  1:2652744/1‑100 (MQ=255)
                                                                gCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                              >  1:75047/1‑100 (MQ=255)
                                                                gCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                              <  2:3013208/100‑1 (MQ=255)
                                                                  tGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGAt                            >  2:1150883/1‑100 (MQ=255)
                                                                   ggTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                           >  1:2363750/1‑100 (MQ=255)
                                                                   ggTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                           <  1:1410019/100‑1 (MQ=255)
                                                                   ggTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                           >  1:2362775/1‑100 (MQ=255)
                                                                    gTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAg                          <  1:98183/100‑1 (MQ=255)
                                                                     tCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGc                         <  2:69911/100‑1 (MQ=255)
                                                                     tCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGc                         <  1:2871931/100‑1 (MQ=255)
                                                                     tCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGc                         >  2:815176/1‑100 (MQ=255)
                                                                         cTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTg                     >  1:1342606/1‑100 (MQ=255)
                                                                          tATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTgg                    <  2:1355625/100‑1 (MQ=255)
                                                                            tGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGgag                  >  2:389997/1‑100 (MQ=255)
                                                                              cGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgat                >  1:887293/1‑100 (MQ=255)
                                                                              cGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgat                >  2:2761812/1‑100 (MQ=255)
                                                                               gTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatg               >  1:904849/1‑100 (MQ=255)
                                                                                 cGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatgat             >  2:2851631/1‑100 (MQ=255)
                                                                                  ggAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatgat             >  1:674748/1‑99 (MQ=255)
                                                                                  ggAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGAtt            >  2:1374504/1‑100 (MQ=255)
                                                                                   gaagaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTg           <  1:1374484/100‑1 (MQ=255)
                                                                                   gaagaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTg           >  2:3286838/1‑100 (MQ=255)
                                                                                   gaagaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTg           <  1:3456522/100‑1 (MQ=255)
                                                                                   gaagaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAACTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTg           <  1:2418034/100‑1 (MQ=255)
                                                                                    aagaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTgg          <  2:1091261/100‑1 (MQ=255)
                                                                                     agaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGc         <  2:1205421/100‑1 (MQ=255)
                                                                                     agaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCGTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGc         <  1:381697/100‑1 (MQ=255)
                                                                                      gaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCt        >  1:1243630/1‑100 (MQ=255)
                                                                                      gaaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCt        >  1:3041436/1‑100 (MQ=255)
                                                                                       aaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTa       >  1:1298703/1‑100 (MQ=255)
                                                                                         ccGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGTATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTAtt     <  1:743497/100‑1 (MQ=255)
                                                                                         ccGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTAtt     >  2:2465335/1‑100 (MQ=255)
                                                                                           gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc   <  1:474659/100‑1 (MQ=255)
                                                                                           gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc   >  1:2999319/1‑100 (MQ=255)
                                                                                           gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc   <  2:2688358/100‑1 (MQ=255)
                                                                                           gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc   <  2:1949567/100‑1 (MQ=255)
                                                                                            cccTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACg  <  1:3450226/100‑1 (MQ=255)
                                                                                                |                                                                                               
GTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACG  >  NC_000913/2386269‑2386460

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