Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,790,276 C→T 7.4% Q1188Q (CAG→CAA pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,790,2760CT 7.4% 200.4 / 4.6 67Q1188Q (CAG→CAApvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (2/3);  total (31/36)
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.62e-01

ATCACGCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCAATGATTCATGCCGTGACACCAGCGCCGAAAA  >  NC_002947/4790130‑4790422
                                                                                                                                                    |                                                                                                                                                   
aTCACGCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAActgctg                                                                                                                                                   <  5:147105/149‑1 (MQ=255)
  cACGCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACtgctgc                                                                                                                                                  >  6:128886/1‑148 (MQ=255)
     gCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTCAGGCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCaa                                                                                                                                                >  1:96405/1‑149 (MQ=255)
     gCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCaa                                                                                                                                                >  4:142035/1‑147 (MQ=255)
     gCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACa                                                                                                                                              <  2:181297/149‑1 (MQ=255)
            gACCACCCGTCGCACCCGTTGTGGGCCAGGGTCCCGAGCAGTAGGTGCCCCTGCCC‑‑GCCACGCTGAAACAGCTTCATCCCAATCAGCGGCCCGTTGGCAAGGGCTAAGGGCTGCTGGATGGCGGCCTGCAACTGTTGCAGCAGTTgcag                                                                                                                                       <  6:197176/149‑3 (MQ=255)
              ccAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCgg                                                                                                                                     <  5:128886/149‑1 (MQ=255)
               cAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTATGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCggg                                                                                                                                    <  6:160128/149‑1 (MQ=255)
               cAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCggg                                                                                                                                    <  7:78147/149‑1 (MQ=255)
               cAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCggg                                                                                                                                    <  4:60748/149‑1 (MQ=255)
                    gTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCa                                                                                                                                >  5:133364/1‑148 (MQ=255)
                    gTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCa                                                                                                                                >  4:256126/1‑148 (MQ=255)
                              aTGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTgcgc                                                                                                                    <  2:96405/149‑1 (MQ=255)
                               tgtgGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGcc                                                                                                                   <  3:256126/149‑1 (MQ=255)
                                   gtgCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCggg                                                                                                                                    <  6:276776/129‑1 (MQ=255)
                                   gtgCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCggg                                                                                                                                    >  5:276776/1‑129 (MQ=255)
                                   gtgCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCctgc                                                                                                                <  1:289780/148‑1 (MQ=255)
                                              cGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAgg                                                                                                    <  3:142035/149‑1 (MQ=255)
                                                aGCAGTAGGTGCTCCTGCTC‑‑TCCACGCTGATACAGCTGCATGCGAATCAGCGGCACGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGt                                                                                                  <  8:249220/149‑1 (MQ=255)
                                                    gTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCtgctgctgc                                                                                                          <  5:111186/137‑1 (MQ=255)
                                                    gTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCtgctgctgc                                                                                                          >  6:111186/1‑137 (MQ=255)
                                                         tGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACcgc                                                                                            >  4:280922/1‑146 (MQ=255)
                                                            tCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGTTGCAGCAGTTgc                                                                                                                                         >  5:156903/1‑99 (MQ=255)
                                                            tCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGTTGCAGCAGTTgc                                                                                                                                         <  6:156903/99‑1 (MQ=255)
                                                            tCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGAGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTTACCGCGAGAtt                                                                                      >  2:156361/1‑149 (MQ=255)
                                                                 cTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGc                                                                                   <  4:173222/147‑1 (MQ=255)
                                                                      gCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGTAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTg                                                                                <  2:115012/147‑1 (MQ=255)
                                                                      gCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGa                                                                              >  8:29012/1‑149 (MQ=255)
                                                                        cACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGc                                                                                                                           <  7:152272/103‑1 (MQ=255)
                                                                        cACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGc                                                                                                                           >  8:152272/1‑103 (MQ=255)
                                                                        cACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGAcg                                                                            >  5:294671/1‑149 (MQ=255)
                                                                         aCGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGAcgc                                                                           <  1:156361/149‑1 (MQ=255)
                                                                         aCGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGAcgc                                                                           <  2:196299/149‑1 (MQ=255)
                                                                          cGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGAcgcg                                                                          <  3:280922/149‑1 (MQ=255)
                                                                                 aCAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGc                                                                     >  3:45455/1‑147 (MQ=255)
                                                                                   aGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAgg                                                                 <  4:45455/149‑1 (MQ=255)
                                                                                         aTGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCaa                                                                       <  5:204109/137‑1 (MQ=255)
                                                                                         aTGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCaa                                                                       >  6:204109/1‑137 (MQ=255)
                                                                                           gCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCt                                                          <  5:125719/148‑1 (MQ=255)
                                                                                           gCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCt                                                          >  6:125719/1‑148 (MQ=255)
                                                                                            cGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTtg                                                        <  3:225763/149‑1 (MQ=255)
                                                                                             gAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTtgc                                                       >  5:143995/1‑149 (MQ=255)
                                                                                               aTCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTtgctg                                                     <  6:294671/149‑1 (MQ=255)
                                                                                                 cAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTgctgct                                                   >  7:8422/1‑149 (MQ=255)
                                                                                                    cGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCtt                                                >  4:91803/1‑149 (MQ=255)
                                                                                                     ggCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTc                                               <  8:8422/149‑1 (MQ=255)
                                                                                                         cGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGca                                                               >  4:89142/1‑129 (MQ=255)
                                                                                                         cGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGca                                                               <  3:89142/129‑1 (MQ=255)
                                                                                                           ttGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACa                                         <  3:91803/149‑1 (MQ=255)
                                                                                                                aaGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTAc                                    >  6:31822/1‑149 (MQ=255)
                                                                                                                 aGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACg                                   <  7:29012/149‑1 (MQ=255)
                                                                                                                    tCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTtgct                                                      <  5:4999/127‑1 (MQ=255)
                                                                                                                    tCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTtgct                                                      >  6:4999/1‑127 (MQ=255)
                                                                                                                    tCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCaa                                >  3:100754/1‑149 (MQ=255)
                                                                                                                      gAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACa                                         >  4:57174/1‑138 (MQ=255)
                                                                                                                      gAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACa                                         >  4:61584/1‑138 (MQ=255)
                                                                                                                      gAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACa                                         <  3:61584/138‑1 (MQ=255)
                                                                                                                      gAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACa                                         <  3:57174/138‑1 (MQ=255)
                                                                                                                       aaTGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGcac                                                              <  8:216993/116‑1 (MQ=255)
                                                                                                                       aaTGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGcac                                                              >  7:216993/1‑116 (MQ=255)
                                                                                                                                    tGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCaa                                                                       <  2:31684/94‑1 (MQ=255)
                                                                                                                                    tGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCaa                                                                       >  1:31684/1‑94 (MQ=255)
                                                                                                                                       cGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCAATGATTCATGCCGTGacac              >  8:199471/1‑148 (MQ=255)
                                                                                                                                         gCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTgg                                                                                    <  1:205714/77‑1 (MQ=255)
                                                                                                                                         gCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTgg                                                                                    >  2:205714/1‑76 (MQ=255)
                                                                                                                                           cTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCAATGATTCATGCCGTGACACCAgcg         >  5:269682/1‑149 (MQ=255)
                                                                                                                                             gCAACTGTTGCAGCAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAg                                        >  3:70757/1‑116 (MQ=255)
                                                                                                                                             gCAACTGTTGCAGCAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAg                                        <  4:70757/116‑1 (MQ=255)
                                                                                                                                                aCTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCTAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCAATGATTCATGCCGTGACACCAGCGCCGa     >  6:301371/1‑148 (MQ=255)
                                                                                                                                                  tgctgcAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCAATGATTCATGCCGTGACACCAGCGCCGaaaa  >  7:234012/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCACGCCCATCGACCAGCCGTCGGAGACGATGTGGTGCAGGGTCACGAGCAGTAGGTGCTCCTGCTC‑‑GCCACGCTGATACAGCTTCATGCGAATCAGCGGCCCGTTGGCAAGGTCGAATGGCTGCTGGATGTCGGCCTGCAACTGCTGCAACAGTTGCGCGGGGTCAAAGGCTTGCGCCTGCTGCTGCTCAAGGGTGACCGCGAGATTGGCGTGGACGCGCAAGCCAGGCACACCTTGCTGCTCTTCGAAACAGGTACGCAATGATTCATGCCGTGACACCAGCGCCGAAAA  >  NC_002947/4790130‑4790422

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