Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,922,769 C→T 7.8% A828A (GCG→GCA)  PP_2563 ← antibiotic biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,922,7690CT 7.8% 186.7 / 4.4 65A828A (GCG→GCA) PP_2563antibiotic biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

AGCGGGTGCCTGAAAACATGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATGGTG  >  NC_002947/2922623‑2922906
                                                                                                                                                   |                                                                                                                                         
aGCGGGTGCCTGAAAACATGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGc                                                                                                                                         <  3:20726/149‑1 (MQ=255)
 gCGGGTGCCTGAAAACATGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCt                                                                                                                                        <  6:104332/149‑1 (MQ=255)
      tGCCTGAAAACATGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCGATGCTGGCac                                                                                                                                   <  6:248418/149‑2 (MQ=255)
         cTGAAAACAGGCTCCGCTACCGACTCGACCTGGAAATGCGTCAGGTAGATATACAAGGGTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGc                                                                                                                                <  5:38585/149‑1 (MQ=255)
           gAAAACATGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGa                                                                                                                              >  3:113432/1‑149 (MQ=255)
                  tGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGAGGGGGAACCAGACAAGTGCCAGCAAAGCTGGCAGGGCGATGTCg                                                                                                                         >  1:114674/1‑148 (MQ=255)
                       cacTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAgt                                                                                                                  <  6:313518/147‑1 (MQ=255)
                       cGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAgt                                                                                                                  <  2:298574/149‑1 (MQ=255)
                            ccGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTa                                                                                                             <  4:113432/149‑1 (MQ=255)
                                   gAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCtt                                                                                                      >  3:291871/1‑149 (MQ=255)
                                       tGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCg                                                                                                                         >  1:125149/1‑126 (MQ=255)
                                       tGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCg                                                                                                                         <  2:125149/126‑1 (MQ=255)
                                         gAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTgcca                                                                                                  <  1:67630/147‑1 (MQ=255)
                                           aaTGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTgccagcca                                                                                              <  3:17416/149‑1 (MQ=255)
                                                gTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACAcc                                                                                         >  3:313346/1‑149 (MQ=255)
                                                        gATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAgc                                                                                  <  3:242047/148‑1 (MQ=255)
                                                          tataCAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAgcgc                                                                                <  3:53363/148‑1 (MQ=255)
                                                          tataCAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAgcgc                                                                                >  6:229883/1‑148 (MQ=255)
                                                              cAAGGTTGACGAAGCGATCACCTCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCa                                                                           <  3:199652/149‑1 (MQ=255)
                                                              cAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATc                                                                            <  4:313346/148‑1 (MQ=255)
                                                              cAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACCGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCAt                                                                             >  7:233092/1‑147 (MQ=255)
                                                                 ggTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAg                                                                          <  6:155453/147‑1 (MQ=255)
                                                                 ggTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAg                                                                          >  5:155453/1‑147 (MQ=255)
                                                                  gTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTc                                                                       >  7:92431/1‑149 (MQ=255)
                                                                        gAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTc                                                                 >  3:13163/1‑149 (MQ=255)
                                                                        gAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTc                                                                 >  7:71045/1‑149 (MQ=255)
                                                                        gAAGCGATCACCGCACCTGCCCCTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGCCGATGTCGAGATACTGAGAAAACGCTTGCCAGCCAACACAGGCAAGCGCCATCAGGTCCAACTc                                                                 >  7:205074/1‑149 (MQ=255)
                                                                            cGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGcc                                                              >  7:44534/1‑148 (MQ=255)
                                                                              aTCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAac                                                           >  5:256823/1‑149 (MQ=255)
                                                                                      aCCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCa                                                   >  2:274906/1‑149 (MQ=255)
                                                                                       ccTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCAc                                                  >  2:198490/1‑149 (MQ=255)
                                                                                          gCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCg                                                                <  7:56646/132‑1 (MQ=255)
                                                                                          gCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCg                                                                >  8:56646/1‑132 (MQ=255)
                                                                                          gCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGc                                               >  4:222144/1‑149 (MQ=255)
                                                                                          gCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGc                                               >  6:202933/1‑149 (MQ=255)
                                                                                           cccGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGc                                               >  2:322103/1‑148 (MQ=255)
                                                                                           cccGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGc                                               <  6:256823/148‑1 (MQ=255)
                                                                                            ccGTGCAACCAGGGCGGGTACCGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCca                                             <  8:233092/149‑1 (MQ=255)
                                                                                              gTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCcact                                           <  4:291871/149‑1 (MQ=255)
                                                                                               tGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCcactg                                          <  1:198490/149‑1 (MQ=255)
                                                                                               tGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCcactg                                          <  2:114674/149‑1 (MQ=255)
                                                                                                 cAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCg                                                                                        <  2:135210/101‑1 (MQ=255)
                                                                                                 cAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCg                                                                                        >  1:135210/1‑101 (MQ=255)
                                                                                                  aaCCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGc                                               <  8:130960/141‑1 (MQ=255)
                                                                                                  aaCCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCcagatcgg                                       >  7:130960/1‑143 (MQ=255)
                                                                                                    ccAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCcactgcactg                                     <  8:44534/149‑1 (MQ=255)
                                                                                                                 tGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCAGAAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTCGCCACCCATTTGc                        >  8:144920/1‑149 (MQ=255)
                                                                                                                   gCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGt                      <  4:13163/149‑1 (MQ=255)
                                                                                                                     aCCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGt                      >  6:205239/1‑147 (MQ=255)
                                                                                                                            cGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACGAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCgt              >  8:226631/1‑149 (MQ=255)
                                                                                                                            cGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCgt              <  8:71045/149‑1 (MQ=255)
                                                                                                                             ggAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCgtg             >  7:75134/1‑149 (MQ=255)
                                                                                                                              gAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCgtgt            <  1:322103/149‑1 (MQ=255)
                                                                                                                               aaCCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTc           <  1:274906/149‑1 (MQ=255)
                                                                                                                                aCCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAg                                                                          <  7:268060/85‑1 (MQ=255)
                                                                                                                                aCCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAg                                                                          >  8:268060/1‑85 (MQ=255)
                                                                                                                                aCCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCg          <  3:222144/149‑1 (MQ=255)
                                                                                                                                      cAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATgg    <  8:92431/149‑1 (MQ=255)
                                                                                                                                       aaGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATGGt   >  1:336446/1‑149 (MQ=255)
                                                                                                                                        cctGCCAGCGATGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTga                    <  5:312003/129‑1 (MQ=255)
                                                                                                                                        cctGCCAGCGATGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTga                    >  6:312003/3‑131 (MQ=255)
                                                                                                                                        cctGCCAGCGATGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATGGTg  >  2:304843/3‑149 (MQ=255)
                                                                                                                                         gTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATGGTg  >  8:168376/1‑148 (MQ=255)
                                                                                                                                         gTGCCAGCATCGCTGGAAGGGCGATGTCGACATAGTGTGTAAACGCTTGACAGCCAACAACGGCAAGCGCCATCAGGTCCAACTCGCCAACCACAAGCACGGCCACTGAACTGGCCACCCATTTGCGTGAGATCGTGTCGGCaggttg  >  2:325873/1‑143 (MQ=255)
                                                                                                                                         ctGCCAGCGATGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATGGTg  >  7:310288/2‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
AGCGGGTGCCTGAAAACATGCTCCGCTACCGACTCGAACTGGAAATGCGTCAGGTAGATATACAAGGTTGACGAAGCGATCACCGCACCTGCCCGTGCAACCAGGGCGGGTACTGGCACCGA‑GCGGAACCAGACAAGTGCCAGCATCGCTGGCAGGGCGATGTCGACATAGTGTGTAAACGCTTGCCAGCCAACACCGGCAAGCGCCATCAGGTCCAACTCGCCAACCACCAGCACGGCCACTGCACTGGCCACCCATTTGCGTGAGATCGTGTCGGCATGGTG  >  NC_002947/2922623‑2922906

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