Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 516,279 A→T 8.1% L238Q (CTG→CAG)  PP_0429 ← short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947516,2790AT 8.1% 214.9 / 7.1 74L238Q (CTG→CAG) PP_0429short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (38/30);  new base T (4/2);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

ATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGC  >  NC_002947/516133‑516426
                                                                                                                                                  |                                                                                                                                                   
aTAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgcag                                                                                                                                                    >  7:40665/1‑148 (MQ=255)
       gTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCcgcg                                                                                                                                            <  7:42730/149‑1 (MQ=255)
         ccACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                          >  2:83504/1‑149 (MQ=255)
          cacaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTg                                                                                                                                         >  2:212725/1‑149 (MQ=255)
           acaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGc                                                                                                                                        >  1:34498/1‑149 (MQ=255)
            caGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCg                                                                                                                                       <  2:34498/149‑1 (MQ=255)
               cACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCTGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGc                                                                                                                                        >  5:277769/1‑145 (MQ=255)
               cACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGc                                                                                                                                        <  6:277769/145‑1 (MQ=255)
                     cTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgccagc                                                                                                                              >  7:271125/1‑149 (MQ=255)
                            gTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGt                                                                                                                       >  6:240401/1‑149 (MQ=255)
                                 ggTCGTGCCGAGTGTACTCAGCCGGGAGCCGCTCGCCCAGTGTCTATTCCAAGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgcgc                                                                                                                  <  4:232905/149‑1 (MQ=255)
                                       gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgcca                                                                                                                                <  8:273101/129‑1 (MQ=255)
                                       gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgcca                                                                                                                                >  7:273101/1‑129 (MQ=255)
                                       gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCgccggccg                                                                                                            >  5:127045/1‑149 (MQ=255)
                                        ccgcgctgcctcagccccgaggaagtggcgctgtacaccgtgcgggcactggacaagaaccgcgcactgatcatgCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGccggccgg                                                                                                           >  3:130087‑M1/76‑149 (MQ=255)
                                          gAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCggg                                                                                                          >  7:141240/1‑148 (MQ=255)
                                          gAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCggg                                                                                                          >  2:52508/1‑148 (MQ=255)
                                           aGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCCGCGCTggctgg                                                                                                                                    >  1:290803/1‑116 (MQ=255)
                                           aGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCCGCGCTggctgg                                                                                                                                    <  2:290803/121‑6 (MQ=255)
                                            gctgcctcagccccgaggaagtggcgctgtacaccgtgcgggcactggacaagaaccgcgcactgatcatgCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCAt                                                                                                       >  8:289306‑M1/72‑149 (MQ=255)
                                               tACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGAt                                                                                                    <  8:135406/149‑1 (MQ=255)
                                                cctcagccccgaggaagtggcgctgtacaccgtgcgggcactggacaagaaccgcgcactgatcatgCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCCGCGGTTGCGCCGGCCGGGCATGAt                                                                                                    >  1:146444‑M1/68‑148 (MQ=255)
                                                 cTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCa                                                                                                  >  3:205055/1‑149 (MQ=255)
                                                  tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAg                                                                                                 >  7:245074/1‑149 (MQ=255)
                                                  tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAg                                                                                                 >  7:181941/1‑149 (MQ=255)
                                                  tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAg                                                                                                 >  8:35895/1‑149 (MQ=255)
                                                             ccGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCt                                                                                      >  8:278845/1‑149 (MQ=255)
                                                               gCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTg                                                                                    <  1:83504/149‑1 (MQ=255)
                                                                   gCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                  <  4:205055/147‑1 (MQ=255)
                                                                   gCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                  <  5:240401/147‑1 (MQ=255)
                                                                   gCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                  >  4:54398/1‑147 (MQ=255)
                                                                           gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCg                                                                         >  2:226128/1‑148 (MQ=255)
                                                                            cTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCa                                                                       >  8:54963/1‑149 (MQ=255)
                                                                            cTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCa                                                                       >  4:19681/1‑149 (MQ=255)
                                                                               ttCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCAc                                                                      <  8:271125/147‑1 (MQ=255)
                                                                               ttCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCAc                                                                      <  1:183961/147‑1 (MQ=255)
                                                                                tCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGccc                                                                          >  2:155884/1‑142 (MQ=255)
                                                                                tCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGccc                                                                          <  1:155884/142‑1 (MQ=255)
                                                                                 ccATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCCGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgt                                                                   <  1:226128/148‑1 (MQ=255)
                                                                                  cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgtg                                                                  >  3:305727/1‑148 (MQ=255)
                                                                                  cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgtg                                                                  <  8:141240/148‑1 (MQ=255)
                                                                                     gCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCa                                                                                                  >  5:285719/1‑113 (MQ=255)
                                                                                     gCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCa                                                                                                  <  6:285719/113‑1 (MQ=255)
                                                                                      ccTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAg                                                             <  7:54963/149‑1 (MQ=255)
                                                                                          aGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcg                                                           <  8:245074/147‑1 (MQ=255)
                                                                                                 aTAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCt                                                  >  7:45327/1‑149 (MQ=255)
                                                                                                 aTAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCt                                                  <  1:212725/149‑1 (MQ=255)
                                                                                                    gcggcgGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCg                                                <  6:289427/148‑1 (MQ=255)
                                                                                                        cggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCt                                           >  4:101025/1‑149 (MQ=255)
                                                                                                             gATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGa                                         <  2:216279/146‑1 (MQ=255)
                                                                                                             gATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGa                                         >  1:216279/1‑146 (MQ=255)
                                                                                                              aTGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCa                                     <  1:52508/149‑1 (MQ=255)
                                                                                                                 gCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCgg                                               >  6:270095/1‑136 (MQ=255)
                                                                                                                  ccccGGCCAGGTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAgcgc                                 <  6:302071/149‑1 (MQ=255)
                                                                                                                  ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGa                                         >  3:248369/1‑141 (MQ=255)
                                                                                                                  ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGa                                         <  4:248369/141‑1 (MQ=255)
                                                                                                                    CCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGTGCATGATCAGTGCGCGGTTCTTGTCCAGt                                                                              >  1:78846‑M1/1‑102 (MQ=255)
                                                                                                                    CCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGCGCATGATCAGTGCGCGGTTCTTGTCCAGt                                                                              <  2:78846‑M1/102‑1 (MQ=255)
                                                                                                                        CCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTg                           >  3:54398‑M1/1‑149 (MQ=255)
                                                                                                                          aGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGc                          >  5:247130/1‑148 (MQ=255)
                                                                                                                               gcgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATc                      >  6:278602/1‑147 (MQ=255)
                                                                                                                               gcgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATc                      >  2:146521/1‑147 (MQ=255)
                                                                                                                                  cACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGcc                 <  8:181941/149‑1 (MQ=255)
                                                                                                                                   aCCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcg                >  1:116810/1‑149 (MQ=255)
                                                                                                                                     cagccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgcg              <  8:45327/149‑1 (MQ=255)
                                                                                                                                      agccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgcg              >  8:82449/1‑148 (MQ=255)
                                                                                                                                        ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAgg                                       <  7:118367/121‑1 (MQ=255)
                                                                                                                                        ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAgg                                       >  8:118367/1‑121 (MQ=255)
                                                                                                                                         cagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAAt          <  6:305455/149‑1 (MQ=255)
                                                                                                                                         cagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGGGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAAt          <  6:310029/149‑1 (MQ=255)
                                                                                                                                         caccgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                            >  6:212847/4‑131 (MQ=255)
                                                                                                                                         caccgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                            <  5:212847/128‑1 (MQ=255)
                                                                                                                                           gcgcGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATgc        <  4:305727/149‑1 (MQ=255)
                                                                                                                                           ccgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATgc        >  7:278830/2‑149 (MQ=255)
                                                                                                                                               ggCTGCTGCCGCGCTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCg    >  1:290042/1‑149 (MQ=255)
                                                                                                                                                 cagcagCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGc  <  6:247130/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGC  >  NC_002947/516133‑516426

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