Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,181,223 T→C 5.9% H140H (CAT→CAC)  PP_2791 → aminoglycoside phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,181,2230TC 5.9% 199.7 / 2.0 68H140H (CAT→CAC) PP_2791aminoglycoside phosphotransferase
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/2);  total (33/35)
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.55e-01

CGGCGTTGGGTACTGCGCAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCG  >  NC_002947/3181077‑3181371
                                                                                                                                                  |                                                                                                                                                    
cGGCGTTGGGTACTGCGCAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATgg                                                                                                                                                    >  1:233056/1‑149 (MQ=255)
       gggTACTGCGCAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGgt                                                                                                                                              <  3:326967/148‑1 (MQ=255)
         gTACTGCGCAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTTGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGtt                                                                                                                                           <  1:239583/149‑1 (MQ=255)
          tACTGCGCAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGttc                                                                                                                                          >  7:90216/1‑149 (MQ=255)
              gcgcAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTcc                                                                                                                                      >  3:233791/1‑149 (MQ=255)
              gcgcAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTcc                                                                                                                                      >  8:191814/1‑149 (MQ=255)
               cgcAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGttcttc                                                                                                                                       <  2:113373/147‑1 (MQ=255)
               cgcAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGttcttc                                                                                                                                       <  7:60478/147‑1 (MQ=255)
               cgcAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCt                                                                                                                                     <  6:55644/149‑1 (MQ=255)
                  aGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGAt                                                                                                                                  <  7:339664/149‑1 (MQ=255)
                  aGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCACGGGCGGGTGTTCTTCCTGAt                                                                                                                                  >  7:133799/1‑149 (MQ=255)
                  aGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCACGGGCGGGTGTTCTTCCTGAt                                                                                                                                  >  7:128921/1‑149 (MQ=255)
                      aTGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAg                                                                                                                              >  6:123403/1‑149 (MQ=255)
                      aTGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAg                                                                                                                              >  8:241462/1‑149 (MQ=255)
                             ttCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGc                                                                                                                       >  1:63122/1‑149 (MQ=255)
                                 gcgcTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCggg                                                                                                                   <  1:141815/149‑1 (MQ=255)
                                   gcTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTa                                                                                                                 <  1:67494/149‑1 (MQ=255)
                                      aCCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAgcg                                                                                                              >  3:288745/1‑149 (MQ=255)
                                       cccGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAgcg                                                                                                              <  1:315867/148‑1 (MQ=255)
                                       cccGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAgcg                                                                                                              <  7:128916/148‑1 (MQ=255)
                                           ccAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTgcc                                                                                                         <  7:191814/149‑1 (MQ=255)
                                           ccAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTgcc                                                                                                         <  4:253647/149‑1 (MQ=255)
                                           ccAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTgcc                                                                                                         <  7:241462/149‑1 (MQ=255)
                                               ccTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTgccggcc                                                                                                     <  5:123403/149‑1 (MQ=255)
                                                 ttGACAGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGACACGCCAGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGAGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCcgc                                                                                                   >  3:250494/1‑149 (MQ=255)
                                                   gACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCcgc                                                                                                   <  2:51732/147‑1 (MQ=255)
                                                      cGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCt                                                                                              <  6:150015/149‑1 (MQ=255)
                                                       gTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCt                                                                                              >  6:40056/1‑148 (MQ=255)
                                                          aaCAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAg                                                                                           >  4:28663/1‑148 (MQ=255)
                                                          aaCAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAg                                                                                           <  5:40056/148‑1 (MQ=255)
                                                          aaCAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAg                                                                                           >  7:98780/1‑148 (MQ=255)
                                                          aaCAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCACGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAg                                                                                           <  8:128921/148‑1 (MQ=255)
                                                          aaCAGGAGTTCGCCGGGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCACGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAg                                                                                           <  8:133799/148‑1 (MQ=255)
                                                            cAGGAGTTCGCCGTGATGCAGGTGGTTAACCAGCCCGGCGTGAAAGTGCAACGCCCGCTCTGGCTGTGCCGCGATGTCCGCGTGCATGGGCGGGTGTTCTTCCTGATCGAGTATGTGCCGGCTAGCGCTGCCGGCCGCGCGCTCAGCAc                                                                                        >  3:250282/1‑149 (MQ=255)
                                                             aGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTAt                                                                                                                           >  1:108654/1‑113 (MQ=255)
                                                             aGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTAt                                                                                                                           <  2:108654/113‑1 (MQ=255)
                                                             aGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCAGGGTAGCGCTGCCGGCCGCGCGCTCAGCAcc                                                                                       >  1:237201/1‑149 (MQ=255)
                                                                     gCCGTGCTGCAGGTGGTTTACCAGACCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCgg                                                                               >  3:67267/1‑149 (MQ=255)
                                                                      ccGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGt                                                                              >  4:31931/1‑149 (MQ=255)
                                                                       cGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTc                                                                             <  7:95634/149‑1 (MQ=255)
                                                                              cAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAggg                                                                      <  4:233791/149‑1 (MQ=255)
                                                                              cAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAg                                                                        <  5:13631/147‑1 (MQ=255)
                                                                                  tggtTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCgg                                                                  <  3:28663/149‑1 (MQ=255)
                                                                                    gtTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTc                                                                             >  8:164991/1‑136 (MQ=255)
                                                                                    gtTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTc                                                                             <  7:164991/136‑1 (MQ=255)
                                                                                    gtTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGgc                                                                <  2:233056/149‑1 (MQ=255)
                                                                                        aCCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTgccgg                                                                                                       >  6:150950/1‑106 (MQ=255)
                                                                                         ccAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCaa                                                            <  2:237201/148‑1 (MQ=255)
                                                                                          cAGGCCGGCATGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACt                                                          >  5:37528/1‑149 (MQ=255)
                                                                                              ccGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCg                                                      >  8:261697/1‑149 (MQ=255)
                                                                                                   gTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACg                                                   >  1:239295/1‑147 (MQ=255)
                                                                                                            ccACGCCCTCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTACTGATGGAGTATGTGCCGGGTAGCGCGGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGAAGCTTgtc                                          <  4:250494/147‑3 (MQ=255)
                                                                                                            ccACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGGGCTCAGGGCCGGGCGCAACTGGCGACGCGGCTTGgcgg                                        >  4:261976/1‑148 (MQ=255)
                                                                                                                          gCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACTGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTc                          <  8:90216/149‑1 (MQ=255)
                                                                                                                          gCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGgcgc                                                                                  >  3:238748/1‑93 (MQ=255)
                                                                                                                          gCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGgcgc                                                                                  <  4:238748/93‑1 (MQ=255)
                                                                                                                          gCTGTGCCGCGATGGGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGt                           <  3:261976/148‑1 (MQ=255)
                                                                                                                            tgtgCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTc                                                                          >  7:185955/1‑99 (MQ=255)
                                                                                                                            tgtgCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTc                                                                          <  8:185955/99‑1 (MQ=255)
                                                                                                                                ccGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCAt                     >  6:100167/1‑148 (MQ=255)
                                                                                                                                  gcgATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAg                  <  8:98780/149‑1 (MQ=255)
                                                                                                                                    gATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTAGCGCGTCTGCATCAGGt                >  3:298883/1‑149 (MQ=255)
                                                                                                                                      tgtgCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACt                                                          <  6:149312/105‑1 (MQ=255)
                                                                                                                                      tgtgCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACt                                                          >  5:149312/1‑105 (MQ=255)
                                                                                                                                      tgtgCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCCACCTGGCGCGTCTGCATCAGGTcc              >  5:245084/1‑149 (MQ=255)
                                                                                                                                       gtgCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCgg                                                                  >  6:221228/1‑96 (MQ=255)
                                                                                                                                       gtgCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCgg                                                                  <  5:221228/96‑1 (MQ=255)
                                                                                                                                               gCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGt     >  8:213189/1‑149 (MQ=255)
                                                                                                                                                  tGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTgcg  <  3:31931/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGCGTTGGGTACTGCGCAGCGATGCGCTTTCAGCGCTACCCGCCAGCCTTGACCGTGAACAGGAGTTCGCCGTGCTGCAGGTGGTTTACCAGGCCGGCGTGAAAGTGCCACGCCCGCTCTGGCTGTGCCGCGATGTGCGCGTGCATGGGCGGGTGTTCTTCCTGATGGAGTATGTGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCAACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCG  >  NC_002947/3181077‑3181371

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