Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 874,698 A→G 7.0% I93V (ATC→GTC) ‡ PP_0759 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947874,6980AG 7.0% 205.1 / 4.7 71I93V (ATC→GTC) ‡PP_0759hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (2/3);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCAGCTCAAGCT  >  NC_002947/874562‑874840
                                                                                                                                        |                                                                                                                                              
cATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcgg                                                                                                                                    <  2:126750/149‑2 (MQ=255)
cATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGgcgc                                                                                                                                    <  2:20717/149‑1 (MQ=255)
cATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGgcgc                                                                                                                                    <  2:255730/149‑1 (MQ=255)
       tcGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggt                                                                                                                                   <  6:39452/143‑3 (MQ=255)
       tcGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggt                                                                                                                                   >  5:39452/1‑141 (MQ=255)
            cGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGcc                                                                                                                        >  7:9177/1‑149 (MQ=255)
             gATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGccc                                                                                                                       >  1:78246/1‑149 (MQ=255)
             gATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGccc                                                                                                                       <  2:78246/149‑1 (MQ=255)
                 aGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGAt                                                                                                                   <  1:51059/149‑1 (MQ=255)
                  ggTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATg                                                                                                                  >  5:114177/1‑149 (MQ=255)
                   gTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGccc                                                                                                                                          <  4:177372/124‑1 (MQ=255)
                   gTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGccc                                                                                                                                          >  3:177372/1‑124 (MQ=255)
                      aaCCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTTTGATCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGACATTGCCAAAACCCTGATCGCCAAGTCGCCCAAATTACGCCCGGATGCagc                                                                                                              >  3:196598/1‑149 (MQ=255)
                      aaCCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCcgc                                                                                                              >  5:266547/1‑149 (MQ=255)
                       aCCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCcgcg                                                                                                             <  3:148179/149‑1 (MQ=255)
                             tGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTgg                                                                                                       <  5:288169/149‑1 (MQ=255)
                               tACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTggg                                                                                                      <  6:242395/148‑1 (MQ=255)
                                    gcAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTgg                                                                                                 >  5:174766/1‑148 (MQ=255)
                                    gcAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGc                                                                                                <  2:191721/149‑1 (MQ=255)
                                      aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggttg                                                                                                                                 >  4:273158/1‑110 (MQ=255)
                                      aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGGGCGATCAGgcggttg                                                                                                                                 <  3:273158/114‑5 (MQ=255)
                                      aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGc                                                                                                >  4:49379/1‑147 (MQ=255)
                                      aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCaa                                                                                              <  4:49192/149‑1 (MQ=255)
                                      aGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCaa                                                                                              <  8:9177/149‑1 (MQ=255)
                                               agcagACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGt                                                                                      <  4:46064/148‑1 (MQ=255)
                                                gcagACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGa                                                                                    <  3:138910/149‑1 (MQ=255)
                                                gcagACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGa                                                                                    >  4:99921/1‑149 (MQ=255)
                                                      ccTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATc                                                                                  >  3:154182/1‑145 (MQ=255)
                                                      ccTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCaa                                                                                >  4:282091/1‑147 (MQ=255)
                                                          ggTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtg                                                                           <  3:49379/148‑1 (MQ=255)
                                                          ggTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtg                                                                           <  6:174766/148‑1 (MQ=255)
                                                          ggTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtg                                                                           <  3:282091/148‑1 (MQ=255)
                                                           gTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGtgat                                                                         >  8:72987/1‑149 (MQ=255)
                                                                gCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGAtgaact                                                                    >  4:88385/1‑149 (MQ=255)
                                                                   tCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGAtgaactgaa                                                                 <  6:266547/149‑1 (MQ=255)
                                                                   tCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGAtgaactgaa                                                                 >  8:196094/1‑149 (MQ=255)
                                                                      aGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGc                                                              >  2:273331/1‑149 (MQ=255)
                                                                       gCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGcc                                                             <  6:114177/149‑1 (MQ=255)
                                                                              tgtgCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCg                                                      >  3:243556/1‑149 (MQ=255)
                                                                              tgtgCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCg                                                      >  7:17713/1‑149 (MQ=255)
                                                                                   tCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGAATGCCGCGCAGTGGGACTgg                                                                                                 >  8:102088/1‑101 (MQ=255)
                                                                                   tCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGAATGCCGCGCAGTGGGACTgg                                                                                                 <  7:102088/101‑1 (MQ=255)
                                                                                     gACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCggg                                               <  4:243556/149‑1 (MQ=255)
                                                                                       cAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCCGCCCGggcg                                             >  7:191322/1‑149 (MQ=255)
                                                                                          gTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAaga                                                                              <  1:203405/113‑1 (MQ=255)
                                                                                          gTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAaga                                                                              >  2:203405/1‑113 (MQ=255)
                                                                                           tCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                            >  5:138036/1‑146 (MQ=255)
                                                                                           tCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                            <  6:138036/146‑1 (MQ=255)
                                                                                              aGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                            <  4:195606/143‑1 (MQ=255)
                                                                                              aGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGggcgg                                            >  3:195606/1‑143 (MQ=255)
                                                                                              aGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCATGa                                       >  2:137627/1‑148 (MQ=255)
                                                                                               gcgcCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCACg                                                 >  4:185916/1‑137 (MQ=255)
                                                                                                          aaGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTacac                          >  5:128063/1‑149 (MQ=255)
                                                                                                          aaGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGCCGCCCAAATTCCGCACGGATGCCGCGCACTGGGACTGGAAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTacac                          >  8:185995/1‑149 (MQ=255)
                                                                                                             cGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCgg                       >  8:219760/1‑149 (MQ=255)
                                                                                                             cGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCgg                       <  3:88385/149‑1 (MQ=255)
                                                                                                              gtgtGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTac                            <  1:303034/143‑1 (MQ=255)
                                                                                                              gtgtGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTac                            >  2:303034/1‑143 (MQ=255)
                                                                                                                   cAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTgat                 >  5:131503/1‑149 (MQ=255)
                                                                                                                   cAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTga                  <  7:105085/148‑1 (MQ=255)
                                                                                                                      gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCCACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGc                                                    >  5:279995/1‑111 (MQ=255)
                                                                                                                      gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGc                                                    <  6:279995/111‑1 (MQ=255)
                                                                                                                      gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGa                                       >  6:259880/1‑124 (MQ=255)
                                                                                                                      gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGa                                       <  5:259880/124‑1 (MQ=255)
                                                                                                                      gCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTgatcg               >  3:90799/1‑148 (MQ=255)
                                                                                                                        cATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTgatcgat             >  4:87702/1‑148 (MQ=255)
                                                                                                                         aTTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCATGATCATCGTCTACACCGGGCTGATCGATCa           <  1:137627/149‑1 (MQ=255)
                                                                                                                             ccAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctc       <  4:154182/149‑1 (MQ=255)
                                                                                                                                aCCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctcaag    <  1:273331/149‑1 (MQ=255)
                                                                                                                                 ccGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctcaagc   >  3:159230/1‑149 (MQ=255)
                                                                                                                                   gCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCagctcaagct  >  8:131415/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                              
CATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGCCAACCGCCTGATCGCCCAGGCGCCCAAATTCCGCCCGGATGCCGCGCAGTGGGACTGGCAAGTCAACGTGATCAAGAGTGATGAACTGAACGCCAACTGCGGCCCGGGCGGCAAGATCATCGTCTACACCGGGCTGATCGATCAGCTCAAGCT  >  NC_002947/874562‑874840

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