Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 874,538 C→T 6.5% V39V (GTC→GTT)  PP_0759 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947874,5380CT 6.5% 227.8 / 2.7 76V39V (GTC→GTT) PP_0759hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (2/3);  total (38/38)
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.30e-01

CACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCAT  >  NC_002947/874393‑874684
                                                                                                                                                 |                                                                                                                                                  
cACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGaa                                                                                                                                                 <  5:213433/149‑1 (MQ=255)
  cTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCTTCCTTCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACg                                                                                                                                               <  5:251551/149‑1 (MQ=255)
     gCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCa                                                                                                                                            <  2:21381/149‑1 (MQ=255)
       gACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTTCGACGCCAAc                                                                                                                                          <  1:69764/149‑1 (MQ=255)
               cTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCAGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATg                                                                                                                                  >  5:309129/1‑149 (MQ=255)
                    aTACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATg                                                                                                                                  >  4:230302/1‑144 (MQ=255)
                    aTACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATg                                                                                                                                  <  3:230302/144‑1 (MQ=255)
                     tACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGc                                                                                                                            <  3:327299/149‑1 (MQ=255)
                      aCGGCGCATGAGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCa                                                                                                                           >  6:110575/1‑149 (MQ=255)
                           gcATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCAGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGctc                                                                                                                      <  6:309129/149‑1 (MQ=255)
                             aTGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCtctc                                                                                                                    >  1:110730/1‑149 (MQ=255)
                             aTGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGAGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCtctc                                                                                                                    <  1:296587/149‑1 (MQ=255)
                                 ggATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCAt                                                                                                                          >  8:290374/1‑139 (MQ=255)
                                 ggATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCAt                                                                                                                          <  7:290374/139‑1 (MQ=255)
                                       ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                          <  5:110575/149‑1 (MQ=255)
                                       ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                          >  3:119112/1‑149 (MQ=255)
                                        gAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGctct                                                                                                                     >  1:238790/1‑137 (MQ=255)
                                        gAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGctct                                                                                                                     <  2:238790/137‑1 (MQ=255)
                                        gAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGATGGCCGGCTGCCAGGCGGTGAATACAACCAGAGGAGGTGCATTTGGCGTCCAACGCCAACAGTACATGTTCAACATGCACTCGACCGATGAg                                                                                                          >  5:59895/1‑148 (MQ=255)
                                         aGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                             <  2:322767/144‑1 (MQ=255)
                                         aGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                             >  1:322767/1‑144 (MQ=255)
                                         aGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGt                                                                                                        >  4:2230/1‑149 (MQ=255)
                                          gtgtTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAAc                                                                                                                                          <  7:290351/114‑1 (MQ=255)
                                          gtgtTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAAc                                                                                                                                          >  8:290351/1‑114 (MQ=255)
                                                tCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCa                                                                                                  <  8:168710/148‑1 (MQ=255)
                                                  aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCTATGAGGTCAACCAGa                                                                                                <  7:274921/148‑1 (MQ=255)
                                                  aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGAt                                                                                               >  6:143706/1‑149 (MQ=255)
                                                   tGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGAc                                                                                                                 >  8:157896/1‑130 (MQ=255)
                                                   tGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGAc                                                                                                                 <  7:157896/130‑1 (MQ=255)
                                                      gTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGACAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTAc                                                                                           >  4:298506/1‑149 (MQ=255)
                                                             ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGt                                                                                    >  1:349268/1‑149 (MQ=255)
                                                             ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGt                                                                                    >  5:74798/1‑149 (MQ=255)
                                                             ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATAACACCAGCGGCGGTGCTGTTTGCGTCGAACGCCAACAGTACATGTTAAGCAAGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGt                                                                                    >  1:285768/1‑149 (MQ=255)
                                                                  tcatcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                             >  4:179025/1‑119 (MQ=255)
                                                                  tcatcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                             <  3:179025/119‑1 (MQ=255)
                                                                    atcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGa                                                                                                           >  3:312903/1‑119 (MQ=255)
                                                                    atcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGa                                                                                                           <  4:312903/119‑1 (MQ=255)
                                                                          ctgctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAc                                                                       <  3:2230/149‑1 (MQ=255)
                                                                                    cTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTg                                                               >  8:75377/1‑147 (MQ=255)
                                                                                    cTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGaa                                                             <  5:143706/149‑1 (MQ=255)
                                                                                     tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGa                                                                                                                  <  1:126017/95‑1 (MQ=255)
                                                                                     tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGa                                                                                                                  >  2:126017/1‑95 (MQ=255)
                                                                                     tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAg                                                            >  2:175976/1‑149 (MQ=255)
                                                                                     tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAg                                                            >  6:100476/1‑149 (MQ=255)
                                                                                               ctgGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCaa                                                  <  6:74798/149‑1 (MQ=255)
                                                                                                  gCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAgg                                                <  3:298506/148‑1 (MQ=255)
                                                                                                         gCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATg                                                                                                            >  2:300895/1‑81 (MQ=255)
                                                                                                         gCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATg                                                                                                            <  1:300895/81‑1 (MQ=255)
                                                                                                          ccAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAgg                                                <  8:20535/140‑1 (MQ=255)
                                                                                                          ccAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAgg                                                >  7:20535/1‑140 (MQ=255)
                                                                                                          ccAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTc                                        >  1:319990/1‑148 (MQ=255)
                                                                                                            aGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGAc                                     <  1:175976/149‑1 (MQ=255)
                                                                                                              gCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGTATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGAc                                     >  4:100213/1‑147 (MQ=255)
                                                                                                                ggTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGt                                 <  7:75377/149‑1 (MQ=255)
                                                                                                                 gTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTc                                >  6:224607/1‑149 (MQ=255)
                                                                                                                 gTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTc                                >  6:229796/1‑149 (MQ=255)
                                                                                                                  tGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTc                                >  3:100270/1‑148 (MQ=255)
                                                                                                                     aTACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTgg                                                                  <  4:151948/111‑1 (MQ=255)
                                                                                                                     aTACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTgg                                                                  >  3:151948/1‑111 (MQ=255)
                                                                                                                     aTACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAgc                            >  7:81588/1‑149 (MQ=255)
                                                                                                                         caccaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGcc                         >  4:323326/1‑148 (MQ=255)
                                                                                                                          accaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGa                       <  5:224607/149‑1 (MQ=255)
                                                                                                                          accaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGa                       <  5:229796/149‑1 (MQ=255)
                                                                                                                                ggcggTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGTATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCaa                 <  3:100213/149‑1 (MQ=255)
                                                                                                                                    gTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgt             <  4:100270/149‑1 (MQ=255)
                                                                                                                                    gTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgt             <  3:323326/149‑1 (MQ=255)
                                                                                                                                     tacTGTTGGTGTTCGACGCCAACAGCACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGgtg                                             <  7:2214/114‑1 (MQ=255)
                                                                                                                                     tacTGTTGGTGTTCGACGCCAACAGCACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGgtg                                             >  8:2214/3‑116 (MQ=255)
                                                                                                                                     tacTGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAgg                                                <  7:182594/111‑1 (MQ=255)
                                                                                                                                     tacTGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAgg                                                >  8:182594/3‑113 (MQ=255)
                                                                                                                                           tGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTAAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGTGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAg       <  6:59895/148‑1 (MQ=255)
                                                                                                                                              cGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCaa                 >  5:170467/1‑135 (MQ=255)
                                                                                                                                              cGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCaa                 <  6:170467/135‑1 (MQ=255)
                                                                                                                                               gTCGAACGCCAACAGTACATGTTCGGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCAt  >  4:343306/1‑149 (MQ=255)
                                                                                                                                               gTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGAt                                                                                               >  8:152618/1‑56 (MQ=255)
                                                                                                                                               gTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGAt                                                                                               <  7:152618/56‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCAT  >  NC_002947/874393‑874684

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