Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,175,196 C→T 8.3% L91F (CTC→TTC)  PP_1030 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,175,1960CT 8.3% 166.2 / 4.3 60L91F (CTC→TTC) PP_1030hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (2/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTG  >  NC_002947/1175048‑1175342
                                                                                                                                                     |                                                                                                                                                  
tgGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTcc                                                                                                                                                    <  7:61573/149‑1 (MQ=255)
tgGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTcc                                                                                                                                                    <  7:350838/149‑1 (MQ=255)
 ggCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCt                                                                                                                                                   >  8:258126/1‑149 (MQ=255)
           gTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGcc                                                                                                                                         >  3:44714/1‑149 (MQ=255)
            tAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCGAGGAAGcct                                                                                                                                        <  5:78966/149‑2 (MQ=255)
              aCCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCTTCGAGGAAGccttc                                                                                                                                      <  6:209129/149‑4 (MQ=255)
               ccGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCg                                                                                                                                        >  5:67335/1‑146 (MQ=255)
               ccGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCg                                                                                                                                        <  6:67335/146‑1 (MQ=255)
                cGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTg                                                                                                                                    >  1:160329/1‑149 (MQ=255)
                   aTGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGccc                                                                                                                                 >  6:370035/1‑149 (MQ=255)
                                  gcggcAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGACCCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGa                                                                                                                  <  5:382401/149‑1 (MQ=255)
                                  gcggcAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGa                                                                                                                  <  6:13586/149‑1 (MQ=255)
                                  gcggcAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGa                                                                                                                  <  7:258126/149‑1 (MQ=255)
                                    ggcAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGagg                                                                                                                <  2:279920/149‑1 (MQ=255)
                                       aGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGaggtag                                                                                                             <  5:9872/149‑1 (MQ=255)
                                       aGGGGCGCGCAAAGTGGCCCTGCACACGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACACCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGaggtag                                                                                                             >  2:45826/1‑149 (MQ=255)
                                           gcgcgcAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAACCGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGa                                                                                                         <  5:370035/149‑1 (MQ=255)
                                             gcgcAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGc                                                                                                       >  5:347653/1‑149 (MQ=255)
                                                 aaaGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCgtt                                                                                                   >  7:220243/1‑149 (MQ=255)
                                                        cccTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGGCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGaa                                                                                            <  4:19199/149‑1 (MQ=255)
                                                        cccTGCACCAGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGaa                                                                                            <  4:14128/149‑1 (MQ=255)
                                                         ccTGCACCCGGAAAAACTGCCCTTCGGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGCTCCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCg                                                                                              >  1:201053/1‑147 (MQ=255)
                                                         ccTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGCTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCg                                                                                              <  2:201053/146‑1 (MQ=255)
                                                                    aaaaaCTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCACACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAagag                                                                                >  1:297765/1‑149 (MQ=255)
                                                                     aaaaCTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCTGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGt                                                                               <  1:275252/149‑1 (MQ=255)
                                                                      aaaCTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTg                                                                              <  2:314696/149‑1 (MQ=255)
                                                                                    gCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTg                                                                 <  2:160329/149‑1 (MQ=255)
                                                                                     cGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGTAGt                                                                  >  7:145945/1‑147 (MQ=255)
                                                                                      gAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAg                                                               >  5:147885/1‑149 (MQ=255)
                                                                                           ggTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGAc                                                           >  5:122046/1‑148 (MQ=255)
                                                                                            gTCAACGGCCTGGACGTGATGATCAACCGCTGCATCTATCCCCCCACCGCAGGCTTCCTCGAAGCAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAACCCTGGAAGAGTGGATGCCGGCAGTGAGCGACaa                                                         <  6:303170/149‑1 (MQ=255)
                                                                                            gTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTg                                                                                      >  3:303389/1‑120 (MQ=255)
                                                                                            gTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTg                                                                                      <  4:303389/120‑1 (MQ=255)
                                                                                            gTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGGTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGAcga                                                         >  4:342053/1‑147 (MQ=255)
                                                                                             tCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAAc                                                        <  6:347653/149‑1 (MQ=255)
                                                                                              cAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTAAACCCCCACCGAAGGCTTCCTCGAAGAAGCCAGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACt                                                       >  1:248463/1‑149 (MQ=255)
                                                                                               aaCGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACtt                                                      >  6:62294/1‑149 (MQ=255)
                                                                                                 cGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTa                                                    >  2:341971/1‑149 (MQ=255)
                                                                                                   gCCTGTAAGTGATGCTCAAACGCTGCATCTATACCCACACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGAATGCCGGCAGTGAGCGACAACTTTAcc                                                  <  4:100797/149‑1 (MQ=255)
                                                                                                      tGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGc                                               >  5:202384/1‑149 (MQ=255)
                                                                                                                 cTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCgg                                    >  3:307896/1‑149 (MQ=255)
                                                                                                                 cTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCgg                                    >  8:178/1‑149 (MQ=255)
                                                                                                                     aaCGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTc                                <  7:178/149‑1 (MQ=255)
                                                                                                                        gCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGaa                             <  2:297765/149‑1 (MQ=255)
                                                                                                                         ctacATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCg                               >  3:104177/4‑146 (MQ=255)
                                                                                                                            cATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACCACTTTACCTGCCCACTGTGCGGCTTCGAAGAt                          <  1:45826/148‑1 (MQ=255)
                                                                                                                             ctccataCCCCCACCGCAGGCTTCCTCGCAGACGCCGGTTGCCCGCAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGa                        <  4:223250/145‑1 (MQ=255)
                                                                                                                                 ataCCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGc                                      >  4:301355/1‑131 (MQ=255)
                                                                                                                                 ataCCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGc                                      <  3:301355/131‑1 (MQ=255)
                                                                                                                                        cACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCt             >  4:300475/1‑149 (MQ=255)
                                                                                                                                          ccGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTc           <  4:338899/149‑1 (MQ=255)
                                                                                                                                            gAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCAt         <  3:300475/149‑1 (MQ=255)
                                                                                                                                            gAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCAt         >  6:337483/1‑149 (MQ=255)
                                                                                                                                            accggCTTCTTCGAGGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTAc                                                   <  3:115662/104‑1 (MQ=255)
                                                                                                                                            accggCTTCTTCGAGGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTAc                                                   >  4:115662/4‑107 (MQ=255)
                                                                                                                                            accggCTTCTTCGAGGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCAt         >  3:47656/4‑149 (MQ=255)
                                                                                                                                                gCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTAc     <  6:147885/149‑1 (MQ=255)
                                                                                                                                                gCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTAc     <  4:44714/149‑1 (MQ=255)
                                                                                                                                                gCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTAc     <  4:104177/149‑1 (MQ=255)
                                                                                                                                                 cTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTAcc    >  8:164515/1‑149 (MQ=255)
                                                                                                                                                  ttCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTAc                                                   <  2:229190/101‑1 (MQ=255)
                                                                                                                                                  ttCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGACGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTAc                                                   >  1:229190/1‑101 (MQ=255)
                                                                                                                                                    ccTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTg  >  5:360326/1‑148 (MQ=255)
                                                                                                                                                     cTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTg  >  2:314802/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTC‑GGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGATGACATCAATGGCTTCATTTACCTG  >  NC_002947/1175048‑1175342

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