Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,646,266 A→T 5.6% D540V (GAC→GTC)  PP_4954 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,646,2660AT 5.6% 218.4 / 2.7 71D540V (GAC→GTC) PP_4954membrane protein
Reads supporting (aligned to +/- strand):  ref base A (40/27);  new base T (2/2);  total (42/29)
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.85e-01

ACCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGC  >  NC_002947/5646121‑5646410
                                                                                                                                                 |                                                                                                                                                
acCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCCGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTg                                                                                                                                               >  5:69370/1‑149 (MQ=255)
acCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCCGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTg                                                                                                                                               >  3:29154/1‑149 (MQ=255)
 cccTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                              <  7:249206/149‑1 (MQ=255)
  ccTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                              >  1:68833/1‑148 (MQ=255)
    tGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAg                                                                                                                                           >  2:79595/1‑149 (MQ=255)
      gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAggg                                                                                                                                         <  1:79595/149‑1 (MQ=255)
      gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCAGTcaggcagg                                                                                                                                          <  8:19118/148‑1 (MQ=255)
        cACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGcc                                                                                                                                       >  5:75250/1‑149 (MQ=255)
        cACTGCCGCGATCACCCTGAAGGTGCTGTTCTGCTTCAACCAGATCGGAGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCAGTcaggcaggaac                                                                                                                                       <  2:299952/149‑4 (MQ=255)
               gcgATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCt                                                                                                                                >  5:87011/1‑149 (MQ=255)
               gcgATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCt                                                                                                                                <  6:75250/149‑1 (MQ=255)
               gcgATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCt                                                                                                                                <  8:13741/149‑1 (MQ=255)
                     aCCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAAc                                                                                                                            >  3:257690/1‑147 (MQ=255)
                                         cTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTgg                                                                                                      >  7:194312/1‑149 (MQ=255)
                                           tCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGcc                                                                                                    >  7:111326/1‑149 (MQ=255)
                                           tCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGcc                                                                                                    >  8:154108/1‑149 (MQ=255)
                                               ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTgcg                                                                                                >  2:265884/1‑149 (MQ=255)
                                               ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTgcg                                                                                                >  1:312203/1‑149 (MQ=255)
                                                        cGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCAAATACCCTGGCCTGCGCCAGCgtg                                                                                        >  7:102800/1‑148 (MQ=255)
                                                          aTGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCAAATACCCTGGCCTGCGCCAGCGTGTAc                                                                                     <  8:102800/149‑1 (MQ=255)
                                                                        ttCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCAt                                                                       <  1:265884/149‑1 (MQ=255)
                                                                          ccTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      <  3:27173/148‑1 (MQ=255)
                                                                           cTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      >  2:340409/1‑147 (MQ=255)
                                                                              ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATgcagc                                                                  <  6:87011/148‑1 (MQ=255)
                                                                               cgcgCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGACGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGcagcag                                                                >  4:166708/1‑149 (MQ=255)
                                                                                     tGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGcc                                                          >  7:163839/1‑149 (MQ=255)
                                                                                           aTACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGt                                                               >  8:13083/1‑138 (MQ=255)
                                                                                           aTACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGt                                                               <  7:13083/138‑1 (MQ=255)
                                                                                            tACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCgtgt                                                                                       >  8:113544/1‑113 (MQ=255)
                                                                                            tACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCgtgt                                                                                       <  7:113544/113‑1 (MQ=255)
                                                                                              cccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAg                                                 >  7:248456/1‑149 (MQ=255)
                                                                                               ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGTAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                                >  5:248089/1‑149 (MQ=255)
                                                                                               ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                                >  4:104470/1‑149 (MQ=255)
                                                                                               ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                                >  2:175239/1‑149 (MQ=255)
                                                                                               ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCAGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                                <  2:312203/149‑1 (MQ=255)
                                                                                                 tGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGc                                                           >  8:212573/1‑136 (MQ=255)
                                                                                                 tGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGc                                                           <  7:212573/136‑1 (MQ=255)
                                                                                                      ggCATCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGcgacga                                         <  4:29154/149‑1 (MQ=255)
                                                                                                      ggCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGcgacga                                         <  6:69370/149‑1 (MQ=255)
                                                                                                        cAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGAcc                                       >  1:193065/1‑149 (MQ=255)
                                                                                                        cAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGAcc                                       >  1:333021/1‑149 (MQ=255)
                                                                                                          gCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                     <  1:340409/149‑1 (MQ=255)
                                                                                                           ccTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTgg                                    <  8:124900/149‑1 (MQ=255)
                                                                                                            cTGATTGCCATTCTTGCAGTGTTCCTGATCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGc                                   >  5:83067/1‑149 (MQ=255)
                                                                                                                 tGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTAcc                              >  4:77302/1‑149 (MQ=255)
                                                                                                                  ttgaTTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTgg                                    <  8:57308/139‑1 (MQ=255)
                                                                                                                  gCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGACCACGGCAAGCGCGACGACCTGGcagatcg                             >  7:57308/1‑143 (MQ=255)
                                                                                                                   ccATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTa                                                                                      <  2:234839/91‑1 (MQ=255)
                                                                                                                   ccATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTa                                                                                      >  1:234839/1‑91 (MQ=255)
                                                                                                                   ccATTATTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCAAATACCCTCGCCTGCGCCAGCGTGTACCTGCGCCAGATCAAGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGc                            >  3:30132/1‑149 (MQ=255)
                                                                                                                    cATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTgg                                    <  2:221079/140‑1 (MQ=255)
                                                                                                                    cATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTgg                                    >  1:221079/1‑140 (MQ=255)
                                                                                                                    cATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                     >  4:108884/1‑139 (MQ=255)
                                                                                                                    cATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                     <  3:108884/139‑1 (MQ=255)
                                                                                                                     aTTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGTGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCt                          >  3:281390/1‑149 (MQ=255)
                                                                                                                     aTTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCt                          <  8:194312/149‑1 (MQ=255)
                                                                                                                      ttCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTg                         <  8:248456/149‑1 (MQ=255)
                                                                                                                        cTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGc                       >  6:100668/1‑149 (MQ=255)
                                                                                                                        cTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGc                       >  6:95791/1‑149 (MQ=255)
                                                                                                                             aGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgc                  >  8:301123/1‑149 (MQ=255)
                                                                                                                             aGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgc                  <  8:111326/149‑1 (MQ=255)
                                                                                                                                   cctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCGGGCCAGgcggtag             <  1:175239/148‑5 (MQ=255)
                                                                                                                                    ctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGc           >  6:251024/1‑149 (MQ=255)
                                                                                                                                     tgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCGAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGcc          >  7:221434/1‑149 (MQ=255)
                                                                                                                                       gccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca        >  5:212155/3‑149 (MQ=255)
                                                                                                                                        ccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca        >  8:204251/2‑148 (MQ=255)
                                                                                                                                           tgcctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                     >  8:208302/1‑116 (MQ=255)
                                                                                                                                           tgcctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTg                                     <  7:208302/116‑1 (MQ=255)
                                                                                                                                             cctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACg   <  7:301123/148‑1 (MQ=255)
                                                                                                                                             cctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGc  <  2:193065/149‑1 (MQ=255)
                                                                                                                                             cctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGc  >  6:61688/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGC  >  NC_002947/5646121‑5646410

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