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.9% 535.8 / ‑5.0 180S17P (TCC→CCC) nhaBNa+/H+ antiporter, pH independent
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (91/87);  total (92/88)
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.88e-01

Only 100 of 194 total aligned reads displayed.

GCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTA  >  NC_000913/1233802‑1233981
                                                                                          |                                                                                         
gCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGGTTGTACCAGTCGGGGGGCt                                                                                         >  2:3177989/1‑93 (MQ=255)
 cAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTg                                                                                        <  2:7083559/93‑1 (MQ=255)
  aaCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTgg                                                                                       >  2:3033930/1‑93 (MQ=255)
   aCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGc                                                                                      >  1:7172374/1‑93 (MQ=255)
   aCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGc                                                                                      <  1:422337/93‑1 (MQ=255)
   aCCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGc                                                                                      >  2:4174593/1‑93 (MQ=255)
    ccAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGcc                                                                                     >  2:1254428/1‑93 (MQ=255)
     cAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGccc                                                                                    <  2:6752601/93‑1 (MQ=255)
      aGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCa                                                                                   >  1:7559283/1‑93 (MQ=255)
       gCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCaa                                                                                  >  2:7936791/1‑93 (MQ=255)
             gacgaAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAAtt                                                                            <  2:5672787/93‑1 (MQ=255)
             gacgaAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAAtt                                                                            <  2:4375067/93‑1 (MQ=255)
                 aaaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTgcgc                                                                        >  1:6704764/1‑93 (MQ=255)
                 aaaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTgcgc                                                                        <  2:4791223/93‑1 (MQ=255)
                  aaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGcc                                                                       >  1:3682488/1‑93 (MQ=255)
                  aaGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGcc                                                                       >  2:2101478/1‑93 (MQ=255)
                    gggCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCAt                                                                     >  1:7748341/1‑93 (MQ=255)
                    gggCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCAt                                                                     <  2:716495/93‑1 (MQ=255)
                     ggCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATa                                                                    <  1:5636355/93‑1 (MQ=255)
                     ggCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATa                                                                    <  1:6406335/93‑1 (MQ=255)
                     ggCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATa                                                                    <  1:4947657/93‑1 (MQ=255)
                      gCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAg                                                                   <  1:4405705/93‑1 (MQ=255)
                      gCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAg                                                                   >  1:4712234/1‑93 (MQ=255)
                       cTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgc                                                                  <  2:3143984/93‑1 (MQ=255)
                       cTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgc                                                                  <  2:767804/93‑1 (MQ=255)
                         gatgaGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgcgc                                                                >  1:625414/1‑93 (MQ=255)
                           tgaGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAgcgcgc                                                              <  1:3364138/93‑1 (MQ=255)
                             aGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCgg                                                            <  2:2793233/93‑1 (MQ=255)
                             aGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCgg                                                            >  1:2978426/1‑93 (MQ=255)
                              ggAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGc                                                           >  1:1022377/1‑93 (MQ=255)
                              ggAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGc                                                           >  1:3277784/1‑93 (MQ=255)
                               gAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCggtc                                                          >  2:7815277/1‑91 (MQ=255)
                                 aaaGTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGcccc                                                        <  2:1095339/93‑1 (MQ=255)
                                  aaTTAAGGGTTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCa                                                       <  2:7135814/93‑1 (MQ=255)
                                  aaTTAAGGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCAAAAAAATTGCGCCATAGCGCGCGGCCCCa                                                       <  2:1606538/93‑1 (MQ=255)
                                  aaTTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCa                                                       <  1:919276/93‑1 (MQ=255)
                                   ataaaCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAg                                                      <  2:2899149/90‑1 (MQ=255)
                                   aTTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAg                                                      >  2:2628731/1‑93 (MQ=255)
                                   aTTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAg                                                      >  1:5123122/1‑93 (MQ=255)
                                    ttAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAgg                                                     >  2:1387799/1‑93 (MQ=255)
                                     tAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGga                                                    <  2:3758581/93‑1 (MQ=255)
                                     tAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGga                                                    <  1:2647548/93‑1 (MQ=255)
                                      aaCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgag                                                   >  2:6241517/1‑93 (MQ=255)
                                      aaCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgag                                                   <  1:1096700/93‑1 (MQ=255)
                                      aaCGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGgag                                                   >  1:6580708/1‑93 (MQ=255)
                                         gggTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGAtc                                                >  1:7873434/1‑93 (MQ=255)
                                            tttACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTcc                                             >  1:225359/1‑93 (MQ=255)
                                             ttACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCa                                            <  1:6443282/93‑1 (MQ=255)
                                             ttACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCa                                            >  1:335898/1‑93 (MQ=255)
                                             ttACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCa                                            >  1:5717041/1‑93 (MQ=255)
                                               aCGATTAAGAAAATGATAAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatg                                          >  2:5432724/1‑93 (MQ=255)
                                                  aTTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatgatg                                       >  1:1890206/1‑93 (MQ=255)
                                                  aTTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCatgatg                                       >  2:322334/1‑93 (MQ=255)
                                                   ttAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATgg                                      >  1:1845179/1‑93 (MQ=255)
                                                   ttAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATgg                                      >  1:3794615/1‑93 (MQ=255)
                                                    tAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCGAGGAGATCTCCATGATGGt                                     >  1:7698204/1‑93 (MQ=255)
                                                    tAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGt                                     >  1:2899449/1‑93 (MQ=255)
                                                     aaGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGtt                                    <  1:3900829/93‑1 (MQ=255)
                                                     aaGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGtt                                    >  2:2038056/1‑93 (MQ=255)
                                                      aGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGttt                                   <  1:888918/93‑1 (MQ=255)
                                                        aaaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGTCTGGCCCAAAAAATTGCGTCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTcc                                 <  2:6457604/93‑1 (MQ=255)
                                                        aaaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTcc                                 >  1:3971269/1‑93 (MQ=255)
                                                         aaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGGCCCAGGAGATCTCCATGATGGTTTccc                                >  1:7420837/1‑93 (MQ=255)
                                                         aaaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTccc                                <  2:2286582/93‑1 (MQ=255)
                                                          aaTGATGAGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCcct                               >  2:2089389/1‑93 (MQ=255)
                                                            tgatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCcctta                             >  2:2379066/1‑93 (MQ=255)
                                                            tgatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCcctta                             >  1:7086223/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            >  2:6827902/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            >  2:5170703/1‑93 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            <  1:6096697/93‑1 (MQ=255)
                                                             gatgaGGGCGAGTTTGTACCAGTAGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttac                            <  2:4281325/93‑1 (MQ=255)
                                                              atgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacc                           <  2:7404920/93‑1 (MQ=255)
                                                               tgaGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacct                          >  2:4909787/1‑93 (MQ=255)
                                                                 aGGGCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCccttacctta                        >  2:3298041/1‑93 (MQ=255)
                                                                  gggCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCcttaccttac                       >  1:4142068/1‑93 (MQ=255)
                                                                  gggCGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCcttaccttac                       <  1:2797878/93‑1 (MQ=255)
                                                                  gggCGAGTGTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCcttgccttac                       <  2:7513538/93‑1 (MQ=255)
                                                                     cGAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACaaa                    <  1:6798359/93‑1 (MQ=255)
                                                                      gAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaat                   >  2:3460213/1‑93 (MQ=255)
                                                                      gAGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaat                   <  1:5608090/93‑1 (MQ=255)
                                                                       aGTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaata                  <  2:3957152/93‑1 (MQ=255)
                                                                        gTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataa                 >  1:2889053/1‑93 (MQ=255)
                                                                        gTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataa                 <  1:2150036/93‑1 (MQ=255)
                                                                        gTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataa                 >  2:5514762/1‑93 (MQ=255)
                                                                        gTTTGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataa                 >  2:5714181/1‑93 (MQ=255)
                                                                         tttGTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataat                >  2:6176987/1‑93 (MQ=255)
                                                                            gTACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCaa             <  1:2238816/93‑1 (MQ=255)
                                                                             tACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAat            >  1:4817775/1‑93 (MQ=255)
                                                                             tACCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAat            >  2:2710571/1‑93 (MQ=255)
                                                                              aCCAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatg           <  1:3547883/93‑1 (MQ=255)
                                                                               ccAGTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAatga          <  2:403357/93‑1 (MQ=255)
                                                                                 aGTCGGGGGGCTGCACCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTGCAAATAATCAatgatg        <  2:7057021/93‑1 (MQ=255)
                                                                                  gTCGGGGGGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGt       <  1:2220114/93‑1 (MQ=255)
                                                                                   cgggggggCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGtt      >  1:2453299/3‑93 (MQ=255)
                                                                                      ggggGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGttttt   >  2:2723402/1‑93 (MQ=255)
                                                                                      ggggGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGttttt   <  1:5231264/93‑1 (MQ=255)
                                                                                      ggggGCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGttttt   >  1:5713365/1‑93 (MQ=255)
                                                                                       ggggCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTa  <  2:4629355/93‑1 (MQ=255)
                                                                                       ggggCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTa  <  1:6401136/93‑1 (MQ=255)
                                                                                       ggggCTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTa  <  1:1466100/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTA  >  NC_000913/1233802‑1233981

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