Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 516,279 A→T 8.7% 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.7% 202.3 / 5.4 68L238Q (CTG→CAG) PP_0429short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base T (4/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGC  >  NC_002947/516141‑516426
                                                                                                                                          |                                                                                                                                                   
tCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTCGCTATTCCATGCCTGAGGCGCAATAGCGGCGGAGGATGGCCCCGGCCAGCTTGCGCACAAGCCAGCGCGGCTGCAGCCGAgg                                                                                                                                           >  5:171573/1‑149 (MQ=255)
  cacaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTTCGCACCAGCCCCCCCGGCTGCGGCcccgcgg                                                                                                                                         >  5:176493/1‑146 (MQ=255)
  cacaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                          <  6:130327/148‑1 (MQ=255)
  cacaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCcgcgct                                                                                                                                          <  6:282538/148‑3 (MQ=255)
   acaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgca                                                                                                                                                     >  2:101985/1‑136 (MQ=255)
   acaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgca                                                                                                                                                     <  1:101985/136‑1 (MQ=255)
   acaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                          >  6:178031/1‑147 (MQ=255)
   acaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                          <  5:178031/147‑1 (MQ=255)
      gCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                     >  6:275596/1‑149 (MQ=255)
        aCTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACg                                                                                                                                   >  3:123202/1‑149 (MQ=255)
          ttGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCTGCTGCCGCGCTg                                                                                                                                         <  3:193808/141‑1 (MQ=255)
          ttGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCTGCTGCCGCGCTg                                                                                                                                         >  4:193808/1‑141 (MQ=255)
            gCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgccag                                                                                                                               >  5:41427/1‑149 (MQ=255)
            gCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgccag                                                                                                                               >  4:86918/1‑149 (MQ=255)
                         ggTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgcgc                                                                                                                  >  8:239943/1‑149 (MQ=255)
                         ggTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgcgc                                                                                                                  >  8:75014/1‑149 (MQ=255)
                               gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCgccggccg                                                                                                            >  7:236908/1‑149 (MQ=255)
                                  gAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCGTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGc                                                                                                         <  4:203243/149‑1 (MQ=255)
                                       tACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGAt                                                                                                    <  3:15820/149‑1 (MQ=255)
                                         cTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATc                                                                                                   <  6:41427/148‑1 (MQ=255)
                                         cTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCa                                                                                                  <  6:172571/149‑1 (MQ=255)
                                          tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATcgg                                                                                                 >  5:267376/1‑147 (MQ=255)
                                          tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGCCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAg                                                                                                 >  6:299264/1‑149 (MQ=255)
                                              ccGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTgcg                                                                                             >  5:276221/1‑149 (MQ=255)
                                               cGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTgcgc                                                                                            <  5:275596/149‑1 (MQ=255)
                                                           tccatgTGGCTATTGCATGCCTGTGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGGGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCa                                                                                <  6:171573/144‑1 (MQ=255)
                                                           gCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                  <  3:101343/147‑1 (MQ=255)
                                                           gCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                  >  5:268545/1‑147 (MQ=255)
                                                           gCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                                  <  3:67797/147‑1 (MQ=255)
                                                              cAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTgcg                                                                                             <  8:174577/133‑1 (MQ=255)
                                                              cAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTgcg                                                                                             >  7:174577/1‑133 (MQ=255)
                                                                gTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGcc                                                                           >  1:281463/1‑149 (MQ=255)
                                                                     tATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCAc                                                                      >  6:312023/1‑149 (MQ=255)
                                                                       ttCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCAc                                                                      <  6:276221/147‑1 (MQ=255)
                                                                         ccATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgt                                                                   <  6:267376/148‑1 (MQ=255)
                                                                          cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgtg                                                                  >  4:160622/1‑148 (MQ=255)
                                                                           aTGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTa                                                                >  3:74719/1‑149 (MQ=255)
                                                                             gCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTa                                                                >  2:247842/1‑147 (MQ=255)
                                                                              ccTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTaga                                                              <  6:268545/148‑3 (MQ=255)
                                                                               cTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgc                                                                                                                    >  2:208160/1‑93 (MQ=255)
                                                                               cTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgc                                                                                                                    <  1:208160/93‑1 (MQ=255)
                                                                                 gAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                          >  6:99515/1‑149 (MQ=255)
                                                                                      gCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCCCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGACATGATCAGTGCGAGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACtt                                                     >  7:21836/1‑149 (MQ=255)
                                                                                       cAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGt                                                                              >  2:146833/1‑123 (MQ=255)
                                                                                       cAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGt                                                                              <  1:146833/123‑1 (MQ=255)
                                                                                        aaTAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTcc                                                   >  8:83599/1‑149 (MQ=255)
                                                                                              ggcggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCgggg                                             <  7:239943/149‑1 (MQ=255)
                                                                                                 ggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCTGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTg                                          <  5:299264/149‑1 (MQ=255)
                                                                                                      aTGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAgg                                       >  3:231665/1‑147 (MQ=255)
                                                                                                          ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAgcgc                                 <  2:281463/149‑1 (MQ=255)
                                                                                                          ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAg                                    <  5:186731/146‑1 (MQ=255)
                                                                                                          ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAg                                    >  6:186731/1‑146 (MQ=255)
                                                                                                               gCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                            >  3:90551/1‑149 (MQ=255)
                                                                                                                  aGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGcc                         <  4:123202/149‑1 (MQ=255)
                                                                                                                  aGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGCCATGATCAGTGCGCGGTGCTTGTCCAGAGCCCGCACGGTCTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTAGc                          >  3:221120/1‑148 (MQ=255)
                                                                                                                    cTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                            <  1:133460/144‑1 (MQ=255)
                                                                                                                    cTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                            >  2:133460/1‑144 (MQ=255)
                                                                                                                    cTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCAt                       <  5:312023/149‑1 (MQ=255)
                                                                                                                     ttGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATc                      >  2:295899/1‑149 (MQ=255)
                                                                                                                        cgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGa                    <  4:90551/148‑1 (MQ=255)
                                                                                                                         gcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCAACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGc                  >  8:268993/1‑149 (MQ=255)
                                                                                                                          cACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGcc                 >  3:257139/1‑149 (MQ=255)
                                                                                                                            cctgccagcGGTGCAGCAGCCTCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCTGGCATGTTCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGTCTTAGGCAGCGCGGCTTGCCATCGAGCcgc               <  8:218699/146‑1 (MQ=255)
                                                                                                                             cagccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgcg              >  3:121049/1‑149 (MQ=255)
                                                                                                                               gccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCgcgc             >  5:95993/1‑148 (MQ=255)
                                                                                                                                ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCaa           <  3:86918/149‑1 (MQ=255)
                                                                                                                                 cagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAAt          <  7:75014/149‑1 (MQ=255)
                                                                                                                                 caccgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAAt          >  6:267648/4‑149 (MQ=255)
                                                                                                                                     gcgGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATgcgc      >  6:93285/1‑149 (MQ=255)
                                                                                                                                       ggCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCg    <  4:74719/149‑1 (MQ=255)
                                                                                                                                         cagcagCCGCGGTGCGAACGCCAGCCAGCGGTTGCGccggccgg                                                                                                           >  2:81780/1‑44 (MQ=255)
                                                                                                                                         cagcagCCGCGGTGCGAACGCCAGCCAGCGGTTGCGccggccgg                                                                                                           <  1:81780/44‑1 (MQ=255)
                                                                                                                                         cagcagCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGc  <  4:257139/149‑1 (MQ=255)
                                                                                                                                         cagcagCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGc  <  5:93285/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCGCCGGC  >  NC_002947/516141‑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: