Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,080,148 A→T 5.2% Q31L (CAG→CTG)  PP_5332 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,080,1480AT 5.2% 238.4 / 2.3 76Q31L (CAG→CTG) PP_5332hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base T (2/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AACTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACC  >  NC_002947/6080001‑6080291
                                                                                                                                                     |                                                                                                                                               
aaCTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCcagcag                                                                                                                                                <  8:222072/149‑1 (MQ=255)
   tCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGgc                                                                                                                                             >  1:179733/1‑149 (MQ=255)
         cATCTGCCCGGATCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGgcgatc                                                                                                                                         <  8:297013/147‑4 (MQ=255)
         cATCTGCCCGGATCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGgcgatc                                                                                                                                         >  7:297013/1‑144 (MQ=255)
          aTCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGTGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAgg                                                                                                                                      >  4:9665/1‑149 (MQ=255)
             tGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTcc                                                                                                                                   >  8:86008/1‑149 (MQ=255)
               cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGTGCGCCCTGCTGGGCGATCAgga                                                                                                                                     <  2:179733/145‑2 (MQ=255)
               cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                  >  2:255025/1‑148 (MQ=255)
               cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTc                                                                                                                                 <  7:86008/149‑1 (MQ=255)
               cccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTc                                                                                                                                 <  4:195562/149‑1 (MQ=255)
                ccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                  >  1:244999/1‑147 (MQ=255)
                ccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                  >  7:298230/1‑147 (MQ=255)
                ccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCt                                                                                                                                  >  2:139983/1‑147 (MQ=255)
                ccGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                <  3:99804/149‑1 (MQ=255)
                    gTCCAAGGCCCACT‑GCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTAc                                                                                                                           <  7:82739/149‑1 (MQ=255)
                         aGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCATCAGCTACcagc                                                                                                                       <  6:100254/149‑1 (MQ=255)
                           gCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAAAGAGCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCCAcagca                                                                                                                       >  1:282128/1‑145 (MQ=255)
                            cccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAgg                                                                                                                                               <  7:271951/122‑1 (MQ=255)
                            cccACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAgg                                                                                                                                               >  8:271951/1‑122 (MQ=255)
                                   gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCg                                                                                                              >  2:285562/1‑148 (MQ=255)
                                   gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCg                                                                                                              >  3:131347/1‑148 (MQ=255)
                                   gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCg                                                                                                              <  1:161811/148‑1 (MQ=255)
                                   gCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGACCAAAATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGACCAGCAGGGCGCACAGCTCCTCAGCTACCAGCACGTAGGCg                                                                                                              >  7:282516/1‑148 (MQ=255)
                                     aatcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAc                                                                                                           <  6:29680/149‑1 (MQ=255)
                                     aatcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAc                                                                                                           <  1:255025/149‑1 (MQ=255)
                                       tcaatGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGTGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAAccg                                                                                                         >  6:308789/1‑149 (MQ=255)
                                           tGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACcgccgc                                                                                                     >  3:305656/1‑149 (MQ=255)
                                             agaagaACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctg                                                                                                   >  6:97500/1‑149 (MQ=255)
                                               aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGTGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgct                                                                                                 <  5:308789/149‑1 (MQ=255)
                                               aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgct                                                                                                 >  5:246322/1‑149 (MQ=255)
                                               aagaaCGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCgctgct                                                                                                 >  6:124330/1‑149 (MQ=255)
                                                agaaCGCAATGGCTACCTTTCATGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGctgctg                                                                                                <  6:246322/149‑1 (MQ=255)
                                                   aCGCAATGGCTACCTTTCAGGTCGTAACCGTGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTgg                                                                                             <  3:9665/149‑1 (MQ=255)
                                                    cGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGc                                                                                            <  1:139983/149‑1 (MQ=255)
                                                    cGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGc                                                                                            <  4:305656/149‑1 (MQ=255)
                                                     gCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCt                                                                                           <  8:298230/149‑1 (MQ=255)
                                                       aaTGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGa                                                                                         <  5:124330/149‑1 (MQ=255)
                                                       aaTGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGa                                                                                         <  5:97500/149‑1 (MQ=255)
                                                         tGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCGCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGc                                                                                       >  7:198766/1‑149 (MQ=255)
                                                           gCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGa                                                                                     <  1:285562/149‑1 (MQ=255)
                                                              aCCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGAc                                                                                    >  1:308773/1‑147 (MQ=255)
                                                              aCCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGAc                                                                                    >  1:13879/1‑147 (MQ=255)
                                                                     aGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATc                                                                           >  1:263767/1‑149 (MQ=255)
                                                                     aGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATc                                                                           <  2:263767/149‑1 (MQ=255)
                                                                          gTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                       <  4:131347/148‑1 (MQ=255)
                                                                               cGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAgg                                                                 <  2:13879/149‑1 (MQ=255)
                                                                                     aCATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCGCAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGt                                                           <  8:198766/149‑1 (MQ=255)
                                                                                     aCATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                       <  6:114736/137‑1 (MQ=255)
                                                                                     aCATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTcc                                                                       >  5:114736/1‑137 (MQ=255)
                                                                                                  tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAg                                                                  >  5:72272/1‑131 (MQ=255)
                                                                                                  tCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAg                                                                  <  6:72272/131‑1 (MQ=255)
                                                                                                    aaCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCgg                                              <  2:164705/149‑1 (MQ=255)
                                                                                                     aCTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCggg                                             >  5:282181/1‑149 (MQ=255)
                                                                                                      ctgctgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCgggg                                            <  6:282181/149‑1 (MQ=255)
                                                                                                          tgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGcc                                                                              <  5:23826/111‑1 (MQ=255)
                                                                                                          tgGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGcc                                                                              >  6:23826/1‑111 (MQ=255)
                                                                                                            gCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAggg                                                                <  3:41024/123‑1 (MQ=255)
                                                                                                            gCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAggg                                                                >  4:41024/1‑123 (MQ=255)
                                                                                                            gCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCgg                                      >  7:208120/1‑149 (MQ=255)
                                                                                                                  cagcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTAc                                                     <  6:331022/128‑1 (MQ=255)
                                                                                                                  cagcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTAc                                                     >  5:331022/1‑128 (MQ=255)
                                                                                                                   agcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTAc                                                     <  6:73057/127‑1 (MQ=255)
                                                                                                                   agcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGCGCAAGTCGGTAc                                                     >  5:73057/1‑127 (MQ=255)
                                                                                                                    gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                >  4:310728/1‑51 (MQ=255)
                                                                                                                    gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCa                                                                                                                                <  3:310728/51‑1 (MQ=255)
                                                                                                                    gcagTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGAc                                                                                    <  4:58223/95‑1 (MQ=255)
                                                                                                                            cGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAgg                                                                                                                >  6:878/1‑59 (MQ=255)
                                                                                                                              tGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                    <  8:208120/149‑1 (MQ=255)
                                                                                                                              tGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                    <  2:33602/149‑1 (MQ=255)
                                                                                                                                   aaCTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAAc               <  2:308773/149‑1 (MQ=255)
                                                                                                                                          gATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAgc        >  3:27682/1‑149 (MQ=255)
                                                                                                                                          gATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAgc        >  8:98512/1‑149 (MQ=255)
                                                                                                                                           aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                    <  5:148949/136‑1 (MQ=255)
                                                                                                                                           aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTg                    >  6:148949/1‑136 (MQ=255)
                                                                                                                                           aTCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAg         <  2:244999/147‑1 (MQ=255)
                                                                                                                                              gCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCaa    >  5:324097/1‑149 (MQ=255)
                                                                                                                                                ccAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAAcc  >  2:277730/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AACTCATTGCATCTGCCCGGGTCCAAGGCCCACTGGCAATCAATGAGAAGAACGCAATGGCTACCTTTCAGGTCGTAACCGAGCAACATGGCGAA‑‑CTCAACTGCTGGCGTATCAGCAGTGACCGTGCCGAACTCCTGATCGCCCAGCAGGGCGCACAGGTCCTCAGCTACCAGCAGGTAGGCGAACCGCCGCTGCTGTGGCTGAGCGACCAGGCCATCTTCCGCCAGGGCAAGTCGGTACGCGCCGGGGTACCGGTGTGCTGGCCATGGTTTGGCAACCTGCAGCGCAACC  >  NC_002947/6080001‑6080291

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑