Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,080,155 A→T 6.3% A33A (GCA→GCT) ‡ PP_5332 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,080,1550AT 6.3% 270.6 / 4.0 90A33A (GCA→GCT) ‡PP_5332hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base T (2/4);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGCCGTC  >  NC_002947/6080009‑6080302
                                                                                                                                                  |                                                                                                                                                    
gCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGca                                                                                                                                                      >  4:160790/1‑147 (MQ=255)
      gCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                             >  1:385143/1‑149 (MQ=255)
       cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCTAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTc                                                                                                                                               <  4:282234/146‑1 (MQ=255)
       cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGGCGATCAgga                                                                                                                                                <  1:152624/145‑2 (MQ=255)
       cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCTGCTGGGCGATCAgga                                                                                                                                                <  7:279326/145‑2 (MQ=255)
         cGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAg                                                                                                                                          >  6:358898/1‑149 (MQ=255)
            gTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTa                                                                                                                                       <  6:291475/149‑1 (MQ=255)
                aaGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGcac                                                                                                                                                    >  3:259233/1‑132 (MQ=255)
                aaGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGcac                                                                                                                                                    <  4:259233/132‑1 (MQ=255)
                    cccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAgg                                                                                                                               <  4:202075/149‑1 (MQ=255)
                     ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCCCCGGTCCTCAGCTACCAGCAgg                                                                                                                               >  7:341731/1‑148 (MQ=255)
                     ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAg                                                                                                                                          >  2:180412/1‑137 (MQ=255)
                     ccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAg                                                                                                                                          <  1:180412/137‑1 (MQ=255)
                       aCTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAg                                                                                                                            <  8:236565/149‑1 (MQ=255)
                           gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                             >  2:118739/1‑128 (MQ=255)
                           gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                             <  1:118739/128‑1 (MQ=255)
                           gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGc                                                                                                                          <  6:196119/147‑1 (MQ=255)
                             aatcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAc                                                                                                                      >  1:261050/1‑149 (MQ=255)
                              atcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCTAACTCAAATGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCACCCAGCAGGGCGCACAGGTCCTCATCTACCAGCAGGTAGGCGaaac                                                                                                                     >  8:184031/1‑147 (MQ=255)
                              atcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAcc                                                                                                                     >  3:313337/1‑149 (MQ=255)
                                   tGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACcgccgc                                                                                                                >  1:271774/1‑149 (MQ=255)
                                   tGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACcgccgc                                                                                                                >  1:217305/1‑149 (MQ=255)
                                     agaagaACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctg                                                                                                              >  1:183498/1‑149 (MQ=255)
                                      gaagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgc                                                                                                             >  7:327370/1‑149 (MQ=255)
                                       aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgct                                                                                                            >  3:279399/1‑149 (MQ=255)
                                         gaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCTAACTCCTGATCGCCCATCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCtgt                                                                                                          <  6:184707/149‑1 (MQ=255)
                                         gaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCtgt                                                                                                          >  5:189318/1‑149 (MQ=255)
                                          aaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCtgtg                                                                                                         <  4:358397/149‑1 (MQ=255)
                                             gCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCt                                                                                                      <  5:358898/149‑1 (MQ=255)
                                                 tGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGc                                                                                                  <  2:217305/149‑1 (MQ=255)
                                                     tACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGAc                                                                                               >  1:177786/1‑148 (MQ=255)
                                                       ccTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAg                                                                                            <  5:318359/149‑1 (MQ=255)
                                                       ccTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAg                                                                                            <  2:261050/149‑1 (MQ=255)
                                                         tttCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGCACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGc                                                                                          >  7:171797/1‑149 (MQ=255)
                                                           tCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGc                                                                                          <  7:339009/147‑1 (MQ=255)
                                                             aGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCAt                                                                                       <  6:189318/148‑1 (MQ=255)
                                                                  gTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                                  <  4:313337/148‑1 (MQ=255)
                                                                  gTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                                  <  5:12312/148‑1 (MQ=255)
                                                                   tAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGc                                                                                <  3:160790/149‑1 (MQ=255)
                                                                   tAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGc                                                                                >  7:332858/1‑149 (MQ=255)
                                                                    aaCCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                           <  6:140939/89‑1 (MQ=255)
                                                                    aaCCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                           >  5:140939/1‑89 (MQ=255)
                                                                    aaCCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGc                                                                                >  6:101284/1‑148 (MQ=255)
                                                                        gAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAggg                                                                           <  5:101284/149‑1 (MQ=255)
                                                                           cAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCaa                                                                        <  4:279399/149‑1 (MQ=255)
                                                                             aCATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGt                                                                      <  7:184031/149‑1 (MQ=255)
                                                                             aCATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGt                                                                      <  4:356761/149‑1 (MQ=255)
                                                                                tGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCgg                                                                   >  5:242269/1‑149 (MQ=255)
                                                                                  gCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTa                                                                 <  7:326129/149‑1 (MQ=255)
                                                                                     aaCTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTAcgc                                                              <  8:332858/149‑1 (MQ=255)
                                                                                       cTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACgcgc                                                            >  1:139725/1‑149 (MQ=255)
                                                                                        tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTa                                                                 >  3:220568/1‑143 (MQ=255)
                                                                                        tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTa                                                                 <  4:220568/143‑1 (MQ=255)
                                                                                           aCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCggg                                                        <  6:114920/149‑1 (MQ=255)
                                                                                             tgctgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCgggg                                                       >  5:208680/1‑148 (MQ=255)
                                                                                               ctgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTAc                                                    <  2:220665/149‑1 (MQ=255)
                                                                                                       tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCGGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                            >  2:63361/1‑149 (MQ=255)
                                                                                                       tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                            <  8:327370/149‑1 (MQ=255)
                                                                                                       tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                            >  2:117579/1‑149 (MQ=255)
                                                                                                       tCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGc                                            >  1:308313/1‑149 (MQ=255)
                                                                                                          gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTgg                                         >  5:157431/1‑149 (MQ=255)
                                                                                                            agTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGc                                        >  7:260299/1‑148 (MQ=255)
                                                                                                             gTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCa                                      >  5:155226/1‑149 (MQ=255)
                                                                                                              tGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCAt                                     >  4:22037/1‑149 (MQ=255)
                                                                                                                aCCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATg                                    <  4:253747/148‑1 (MQ=255)
                                                                                                                   gTGCCGAACTCCTGATCGCCCAGCAGGGCGCAGCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGtt                                 >  1:242935/1‑149 (MQ=255)
                                                                                                                   gTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                                <  6:157431/149‑1 (MQ=255)
                                                                                                                    tGCCGAAGTCCTGATCGCCCAGCAGGGCGCTCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                               <  2:373112/149‑1 (MQ=255)
                                                                                                                    tGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCTGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                               <  1:63361/149‑1 (MQ=255)
                                                                                                                        gAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCaa                           <  6:242269/149‑1 (MQ=255)
                                                                                                                         aaCTCCTGATCGCCCCGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                          <  8:341731/149‑1 (MQ=255)
                                                                                                                         aaCTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                          <  2:271774/149‑1 (MQ=255)
                                                                                                                           cTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCt                        >  7:287289/1‑149 (MQ=255)
                                                                                                                              ctgtgcGCCCTGCTGGGCGATCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTg                       >  1:212198/6‑147 (MQ=255)
                                                                                                                              cTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTg                       <  2:183498/147‑1 (MQ=255)
                                                                                                                              cTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTg                       <  6:208680/147‑1 (MQ=255)
                                                                                                                                 aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCg                                                  >  4:44831/1‑117 (MQ=255)
                                                                                                                                 aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCg                                                  <  3:44831/117‑1 (MQ=255)
                                                                                                                                    gCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCaa               <  2:308313/149‑1 (MQ=255)
                                                                                                                                     cccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                          <  4:50670/137‑1 (MQ=255)
                                                                                                                                     cccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc                          >  3:50670/1‑137 (MQ=255)
                                                                                                                                     cccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAc              >  5:72187/1‑149 (MQ=255)
                                                                                                                                      ccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc             >  4:325151/1‑149 (MQ=255)
                                                                                                                                      ccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc             <  4:85675/149‑1 (MQ=255)
                                                                                                                                      ccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc             >  3:85675/1‑149 (MQ=255)
                                                                                                                                        tgctgggCGATCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                                <  2:326625/124‑1 (MQ=255)
                                                                                                                                        tgctgggCGATCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGttt                                >  1:326625/5‑128 (MQ=255)
                                                                                                                                        agcagGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCt           <  6:155226/149‑1 (MQ=255)
                                                                                                                                           agGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAg        >  7:175614/1‑149 (MQ=255)
                                                                                                                                            gggCGCACACGTCCTCAGCTACCAGCAGGTAGGCGAACCGCAGCTGATGTGGCTGCGCTACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGATGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCCgg       >  2:239498/1‑149 (MQ=255)
                                                                                                                                             ggCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAgg       <  2:216805/148‑1 (MQ=255)
                                                                                                                                             ggCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGc      >  7:225163/1‑149 (MQ=255)
                                                                                                                                             ggCGCACAGGTCCTCAGCTACCAGCAGGAAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAgg       >  1:216805/1‑148 (MQ=255)
                                                                                                                                                  acaGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGc                                                                                                       >  1:280311/1‑47 (MQ=255)
                                                                                                                                                  acaGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGc                                                                                                       <  2:280311/47‑1 (MQ=255)
                                                                                                                                                  acaGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGCCGTc  >  4:97666/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAACTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACCCTCAGGCCGTC  >  NC_002947/6080009‑6080302

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: