Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,580,944 C→A 15.2% K151N (AAG→AAT ttgB ← efflux pump membrane protein TtgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,580,9440CA15.2% 168.3 / 20.2 66K151N (AAG→AATttgBefflux pump membrane protein TtgB
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base A (7/3);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

CGGGGTCAGCTGGAACTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTCTGTACCTGTACCTGCGC  >  NC_002947/1580797‑1581090
                                                                                                                                                   |                                                                                                                                                  
cGGGGTCAGCTGGAACTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCt                                                                                                                                                   >  1:227917/1‑149 (MQ=255)
cGGGGTCAGCTGGAACTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCt                                                                                                                                                   <  2:225294/149‑1 (MQ=255)
  gggTCAGCTGGAACTTATTCTGCTTTGCCGGATCGAGCCAGATACTGATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCTGGTCGTCCTTg                                                                                                                                                 <  4:362072/149‑1 (MQ=255)
     gcAGCTCCAACTAATTCACCCCGGCCGCCACGAGCCACAGACGCAGGGCGTACTGCGCCCCGAACACCTGGCAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTc                                                                                                                                              <  4:325143/148‑1 (MQ=255)
         cTGGAACTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCAtgc                                                                                                                                          >  1:209521/1‑149 (MQ=255)
            gAACTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCAtgctgc                                                                                                                                       >  7:63220/1‑149 (MQ=255)
              aCTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCg                                                                                                                                     >  4:242447/1‑149 (MQ=255)
                     cAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCg                                                                                                                               <  8:67985/148‑1 (MQ=255)
                      aGCTTGGCCGGATCGAGCCAGATATGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGAATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCg                                                                                                                               >  2:30440/1‑147 (MQ=255)
                      aGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGga                                                                                                                             <  2:121267/149‑1 (MQ=255)
                       gCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGAGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGACATTGGTCATGCTGCCGTCTTCg                                                                                                                               >  7:331590/1‑146 (MQ=255)
                        aggtgCCGGATCGAGCCAGATACGCATGGCGTAGTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTGCGAGACGATAATGTTGGCGAGGTCGTCCTTGGTCATGCTGCCGTCTTCGgag                                                                                                                            <  2:269476/144‑1 (MQ=255)
                           ggCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGAc                                                                                                                          <  8:309732/147‑1 (MQ=255)
                           ggCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGAc                                                                                                                          <  5:139107/147‑1 (MQ=255)
                             ccGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGgatcacagg                                                                                                                      >  6:70720/1‑142 (MQ=255)
                             ccGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAgg                                                                                                                      >  4:24115/1‑149 (MQ=255)
                             ccGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAgg                                                                                                                      >  6:133364/1‑149 (MQ=255)
                              cGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGc                                                                                                                     <  1:303732/149‑1 (MQ=255)
                              cGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGc                                                                                                                     <  5:133364/149‑1 (MQ=255)
                                gATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCAtg                                                                                                                                           >  4:263252/1‑125 (MQ=255)
                                gATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCAtg                                                                                                                                           <  3:263252/125‑1 (MQ=255)
                                          gATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACcagca                                                                                                         <  3:230730/149‑1 (MQ=255)
                                                  tGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGttc                                                                                                 <  3:24115/149‑1 (MQ=255)
                                                    gCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGACGATCACCAGACGGAAGttct                                                                                                >  6:185691/1‑148 (MQ=255)
                                                    gCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGttctt                                                                                               <  7:200975/149‑1 (MQ=255)
                                                     cggaCTGCGCACAGAACACCTGGAAGTGACCTACACCCGCTGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCTGGAAGttctt                                                                                               <  8:331590/145‑1 (MQ=255)
                                                     cGTACTGCGCACCGAACACCTGGAAGTCACATACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGttctt                                                                                               >  3:191403/1‑148 (MQ=255)
                                                      gTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCa                                                                                             <  2:258828/149‑1 (MQ=255)
                                                              cACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTg                                                                                     >  2:256269/1‑149 (MQ=255)
                                                                       gcaggaggtgcagcagcaaggtatccgcgtcaccaaagcagtgaagaacttcctgctggtgatcGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGc                                                                                                                     >  3:381187‑M1/65‑108 (MQ=255)
                                                                       gcaggacgtgcagcagcaaggtatccgcgtcaccaaagcagtgaagaacttcctgctggtgatcGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGc                                                                                                                     <  4:381187‑M1/44‑1 (MQ=255)
                                                                        cTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGa                                                                           <  7:127773/149‑1 (MQ=255)
                                                                          ggAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATa                                                                         >  6:88754/1‑149 (MQ=255)
                                                                            aaGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATAcc                                                                       <  1:30440/149‑1 (MQ=255)
                                                                             gcccACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCGGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCt                                                                      <  2:307729/146‑1 (MQ=255)
                                                                                gcagcagcaaggtatccgcgtcaccaaagcagtgaagaacttcctgctggtgatcGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGa                                                                                                                  <  5:336613‑M1/47‑1 (MQ=255)
                                                                                gcagcagcaaggtatccgcgtcaccaaagcagtgaagaacttcctgctggtgatcGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGa                                                                                                                  >  6:336613‑M1/56‑102 (MQ=255)
                                                                                   cTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTtgctgc                                                                >  3:331140/1‑149 (MQ=255)
                                                                                   cTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTtgctgc                                                                >  3:327531/1‑149 (MQ=255)
                                                                                     acacCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTtgctgctg                                                              >  1:5913/1‑149 (MQ=255)
                                                                                     acacCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTtgctgctg                                                              >  8:128366/1‑149 (MQ=255)
                                                                                      cacCCGCGGTAAGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTtgctgct                                                               <  5:185691/147‑1 (MQ=255)
                                                                                         ccGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACt                                                          >  5:187018/1‑149 (MQ=255)
                                                                                               tACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGc                                                    <  2:209521/149‑1 (MQ=255)
                                                                                               tACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGc                                                    <  6:187018/149‑1 (MQ=255)
                                                                                                 cgcgAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGTCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCgg                                                  <  7:128039/149‑1 (MQ=255)
                                                                                                  gcgAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCtgctgc                                                 >  4:38524/1‑149 (MQ=255)
                                                                                                    gagaTCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGgcag                                               >  2:14981/1‑149 (MQ=255)
                                                                                                          gggTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCg                                          <  2:227917/148‑1 (MQ=255)
                                                                                                              ccTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCTGccgcaggaagtgc                                     <  3:38524/149‑8 (MQ=255)
                                                                                                                tGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGAcgcg                                                                             <  4:206753/107‑1 (MQ=255)
                                                                                                                tGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGAcgcg                                                                             >  3:206753/1‑107 (MQ=255)
                                                                                                                tGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTgg                                    >  7:29296/1‑148 (MQ=255)
                                                                                                                tGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGc                                   >  2:290738/1‑149 (MQ=255)
                                                                                                                   aTGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAg                                <  2:92224/149‑1 (MQ=255)
                                                                                                                           gaCGATATAGTTGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCATGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCt                        >  6:360790/1‑149 (MQ=255)
                                                                                                                             cGATATAGTCGGCCAGGTCGTCATTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTg                      >  7:285120/1‑149 (MQ=255)
                                                                                                                              gATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGt                     >  8:182143/1‑149 (MQ=255)
                                                                                                                               atatAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGtt                    <  7:128366/149‑1 (MQ=255)
                                                                                                                               atatAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGtt                    >  2:230095/1‑149 (MQ=255)
                                                                                                                                    gTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTctgt                <  1:256269/148‑1 (MQ=255)
                                                                                                                                     ttGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTctgtac              <  4:191403/149‑1 (MQ=255)
                                                                                                                                        GCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTctgtac              >  7:128395‑M1/1‑146 (MQ=255)
                                                                                                                                        GCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTATTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCAATTCCTGAGGCATCAGCGGGATGGCCAGGTTCATCTAGTtctgtacca            <  8:128395‑M1/148‑2 (MQ=255)
                                                                                                                                          CAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTctgtacctgta         >  4:331140‑M1/1‑149 (MQ=255)
                                                                                                                                           aGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCGCCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTctgtacctgtac        <  8:140558/149‑1 (MQ=255)
                                                                                                                                                   cTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTCTGTACCTGTACCTgcgc  <  3:26980/147‑1 (MQ=255)
                                                                                                                                                   cTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTCTGTACCTGTACCTgcgc  >  4:26980/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGGGTCAGCTGGAACTTATTCAGCTTGGCCGGATCGAGCCAGATACGCATGGCGTACTGCGCACCGAACACCTGGAAGTCACCTACACCCGCGGTACGCGAGATCGGGTCCTGCATGTTCGAGACGATATAGTTGGCCAGGTCGTCCTTGGTCATGCTGCCGTCTTCGGAGACCAGGCCGATCACCAGCAGGAAGTTCTTCACTGCTTTGGTGACGCGGATACCTTGCTGCTGCACTTCCTGCGGCAGCAGCGGGGTGGCCAGGTTCAGCTTGTTCTGTACCTGTACCTGCGC  >  NC_002947/1580797‑1581090

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: