Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,213,926 T→C T191A (ACG→GCG)  yehW ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,213,9260TC89.1% 661.1 / 58.6 230T191A (ACG→GCG) yehWputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base T (11/14);  new base C (99/106);  total (110/120)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 238 total aligned reads displayed.

CCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTT  >  NC_000913/2213836‑2214016
                                                                                          |                                                                                          
ccaccaGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTg                                                                                          >  1:2790461/1‑93 (MQ=255)
 caccaGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGa                                                                                         <  1:4281873/93‑1 (MQ=255)
  cccaGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAg                                                                                        <  2:5397220/92‑1 (MQ=255)
    caGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGc                                                                                      >  2:6000545/1‑93 (MQ=255)
    caGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGc                                                                                      <  2:1016905/93‑1 (MQ=255)
      gTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGa                                                                                    >  2:6639874/1‑93 (MQ=255)
      gTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGa                                                                                    >  2:3385341/1‑93 (MQ=255)
      gTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGa                                                                                    >  2:4126347/1‑93 (MQ=255)
      gTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGa                                                                                    >  1:481055/1‑93 (MQ=255)
       tAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAt                                                                                   <  2:4545238/93‑1 (MQ=255)
       tAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAt                                                                                   <  2:3264942/93‑1 (MQ=255)
       tAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAt                                                                                   <  1:4727163/93‑1 (MQ=255)
       tAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAt                                                                                   <  2:1534212/93‑1 (MQ=255)
       tAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAt                                                                                   <  2:2161362/93‑1 (MQ=255)
        aaCGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAtc                                                                                  <  1:2432954/93‑1 (MQ=255)
         aCGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAtcg                                                                                 <  2:4800956/93‑1 (MQ=255)
         aCGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAtcg                                                                                 <  2:2959675/93‑1 (MQ=255)
          cGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAtcgt                                                                                <  2:4275917/93‑1 (MQ=255)
          cGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGAtcgt                                                                                <  2:806410/93‑1 (MQ=255)
            cccccTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgtcg                                                                              >  2:1720176/1‑93 (MQ=255)
             ccccTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGCTGATGGGCGTACCCAGCGTGCTGGCCCCTCCCGCTGAGGCGATCGTCGc                                                                             <  2:4011335/93‑1 (MQ=255)
              cccTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGcc                                                                            >  2:3760338/1‑93 (MQ=255)
              cccTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGcc                                                                            <  2:474718/93‑1 (MQ=255)
                cTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGt                                                                          >  2:513383/1‑93 (MQ=255)
                    aTCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCa                                                                      >  1:3086775/1‑93 (MQ=255)
                    aTCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCa                                                                      >  1:2510646/1‑93 (MQ=255)
                        caTACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGt                                                                  <  2:2486855/93‑1 (MQ=255)
                        caTACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGt                                                                  >  1:1872750/1‑93 (MQ=255)
                         aTACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGtt                                                                 <  2:2675971/93‑1 (MQ=255)
                           aCGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGa                                                               <  1:6578806/93‑1 (MQ=255)
                             gcgGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATa                                                             <  1:4986115/93‑1 (MQ=255)
                             gcgGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATa                                                             <  2:3150144/93‑1 (MQ=255)
                              cgGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATaa                                                            >  1:6585773/1‑93 (MQ=255)
                                  aTTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCAc                                                        >  2:2857163/1‑93 (MQ=255)
                                   ttAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCAcc                                                       >  1:3802650/1‑93 (MQ=255)
                                   ttAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCAcc                                                       >  1:6266727/1‑93 (MQ=255)
                                   ttAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCAcc                                                       >  2:574154/1‑93 (MQ=255)
                                    tAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGACCGTCGCCGTACCAATGTTGATAATCACCg                                                      >  1:6528993/1‑93 (MQ=255)
                                     aaaTCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGa                                                     <  1:5033939/93‑1 (MQ=255)
                                     aaaTCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGa                                                     >  2:1913286/1‑93 (MQ=255)
                                      aaTCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGaa                                                    <  1:2651813/93‑1 (MQ=255)
                                      aaTCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGaa                                                    <  2:6055047/93‑1 (MQ=255)
                                       aTCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAg                                                   >  2:4408689/1‑93 (MQ=255)
                                        tCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGt                                                  >  1:4846870/1‑93 (MQ=255)
                                        tCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGt                                                  <  2:5692186/93‑1 (MQ=255)
                                        tCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGt                                                  >  1:1284781/1‑93 (MQ=255)
                                        tCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGt                                                  >  1:1139942/1‑93 (MQ=255)
                                        tCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGt                                                  <  1:1941131/93‑1 (MQ=255)
                                        tCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGt                                                  >  1:3877714/1‑93 (MQ=255)
                                         ccGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGtt                                                 >  2:6181129/1‑93 (MQ=255)
                                         ccGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGtt                                                 <  1:2210695/93‑1 (MQ=255)
                                          cGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTc                                                <  1:537808/93‑1 (MQ=255)
                                           gCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCg                                               <  2:2895869/93‑1 (MQ=255)
                                            cTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGc                                              <  2:3571405/93‑1 (MQ=255)
                                             tAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCa                                             >  1:3488593/1‑93 (MQ=255)
                                              aaGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAc                                            >  1:3256656/1‑93 (MQ=255)
                                                 cccGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAcgcc                                         <  1:5888994/93‑1 (MQ=255)
                                                 cccGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAcgcc                                         <  2:1379757/93‑1 (MQ=255)
                                                  ccGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAcgccc                                        <  2:6533164/93‑1 (MQ=255)
                                                  ccGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAcgccc                                        >  2:2811950/1‑93 (MQ=255)
                                                    gatgatgatGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAcgcccgc                                      >  2:3515032/1‑93 (MQ=255)
                                                     atgatgatgGGCGTACCCAGCGTGCTAGCCCCTACCGCTGAGGCGATCGGCGCCGTACCAATGTTGATAATCACCGAAGTTCGCAcgcccgcc                                     >  1:6029287/1‑93 (MQ=255)
                                                        gtgatgCGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGa                                  <  1:3157710/92‑1 (MQ=255)
                                                        atgatgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGa                                  <  2:4874696/93‑1 (MQ=255)
                                                        atgatgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGa                                  >  2:3334471/1‑93 (MQ=255)
                                                        atgatgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGa                                  <  1:3260061/93‑1 (MQ=255)
                                                          gatgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAAt                                <  1:4824731/93‑1 (MQ=255)
                                                          gatgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAAt                                >  1:3169707/1‑93 (MQ=255)
                                                           atgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATc                               <  1:608532/93‑1 (MQ=255)
                                                            tgGGCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCa                              >  1:1519214/1‑93 (MQ=255)
                                                             gggCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCAc                             <  1:1429196/93‑1 (MQ=255)
                                                             gggCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCAc                             >  2:3441731/1‑93 (MQ=255)
                                                              ggCGTACCCGGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCAcc                            <  1:157967/93‑1 (MQ=255)
                                                              ggCGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCAcc                            <  1:892898/93‑1 (MQ=255)
                                                                cGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCgg                          >  1:1185591/1‑93 (MQ=255)
                                                                cGTACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCgg                          <  1:2883042/93‑1 (MQ=255)
                                                                 ctACCCAGCGTGCTGCCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGa                         <  1:5068630/92‑1 (MQ=255)
                                                                  tACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAg                        >  2:836402/1‑93 (MQ=255)
                                                                  tACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAg                        >  2:508775/1‑93 (MQ=255)
                                                                  tACCCAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAg                        <  1:6394338/93‑1 (MQ=255)
                                                                    cccAGCGTGCGGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAgcc                      <  2:2365582/93‑1 (MQ=255)
                                                                     ccAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAgccg                     >  2:5352199/1‑93 (MQ=255)
                                                                      cAGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTGCCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAgccgc                    <  1:6348392/93‑1 (MQ=255)
                                                                       aGCGTGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAgccgcc                   <  2:5226798/93‑1 (MQ=255)
                                                                          gTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGc                <  1:4116381/93‑1 (MQ=255)
                                                                           tGCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCg               <  1:425584/93‑1 (MQ=255)
                                                                            gCTGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCgg              <  2:6336502/93‑1 (MQ=255)
                                                                            gCTGCCCCCTACCGCTGAGGCGATCGTCGCGGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCgg              <  1:4552935/93‑1 (MQ=255)
                                                                              tGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTa            >  1:6370266/1‑93 (MQ=255)
                                                                              tGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTa            <  1:5223897/93‑1 (MQ=255)
                                                                              tGGCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTa            >  2:6287654/1‑93 (MQ=255)
                                                                               ggCCCCTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAg           >  1:2343062/1‑93 (MQ=255)
                                                                                  cccTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTc        >  2:5774488/1‑93 (MQ=255)
                                                                                  cccTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTc        >  2:1360561/1‑93 (MQ=255)
                                                                                  cccTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTc        >  1:5702554/1‑93 (MQ=255)
                                                                                   ccTACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCg       >  2:5934735/1‑93 (MQ=255)
                                                                                     tACCGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGAc     <  2:2738615/93‑1 (MQ=255)
                                                                                       ccGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCt   >  2:2046355/1‑93 (MQ=255)
                                                                                       ccGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCt   >  2:3713089/1‑93 (MQ=255)
                                                                                       ccGCTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCt   <  1:1758710/93‑1 (MQ=255)
                                                                                          |                                                                                          
CCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCAGAATCACCGGAGCCGCCAGCGGTAGCTCGACCTT  >  NC_000913/2213836‑2214016

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