Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,407,056 C→A 6.6% V121L (GTG→TTG)  PP_4751 ← amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,407,0560CA 6.6% 190.6 / 3.9 61V121L (GTG→TTG) PP_4751amino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base A (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

TGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTC  >  NC_002947/5406914‑5407197
                                                                                                                                              |                                                                                                                                             
tgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCGAACGTGGGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAg                                                                                                                                         >  8:309938/1‑149 (MQ=255)
 ggCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAgg                                                                                                                                        <  7:255367/149‑1 (MQ=255)
    caccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGc                                                                                                                                     >  1:360100/1‑149 (MQ=255)
    caccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAgcctg                                                                                                                                     <  6:211538/149‑5 (MQ=255)
       caGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCggg                                                                                                                                  <  6:127/149‑1 (MQ=255)
        agagCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCggg                                                                                                                                  >  1:85962/1‑148 (MQ=255)
         gagCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGaa                                                                                                                                >  1:61610/1‑149 (MQ=255)
         gagCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGaa                                                                                                                                >  8:163926/1‑149 (MQ=255)
              gcgcgcCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGt                                                                                                                           <  8:68855/149‑1 (MQ=255)
              gcgcgcCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGt                                                                                                                           >  3:231372/1‑149 (MQ=255)
               cgcgcCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGtt                                                                                                                          <  8:216559/149‑1 (MQ=255)
               cgcgcCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGtt                                                                                                                          <  7:151272/149‑1 (MQ=255)
                         aaTGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTc                                                                                                                                           <  4:240873/122‑1 (MQ=255)
                         aaTGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTc                                                                                                                                           >  3:240873/1‑122 (MQ=255)
                         aaTGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTg                                                                                                                 >  5:277816/1‑148 (MQ=255)
                                ggCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACAcc                                                                                                          <  6:277816/148‑1 (MQ=255)
                                 gCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACAGTCAGGTGCGGGAACAGGTTGAACGACTGGAacacaa                                                                                                         >  1:271828/1‑146 (MQ=255)
                                  ccTTGTCAGCCAGACCAAAGCGTTCGAGCAGTTGAAAGCAGCGGGAGTCGGCTTCATTCGCGCTGCGGCCCAAGGTGCGGATTTTCGCCAGGCTGACGTATCGCAGCACCGACAGGTGCGGCAACAGGTTGACCCACTGGAACCCCAtc                                                                                                       >  8:68370/1‑148 (MQ=255)
                                      gTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGaa                                                                                                                                >  2:58704/1‑120 (MQ=255)
                                      gTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGaa                                                                                                                                <  1:58704/120‑1 (MQ=255)
                                         agccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGGCGGCTTCTTTGGCCCTGCGGCCCAAGCTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCACTgc                                                                                                <  4:352410/149‑1 (MQ=255)
                                          gccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATgcg                                                                                               >  5:142375/1‑149 (MQ=255)
                                           ccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACt                                                                                                                  <  8:91987/129‑1 (MQ=255)
                                           ccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACt                                                                                                                  >  7:91987/1‑129 (MQ=255)
                                             agccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCATCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCg                                                                                            >  6:75146/1‑149 (MQ=255)
                                              gccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGACCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGGACAGGTTGAAAGACTGGAACACCAAACCAATGCGCCgg                                                                                           >  6:23157/1‑149 (MQ=255)
                                               cccAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGgc                                                                                          >  1:55731/1‑149 (MQ=255)
                                                ccAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGAGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGgcg                                                                                         >  5:118518/1‑149 (MQ=255)
                                                 cAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATgcgc                                                                                              >  3:332384/1‑143 (MQ=255)
                                                 cAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATgcgc                                                                                              <  4:332384/143‑1 (MQ=255)
                                                 cAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGgcgc                                                                                        <  7:309938/149‑1 (MQ=255)
                                                  aaCGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGgcgg                                                                                        >  4:237903/1‑147 (MQ=255)
                                                       gTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACtt                                                                                   >  4:325175/1‑148 (MQ=255)
                                                       gTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCATCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTc                                                                                  <  5:75146/149‑1 (MQ=255)
                                                        ttCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTca                                                                                 <  2:61610/149‑1 (MQ=255)
                                                        ttCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTAGAACACCATACCAATGCGCCGGCGCAACTTca                                                                                 <  7:187417/149‑1 (MQ=255)
                                                           gAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTcag                                                                                >  2:290434/1‑147 (MQ=255)
                                                             gCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGt                                                                            >  2:231692/1‑149 (MQ=255)
                                                                                    gCTTGTTTGGCGCTGCGGCCCATGGTGGGGATTTGCGCCAGGGTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACAGCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGcc                                                     <  6:118518/149‑1 (MQ=255)
                                                                                    gCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGcc                                                     <  6:142375/149‑1 (MQ=255)
                                                                                      ttcttTGGCGCAGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAg                                                   >  7:216673/1‑149 (MQ=255)
                                                                                             gcgcTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGAt                                            <  2:85962/149‑1 (MQ=255)
                                                                                                      cccAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGcc                                   <  7:68370/149‑1 (MQ=255)
                                                                                                      cccAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGc                                    >  1:257329/1‑148 (MQ=255)
                                                                                                                   tttGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGAc                       <  8:216673/148‑1 (MQ=255)
                                                                                                                        gcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGcaatcc                               >  7:220522/1‑135 (MQ=255)
                                                                                                                        gcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGccc                 <  1:290434/149‑1 (MQ=255)
                                                                                                                        gcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGccc                 >  8:124287/1‑149 (MQ=255)
                                                                                                                          cAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCt               <  3:325175/149‑1 (MQ=255)
                                                                                                                            ggCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTg             <  2:257329/149‑1 (MQ=255)
                                                                                                                             ccTGACGGTGCTGCGCAACGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAg                                                   <  4:170105/109‑1 (MQ=39)
                                                                                                                             ccTGACGGTGCTGCGCAACGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAg                                                   >  3:170105/2‑110 (MQ=39)
                                                                                                                              cTGACGGTGCTGCGCAACGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGAt                                            >  2:1232/1‑116 (MQ=255)
                                                                                                                                 aCGTTGCGCTGGACCGTCAGGTGCGGGACCAGGTGGAACGACTGTAACACCATACCAATCCGCCGGCGCAACTTCAGCAGCTCTGTACGGGGCAACGCCTTGGCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcgg        <  2:271828/149‑1 (MQ=255)
                                                                                                                                   gTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGCCTGGAACACCATACCAATACGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcggcg      >  1:279829/1‑149 (MQ=255)
                                                                                                                                    ttGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcggcgg     <  1:231692/149‑1 (MQ=255)
                                                                                                                                      gcgcAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATg                            >  4:275840/1‑124 (MQ=255)
                                                                                                                                      gcgcAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATg                            <  3:275840/124‑1 (MQ=255)
                                                                                                                                      gcgcAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCt   >  1:340311/1‑149 (MQ=255)
                                                                                                                                        gcagcaCCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTc  >  2:10074/1‑148 (MQ=255)
                                                                                                                                         cagcaCCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTc  >  1:282293/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTC  >  NC_002947/5406914‑5407197

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