Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,233,892 A→G S17P (TCC→CCC)  nhaB ← Na+/H+ antiporter, pH independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,8920AG98.5% 397.5 / ‑2.9 133S17P (TCC→CCC) nhaBNa+/H+ antiporter, pH independent
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (84/47);  total (85/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

Only 100 of 144 total aligned reads displayed.

GCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTA  >  NC_000913/1233802‑1233981
                                                                                          |                                                                                         
gCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCt                                                                                         <  2:3609572/93‑1 (MQ=255)
 cAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTg                                                                                        >  2:302382/1‑93 (MQ=255)
  aaCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTgg                                                                                       >  1:1573992/1‑93 (MQ=255)
   aCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGc                                                                                      >  2:3604184/1‑93 (MQ=255)
   aCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGc                                                                                      >  1:392823/1‑93 (MQ=255)
   aCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGc                                                                                      >  1:2686733/1‑93 (MQ=255)
    ccAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGcc                                                                                     >  2:4230793/1‑93 (MQ=255)
    ccAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGGTTGTACCAGTCGGGGGGCTGGcc                                                                                     >  2:4230794/1‑93 (MQ=255)
    ccAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGAGGAGGGCGAGCTTGTACCAGTCGGGGGGCTGGcc                                                                                     >  1:590573/1‑93 (MQ=255)
     cAGCCCGCGACGAAAGGGCTGATGAGTAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGccc                                                                                    >  1:3088461/1‑93 (MQ=255)
     cAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGccc                                                                                    >  1:2207580/1‑93 (MQ=255)
      aGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCa                                                                                   >  1:1762255/1‑93 (MQ=255)
      aGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCa                                                                                   >  1:3016140/1‑93 (MQ=255)
      aGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCa                                                                                   >  2:731634/1‑93 (MQ=255)
       gCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCaa                                                                                  >  1:370262/1‑93 (MQ=255)
       gCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCaa                                                                                  >  1:747814/1‑93 (MQ=255)
          cgcgACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCaaaaa                                                                               >  2:3857525/1‑93 (MQ=255)
           gcgACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCTAGTTTGTACCAGTCGGGGGGCTGGCCCaaaaaa                                                                              >  1:4581181/1‑93 (MQ=255)
           gcgACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCaaaaaa                                                                              >  1:2064587/1‑93 (MQ=255)
             gacgaAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAAtt                                                                            >  2:2546843/1‑93 (MQ=255)
              acgaAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTg                                                                           <  1:4472323/93‑1 (MQ=255)
               cgaAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGTGGCTGGCCCAAAAAATTgc                                                                          >  1:2374478/1‑93 (MQ=255)
               cgaAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAAttac                                                                          >  1:548051/1‑91 (MQ=255)
                 aaaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTgcgc                                                                        >  1:2177251/1‑93 (MQ=255)
                  aaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGcc                                                                       >  2:4384455/1‑93 (MQ=255)
                  aaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGcc                                                                       >  2:3737509/1‑93 (MQ=255)
                   aGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCa                                                                      >  2:3337056/1‑93 (MQ=255)
                    gggCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCAt                                                                     <  2:2057912/93‑1 (MQ=255)
                      gCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAg                                                                   <  2:4671596/93‑1 (MQ=255)
                       cTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgc                                                                  >  2:2312721/1‑93 (MQ=255)
                         gatgaGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgcgc                                                                >  2:2231548/1‑93 (MQ=255)
                           tgaGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgcgcgc                                                              >  1:2146783/1‑93 (MQ=255)
                             aGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCgg                                                            <  2:2040489/93‑1 (MQ=255)
                             aGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCgg                                                            >  1:2393597/1‑93 (MQ=255)
                              ggAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGc                                                           >  2:758678/1‑93 (MQ=255)
                               gAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGcc                                                          <  1:1307327/93‑1 (MQ=255)
                                 aaaTTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGcccc                                                        >  1:4341738/1‑93 (MQ=255)
                                   aTTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGACCCAg                                                      >  2:1388997/1‑93 (MQ=255)
                                      aaCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgag                                                   >  1:4251609/1‑93 (MQ=255)
                                       aCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgaga                                                  <  1:3231223/93‑1 (MQ=255)
                                       aCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgaga                                                  <  1:3231224/93‑1 (MQ=255)
                                       aCGGGTTAACGATTAAGAAAAGGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgaga                                                  <  1:939443/93‑1 (MQ=255)
                                        cGGGTTTGCGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCAGGCACAAAAAATTGCGCAATAGCGAGCGGCCCCAGAAGAt                                                 >  2:4653114/1‑93 (MQ=255)
                                        cGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGAt                                                 <  1:4403189/93‑1 (MQ=255)
                                         gggTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGAtc                                                >  1:3115521/1‑93 (MQ=255)
                                          ggTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCCTAGCGCGCGGCCCCAGGAGAtct                                               >  2:2587510/1‑93 (MQ=255)
                                          ggTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGAtct                                               >  1:3629594/1‑93 (MQ=255)
                                            tttACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTcc                                             >  2:1665640/1‑93 (MQ=255)
                                            tttACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTcc                                             <  1:3791238/93‑1 (MQ=255)
                                             ttACGTTTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCAAAAAAATTGCGCCATAGCGCGCCGCCCCAGGAGATCTCCa                                            >  2:2989885/1‑93 (MQ=255)
                                              tACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCat                                           >  1:944754/1‑93 (MQ=255)
                                                cGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatga                                         <  1:3375886/93‑1 (MQ=255)
                                                cGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatga                                         >  1:4238816/1‑93 (MQ=255)
                                                cGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatga                                         <  2:3123057/93‑1 (MQ=255)
                                                cGACTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatga                                         <  1:1221342/93‑1 (MQ=255)
                                                 gATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatgat                                        <  1:1785923/93‑1 (MQ=255)
                                                  aTTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatgatg                                       >  1:1773247/1‑93 (MQ=255)
                                                   ttAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATgg                                      >  2:4109728/1‑93 (MQ=255)
                                                    tAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGt                                     <  1:4683423/93‑1 (MQ=255)
                                                    tAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGt                                     >  1:4288381/1‑93 (MQ=255)
                                                     aaGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGtt                                    >  2:1424169/1‑93 (MQ=255)
                                                     aaGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGtt                                    <  2:2536219/93‑1 (MQ=255)
                                                       gAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTc                                  >  1:2275133/1‑93 (MQ=255)
                                                       gAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTc                                  >  1:2275132/1‑93 (MQ=255)
                                                        aaaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTcc                                 >  2:729735/1‑93 (MQ=255)
                                                        aaaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTcc                                 >  1:4040813/1‑93 (MQ=255)
                                                        aaaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTcc                                 <  1:2419256/93‑1 (MQ=255)
                                                          aaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATGGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCcct                               >  2:3935113/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            >  2:468973/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            >  2:1645405/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            >  1:2786929/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            >  1:651295/1‑93 (MQ=255)
                                                               tgaGTGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacct                          >  1:2574978/1‑93 (MQ=255)
                                                               tgaGGGCGAGTGTGTACCAGTCGGGGCGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacct                          <  2:3510091/93‑1 (MQ=255)
                                                               ggggggCGAGTTTGTCCCAGTGGGGGGTCTGGCCAAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacct                          <  2:1712240/90‑1 (MQ=255)
                                                                gaGGGCGAGTTTGTCCCAGTCGGGGGGCGGCCCAAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacctt                         <  2:1731871/93‑1 (MQ=255)
                                                                gaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacctt                         >  1:3545535/1‑93 (MQ=255)
                                                                  gggCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCcttaccttac                       <  2:3953578/93‑1 (MQ=255)
                                                                    gCGAGTTTGTCCCAGGCGGGGGGCTGGCCAAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACaa                     <  2:894499/93‑1 (MQ=255)
                                                                       aGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaata                  >  2:2647775/1‑93 (MQ=255)
                                                                       aGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaata                  <  1:4484995/93‑1 (MQ=255)
                                                                        gTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataa                 >  2:4414624/1‑93 (MQ=255)
                                                                         tttGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataat                <  1:1664326/93‑1 (MQ=255)
                                                                          ttGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATc               <  2:526139/93‑1 (MQ=255)
                                                                          ttGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATc               >  1:2531826/1‑93 (MQ=255)
                                                                            gTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCaa             >  2:4569840/1‑93 (MQ=255)
                                                                             tACCAGTCGGGGTGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAat            <  1:3051997/93‑1 (MQ=255)
                                                                              aCCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatg           >  2:3470388/1‑93 (MQ=255)
                                                                              aCCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatg           >  1:1637437/1‑93 (MQ=255)
                                                                               ccAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatga          <  1:4178554/93‑1 (MQ=255)
                                                                               ccAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatga          <  1:109618/93‑1 (MQ=255)
                                                                               ccAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatga          <  2:3479772/93‑1 (MQ=255)
                                                                                cAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatgat         <  2:1625987/93‑1 (MQ=255)
                                                                                 aGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatgatg        >  1:4343975/1‑93 (MQ=255)
                                                                                  gTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGt       >  2:3533380/1‑93 (MQ=255)
                                                                                  gTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGt       >  1:80072/1‑93 (MQ=255)
                                                                                  gTCGGGGGGCTCGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATGAATGATGt       >  2:2784975/1‑93 (MQ=255)
                                                                                   tCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGtt      <  1:3125181/93‑1 (MQ=255)
                                                                                    cGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGttt     >  2:4176916/1‑93 (MQ=255)
                                                                                       ggggCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTa  >  1:1115735/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTA  >  NC_000913/1233802‑1233981

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