Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,035,840 C→T G88G (GGC→GGT hyaE → HyaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,035,8400CT99.2% 426.0 / NA 127G88G (GGC→GGThyaEHyaE
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (67/58);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

TGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAA  >  NC_000913/1035750‑1035928
                                                                                          |                                                                                        
tGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGgtgt                                                                                        <  2:464121/93‑1 (MQ=255)
 ggCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTc                                                                                       <  2:2589202/93‑1 (MQ=255)
 ggCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTc                                                                                       <  2:3554413/93‑1 (MQ=255)
   cGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCtt                                                                                     >  2:3142591/1‑93 (MQ=255)
    gAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCttt                                                                                    <  1:116680/93‑1 (MQ=255)
     aaTTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTc                                                                                   >  2:4148129/1‑93 (MQ=255)
     aaTTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTc                                                                                   >  2:2597971/1‑93 (MQ=255)
       ttACTGCTCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGc                                                                                 <  1:2353718/93‑1 (MQ=255)
        tACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTTGCt                                                                                >  2:2983516/1‑93 (MQ=255)
        tACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCt                                                                                <  1:2143677/93‑1 (MQ=255)
        tACTGCGCGAGTTTCCCGACTATACATGGCAGGCGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCt                                                                                >  2:1774018/1‑93 (MQ=255)
         aCTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCtt                                                                               >  2:2179805/1‑93 (MQ=255)
         aCTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCtt                                                                               >  1:28249/1‑93 (MQ=255)
            gcgcGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTcc                                                                            <  2:85795/93‑1 (MQ=255)
            gcgcGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTcc                                                                            <  1:3707799/93‑1 (MQ=255)
            gcgcGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTcc                                                                            <  1:4461977/93‑1 (MQ=255)
             cgcgAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCt                                                                           <  2:3627722/93‑1 (MQ=255)
             cgcgAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCt                                                                           <  1:2767517/93‑1 (MQ=255)
              gcgAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTg                                                                          >  1:4636774/1‑93 (MQ=255)
              gcgAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTg                                                                          >  1:1093195/1‑93 (MQ=255)
                 aGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCa                                                                       >  2:1156435/1‑93 (MQ=255)
                  gTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCAc                                                                      <  2:3853205/93‑1 (MQ=255)
                  gTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCAc                                                                      >  1:706437/1‑93 (MQ=255)
                  gTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCAc                                                                      <  1:3424207/93‑1 (MQ=255)
                  gTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCAc                                                                      >  1:944939/1‑93 (MQ=255)
                  gTTTCCCGACTATACAGGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCAc                                                                      <  1:3891801/93‑1 (MQ=255)
                    tccccGACTATACATGGCAGGAGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACtt                                                                    <  2:3995172/91‑1 (MQ=255)
                      cccGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTa                                                                  >  2:735589/1‑93 (MQ=255)
                       ccGACTATGCAGGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAg                                                                 <  1:3249958/93‑1 (MQ=255)
                       ccGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAg                                                                 >  1:3789998/1‑93 (MQ=255)
                       ccGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAg                                                                 >  1:4365811/1‑93 (MQ=255)
                       ccGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAg                                                                 >  2:3284921/1‑93 (MQ=255)
                        cGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAgt                                                                <  2:3110541/93‑1 (MQ=255)
                         gACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAgtg                                                               >  1:4298941/1‑93 (MQ=255)
                          aCTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAgtgt                                                              >  1:1341888/1‑93 (MQ=255)
                                 aTGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACcgg                                                       <  1:3418114/93‑1 (MQ=255)
                                  tGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACcggc                                                      >  1:791965/1‑93 (MQ=255)
                                  tGGCAGGTGGCGAATGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACcggc                                                      >  2:2309752/1‑93 (MQ=255)
                                      aGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaa                                                  >  1:2632964/1‑93 (MQ=255)
                                      aGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaa                                                  >  1:3087588/1‑93 (MQ=255)
                                      aGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGGGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaa                                                  <  2:1974449/93‑1 (MQ=255)
                                       ggTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGTCGGaaa                                                 >  2:912869/1‑93 (MQ=255)
                                       ggTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaaa                                                 >  2:384212/1‑93 (MQ=255)
                                       ggTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaaa                                                 <  1:4605757/93‑1 (MQ=255)
                                       ggTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaaa                                                 <  1:2987337/93‑1 (MQ=255)
                                       ggTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGaaa                                                 >  1:3494724/1‑93 (MQ=255)
                                        gTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAAc                                                >  1:446163/1‑93 (MQ=255)
                                          ggCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTa                                              <  2:3318875/93‑1 (MQ=255)
                                          ggCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTa                                              >  1:2174021/1‑93 (MQ=255)
                                          ggCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTa                                              <  2:2072212/93‑1 (MQ=255)
                                          ggCGATTGCTGACCTGGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTa                                              <  1:1443252/93‑1 (MQ=255)
                                           gCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTAt                                             <  2:4140621/93‑1 (MQ=255)
                                            ggATTGCTGACCTGGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATc                                            <  1:3220896/92‑1 (MQ=255)
                                            cGATTGCTGCCCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATc                                            <  2:234008/93‑1 (MQ=255)
                                            cGATGGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATc                                            <  2:2222659/93‑1 (MQ=255)
                                              atggCTGACCTGGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATcgc                                          <  2:4236872/90‑1 (MQ=255)
                                              aTTGCTGACCTTGAGCAGAGGGAATCAATCGTCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATcgc                                          <  2:1704443/93‑1 (MQ=255)
                                              aTTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATctc                                          >  2:2550679/1‑91 (MQ=255)
                                                tGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGCAACTATCgcgg                                        >  1:805740/1‑93 (MQ=255)
                                                tGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCgcgg                                        <  2:1269991/93‑1 (MQ=255)
                                                   tGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGt                                     >  2:3162972/1‑93 (MQ=255)
                                                   tGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGt                                     <  2:3017952/93‑1 (MQ=255)
                                                       cTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTg                                 <  1:3425987/93‑1 (MQ=255)
                                                        ttGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGa                                >  2:3817297/1‑93 (MQ=255)
                                                        ttGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGa                                >  1:4641517/1‑93 (MQ=255)
                                                         taaGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTTTTTCGCTTCCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGaa                               <  1:241293/91‑1 (MQ=255)
                                                          gAGCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAAt                              <  2:1983870/93‑1 (MQ=255)
                                                            gCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATgg                            >  1:747997/1‑93 (MQ=255)
                                                            gCAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATgg                            >  1:1726690/1‑93 (MQ=255)
                                                             cAGAGCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGt                           <  1:4089914/93‑1 (MQ=255)
                                                              agagCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTa                          <  2:445200/93‑1 (MQ=255)
                                                              agagCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGGTTACCGGCGGAAACTATCGCGGCGTGGTGAATggca                          >  2:1805038/1‑91 (MQ=255)
                                                               gagCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTAt                         <  2:3663400/93‑1 (MQ=255)
                                                               gagCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTAt                         <  1:1456158/93‑1 (MQ=255)
                                                               gagCGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATAGTAt                         >  1:2156351/1‑93 (MQ=255)
                                                                  cGAAGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCa                      <  1:3641976/93‑1 (MQ=255)
                                                                    aaGCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCAcc                    >  2:2197183/1‑93 (MQ=255)
                                                                    aaGCCATCGGCGATCGCTTTGGTGTCTTGCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCAcc                    <  1:2062522/93‑1 (MQ=255)
                                                                      gCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCg                  >  1:2259942/1‑93 (MQ=255)
                                                                      gCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCg                  >  1:2003518/1‑93 (MQ=255)
                                                                      gCCATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCg                  >  2:780542/1‑93 (MQ=255)
                                                                       ccATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGt                 <  2:3341462/93‑1 (MQ=255)
                                                                        cATCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTg                >  1:3850804/1‑93 (MQ=255)
                                                                         aTCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTgg               >  2:1956307/1‑93 (MQ=255)
                                                                         aTCGGCGATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTgg               <  2:228715/93‑1 (MQ=255)
                                                                         aTCGGCGATCGCTTTGGTGCCTTTCGCTTTCCTGCCACTTTAGTGTTTACCAGCGGAAACTATCGCGGCGTGGTGAATGGTATTCACCCGTgg               >  1:1733074/1‑93 (MQ=255)
                                                                           cGGCGATCGTTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGc             <  1:2665712/93‑1 (MQ=255)
                                                                               gATCGTTTTGGTGTTTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGaa         <  1:909833/93‑1 (MQ=255)
                                                                               gATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGaa         >  1:3977842/1‑93 (MQ=255)
                                                                               gATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGaa         >  2:2072221/1‑93 (MQ=255)
                                                                               gATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGaa         >  1:3844409/1‑93 (MQ=255)
                                                                               gATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGaa         >  1:4396380/1‑93 (MQ=255)
                                                                               gATCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGaa         >  1:4495886/1‑93 (MQ=255)
                                                                                aTCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAAc        >  2:495382/1‑93 (MQ=255)
                                                                                aTCGCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAAc        >  1:3704401/1‑93 (MQ=255)
                                                                                 tCGCTGTGGTGTCTTTCGCGTCCCTGCCACTTTAGTGTTGACCGGCGGAACCTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACt       <  1:366573/93‑1 (MQ=255)
                                                                                   gCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGa     <  1:1208135/93‑1 (MQ=255)
                                                                                   gCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGa     >  2:889801/1‑93 (MQ=255)
                                                                                   gCTTTGGTGTCTTTCGCTTTCCTGCCACTTTAGGGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAGCTGa     >  2:4074737/1‑93 (MQ=255)
                                                                                      ttGGTGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATaa  >  1:1761139/1‑93 (MQ=255)
                                                                                          |                                                                                        
TGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGCAGAGCGAAGCCATCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAA  >  NC_000913/1035750‑1035928

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