Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,863,282 T→A 9.7% L309Q (CTG→CAG)  PP_1666 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,863,2820TA 9.7% 158.3 / 7.4 53L309Q (CTG→CAG) PP_1666hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base A (2/3);  total (26/27)
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.18e-01

GGTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCG  >  NC_002947/1863133‑1863430
                                                                                                                                                     |                                                                                                                                                     
ggTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcct                                                                                                                                                       >  1:683108/1‑150 (MQ=255)
 gTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctg                                                                                                                                                      >  1:1423846/1‑150 (MQ=255)
 gTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctg                                                                                                                                                      >  1:265065/1‑150 (MQ=255)
  tGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctgg                                                                                                                                                     >  1:1249497/1‑150 (MQ=255)
    cAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctggcc                                                                                                                                                   >  2:922984/1‑150 (MQ=255)
        tGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCACGc                                                                                                                                               <  2:1091593/150‑1 (MQ=255)
          cAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCTg                                                                                                                                             <  2:859304/150‑1 (MQ=255)
           aGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGc                                                                                                                                            <  2:1341015/150‑1 (MQ=255)
           aGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCTgg                                                                                                                                            <  2:957951/150‑2 (MQ=255)
            gggCATGAACCTGCAGCGCTATGCCGAGCTGGCGCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCACAGCTTGGCCTGGCCAAGCTGCa                                                                                                                                           <  2:1061485/150‑1 (MQ=255)
                aTGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGTCCTATGGCCCACGCTTGCAGGCTGCTGAAGGACCTCCCCTCACCCACGCGGTTAGTCGGAGTTGTGAACAGGTATGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGaa                                                                                                                                       <  1:387783/150‑1 (MQ=255)
                        gCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGt                                                                                                                               <  2:1069327/150‑1 (MQ=255)
                               tATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAg                                                                                                                        <  2:158711/150‑1 (MQ=255)
                                 tGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAgcg                                                                                                                      >  1:498909/1‑148 (MQ=255)
                                   ccGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCAc                                                                                                                    >  2:1200924/1‑150 (MQ=255)
                                         tGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCCGGCACcgccgg                                                                                                              >  2:843441/1‑150 (MQ=255)
                                         tGGCCCGCGTGCTTGAGCCCTATGGACCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACcggcgg                                                                                                              <  1:714995/150‑1 (MQ=255)
                                           gCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGc                                                                                                                                   >  1:401934/1‑127 (MQ=255)
                                           gCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGc                                                                                                                                   <  2:401934/127‑1 (MQ=255)
                                           gCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTg                                                                                                            >  1:219141/1‑150 (MQ=255)
                                                  tGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCCGGCACCGCCGGTGcctgctc                                                                                                     <  1:843441/150‑6 (MQ=255)
                                                     ttGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGc                                                                                                  <  2:219141/150‑1 (MQ=255)
                                                        aGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTc                                                                                                                              >  2:709646/1‑119 (MQ=255)
                                                        aGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTc                                                                                                                              <  1:709646/119‑1 (MQ=255)
                                                           ccTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                            >  1:1297965/1‑150 (MQ=255)
                                                            cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                            <  2:1091483/149‑1 (MQ=255)
                                                            cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                            >  1:1091483/1‑149 (MQ=255)
                                                            cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCAGGTCGAGCAGGCACCGGCTGCGTGGGCCATGCCATTAcc                                                                                           >  2:307635/1‑150 (MQ=255)
                                                                    aatgtttGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTgg                                                                                   >  1:410187/6‑150 (MQ=255)
                                                                     aCGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCTgg                                                                                                                                            <  2:510604/92‑2 (MQ=255)
                                                                     aCGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCTgg                                                                                                                                            >  1:510604/1‑91 (MQ=255)
                                                                       gCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGACATGCCAGTACCGGCAGTGGAcg                                                                                >  2:1038799/1‑150 (MQ=255)
                                                                        cTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGAcgc                                                                               <  2:550470/150‑1 (MQ=255)
                                                                             aGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAg                                                                          >  2:1132907/1‑150 (MQ=255)
                                                                                 ggcTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCGACGGGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCcgc                                                                      <  2:1328415/150‑1 (MQ=255)
                                                                                           aCTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCcgcg                                                            >  2:48340/1‑150 (MQ=255)
                                                                                              aCCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCt                                                         >  1:910661/1‑150 (MQ=255)
                                                                                                   gACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGCCGGTGcc                                                                                                          >  1:950552/1‑96 (MQ=255)
                                                                                                   gACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGCCGGTGcc                                                                                                          <  2:950552/96‑1 (MQ=255)
                                                                                                     ccTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                                  <  2:910661/150‑1 (MQ=255)
                                                                                                     ccTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                                  <  1:307635/150‑1 (MQ=255)
                                                                                                      cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagcc                                                 <  2:1297965/150‑1 (MQ=255)
                                                                                                      cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagcc                                                 >  2:919780/1‑150 (MQ=255)
                                                                                                        aCGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagccag                                               <  2:265065/150‑1 (MQ=255)
                                                                                                                 cGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGc                                      <  2:1423846/150‑1 (MQ=255)
                                                                                                                  gggCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAGGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGcc                                     >  1:803049/1‑150 (MQ=255)
                                                                                                                  gggCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGcc                                     <  1:1200924/150‑1 (MQ=255)
                                                                                                                    gCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                                   >  2:571810/1‑150 (MQ=255)
                                                                                                                    gCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCAGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                                   >  2:551154/1‑150 (MQ=255)
                                                                                                                          gTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACg                            >  1:430291/1‑150 (MQ=255)
                                                                                                                               cAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCG‑TCGAACGGCTTg                       <  1:48340/150‑1 (MQ=255)
                                                                                                                                   tACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGtt                   >  1:895516/1‑150 (MQ=255)
                                                                                                                                        gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGc                      <  2:1166798/142‑1 (MQ=255)
                                                                                                                                        gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGc                      >  1:1166798/1‑142 (MQ=255)
                                                                                                                                        ctgcAGCTTGGCCCGGACACGCCGCCGGACGTCTCTGCCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGGGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGc              <  2:373435/147‑1 (MQ=255)
                                                                                                                                          ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCAAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGTCTTGCGTTTTCACTg            >  2:227054/1‑150 (MQ=255)
                                                                                                                                           cAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGGGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCtag           >  1:373435/1‑148 (MQ=255)
                                                                                                                                               ttGGCCTGGCAAAGCTGGAGGGAGTGCCGGTCTATCATGCACCTGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCAGCGCCGCCCGAACCACCGCAGCTAGCACGCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAg       >  2:1020792/1‑150 (MQ=255)
                                                                                                                                                 ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGa     <  2:895516/150‑1 (MQ=255)
                                                                                                                                                   cctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAgc   >  1:1024745/1‑150 (MQ=255)
                                                                                                                                                    ctggccAGGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAgcg  <  2:803049/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCG  >  NC_002947/1863133‑1863430

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