Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,034,253 A→T 8.0% G312G (GGT→GGA) ‡ PP_2647 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,034,2530AT 8.0% 180.4 / 5.1 62G312G (GGT→GGA) ‡PP_2647MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (24/33);  new base T (2/3);  total (26/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.97e-01

AGCCCACGGCAATGATCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGCCACGCCATACAGC  >  NC_002947/3034105‑3034397
                                                                                                                                                    |                                                                                                                                                
aGCCCACGGCAATGATCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCa                                                                                                                                                  >  4:107019/1‑149 (MQ=255)
 gCCCACGGCAATGATCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGAGCCACGACAAACATACAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCAc                                                                                                                                                 >  3:202491/1‑149 (MQ=255)
     aCGGCAATGATCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGt                                                                                                                                             >  8:102292/1‑149 (MQ=255)
          aaTGATCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGc                                                                                                                                        <  5:44918/149‑1 (MQ=255)
              aTCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGACGGGTGCCGCCATGctga                                                                                                                                     <  7:252138/148‑4 (MQ=255)
               tCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   >  3:140994/1‑149 (MQ=255)
               tCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   >  8:242976/1‑149 (MQ=255)
               tCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   <  8:139649/149‑1 (MQ=255)
               tCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   <  6:31971/149‑1 (MQ=255)
               tCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   <  2:179400/149‑1 (MQ=255)
               tCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   <  3:107019/149‑1 (MQ=255)
                 aCCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGc                                                                                                                                 >  1:76822/1‑149 (MQ=255)
                  ccTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCg                                                                                                                                <  2:49621/149‑1 (MQ=255)
                   cTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGACGGGTGCCGCCATGct                                                                                                                                       >  6:172334/1‑140 (MQ=255)
                   cTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGACGGGTGCCGCCATGct                                                                                                                                       <  5:172334/141‑2 (MQ=255)
                     ggAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTgg                                                                                                                                  >  2:193904/1‑144 (MQ=255)
                     ggAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTgg                                                                                                                                  <  1:193904/144‑1 (MQ=255)
                      gAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGaa                                                                                                                              <  3:10673/147‑1 (MQ=255)
                      gAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGt                                                                                                                            >  4:52928/1‑149 (MQ=255)
                          aTGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCgg                                                                                                                        >  2:213036/1‑149 (MQ=255)
                           tGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGACGGGTGCCGCCATGct                                                                                                                                       <  4:305574/133‑2 (MQ=255)
                           tGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGACGGGTGCCGCCATGct                                                                                                                                       >  3:305574/1‑132 (MQ=255)
                            ggCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCggg                                                                                                                       >  5:288364/1‑148 (MQ=255)
                                     cAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCgg                                                                                                             <  8:7809/149‑1 (MQ=255)
                                        cGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGc                                                                                                          <  1:171535/149‑1 (MQ=255)
                                            aTGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGAc                                                                                                      <  5:311742/149‑1 (MQ=255)
                                               ccACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGc                                                                                                   >  8:123088/1‑149 (MQ=255)
                                                   gTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGAt                                                                                                >  6:223938/1‑148 (MQ=255)
                                                       gACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   >  1:202117/1‑109 (MQ=255)
                                                       gACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTg                                                                                                                                   <  2:202117/109‑1 (MQ=255)
                                                         cATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAgc                                                                                         >  3:98915/1‑149 (MQ=255)
                                                          aTCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAgcg                                                                                        >  2:96046/1‑149 (MQ=255)
                                                               tGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAgcgc                                                                                       <  7:102292/145‑1 (MQ=255)
                                                                  gcgcCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGaa                                                                                                                              >  2:323166/1‑103 (MQ=255)
                                                                  gcgcCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGaa                                                                                                                              <  1:323166/103‑1 (MQ=255)
                                                                  gcgcCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGc                                                                                <  4:140994/149‑1 (MQ=255)
                                                                  gcgcCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGc                                                                                <  4:98915/149‑1 (MQ=255)
                                                                    gcCACGACAAACATCCAGATATTGGCGCTAGCCGGGAACGCTGCGAAGGCAAAGCCCCATAGTGCTATCAGCATGTCGGCAACCGTCATGCCGGTGGCTAAGTGCTGGAACAGTGCTGTGCTGAGTGCTATCATGTGCGCGACCAGCAt                                                                              <  8:158231/149‑1 (MQ=255)
                                                                          aCAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAggg                                                                         >  5:66656/1‑148 (MQ=255)
                                                                          aCAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGt                                                                        <  3:52928/149‑1 (MQ=255)
                                                                            aaaCATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAgg                                                                          >  6:90500/1‑145 (MQ=255)
                                                                            aaaCATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAgg                                                                          <  5:90500/145‑1 (MQ=255)
                                                                             aaCATCCAGATTTTGGCGCAAGCCGGGAACGCACCGTAGGCATAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTTCGTGAACAGGGCGGTGGTGACTGGGATCATGAGCGCGACCAGCATCAGGGTATg                                                                     <  4:202491/149‑1 (MQ=255)
                                                                                  ccAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCa                                                                >  1:154288/1‑149 (MQ=255)
                                                                                   cAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCAc                                                               >  7:234880/1‑149 (MQ=255)
                                                                                                aaGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGgca                                                  <  6:288364/149‑1 (MQ=255)
                                                                                                     gggAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCaaa                                             <  3:45226/149‑1 (MQ=255)
                                                                                                             aacgcaGGAAAAGCCCCATAGTGAGATCAGCATGGCGGCACCCGTCATGCCGGTGCCGAAGTGCGGGAACAGGGCGGTGCTCACTGCGATCATGAGCGCGCCCCGCATCAGGGTATGCCGCACCCTGCGGTTGGCAGAAAAGCCGGCg                                      >  2:153652/6‑148 (MQ=255)
                                                                                                               cGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCg                                      >  6:107715/1‑146 (MQ=255)
                                                                                                                 aaGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGAt                                 <  2:76822/149‑1 (MQ=255)
                                                                                                                      aaaGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGc                            <  1:213036/149‑1 (MQ=255)
                                                                                                                           cccATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGAAGTTGCCCAt                        >  3:82915/1‑148 (MQ=255)
                                                                                                                            ccATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGt                                                       >  5:43320/1‑116 (MQ=255)
                                                                                                                            ccATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGt                                                       <  6:43320/116‑1 (MQ=255)
                                                                                                                            ccATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGa                      <  7:242976/149‑1 (MQ=255)
                                                                                                                             cATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGaa                     <  7:139974/149‑1 (MQ=255)
                                                                                                                                 gTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGTATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGcc                  <  1:96046/148‑1 (MQ=255)
                                                                                                                                    cGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGcc              <  5:107715/149‑1 (MQ=255)
                                                                                                                                    cGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGc               <  6:66656/148‑1 (MQ=255)
                                                                                                                                    cGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGGTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGc               <  8:234880/148‑1 (MQ=255)
                                                                                                                                            aTGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGCCACGCCATa      >  1:263621/1‑149 (MQ=255)
                                                                                                                                                cggcACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCTGCGAAGATGTTGCCCATGAAGCCGGCCACGCCATAcagc  <  5:223938/149‑1 (MQ=255)
                                                                                                                                                cggcACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGCCACGCCATAcagc  <  2:263621/149‑1 (MQ=255)
                                                                                                                                                    aCCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGccc                          >  2:253105/1‑121 (MQ=255)
                                                                                                                                                    aCCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGccc                          <  1:253105/121‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGCCCACGGCAATGATCACCTGGAACATGGCCACGAACAGCGGCATGCCACGTTCGACATCTTTGGGCGCCACGACAAACATCCAGATATTGGCGCAAGCCGGGAACGCACCGAAGGCAAAGCCCCATAGTGCGATCAGCATGGCGGCACCCGTCATGCCGGTGGCGAAGTGCGGGAACAGGGCGGTGCTGACTGCGATCATGAGCGCGACCAGCATCAGGGTATGCCGCACGCTGCGGTTGGCAGCAAAGCCGGCGAAGATGTTGCCCATGAAGCCGGCCACGCCATACAGC  >  NC_002947/3034105‑3034397

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: ‑