Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,779,044 C→A 5.1% Q177H (CAG→CAT)  PP_5065 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,779,0440CA 5.1% 259.4 / 2.3 78Q177H (CAG→CAT) PP_5065hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (2/2);  total (40/38)
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.78e-01

TCAGCTTGTGCAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAACTGCCAGAACGC  >  NC_002947/5778896‑5779190
                                                                                                                                                    |                                                                                                                                                  
tCAGCTTGTGCAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTcc                                                                                                                                                    >  7:248328/1‑149 (MQ=255)
   gCTTGTGCAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTgg                                                                                                                                                 <  8:248328/149‑1 (MQ=255)
        tgCAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTc                                                                                                                                            <  6:113975/149‑1 (MQ=255)
         gCAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCAGGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCt                                                                                                                                             >  4:69158/1‑147 (MQ=255)
          cAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGTCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCAGGAACTCATGgcgg                                                                                                                                              <  7:166775/145‑2 (MQ=255)
          cAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAg                                                                                                                                          >  4:193838/1‑149 (MQ=255)
          cAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCAGGAACTCATGgcgg                                                                                                                                              >  8:166775/1‑144 (MQ=255)
              aaCATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTg                                                                                                                                      >  3:81916/1‑149 (MQ=255)
                cATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGc                                                                                                                                    >  1:189475/1‑149 (MQ=255)
                         ggCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAgcgc                                                                                                                             <  5:110511/147‑1 (MQ=255)
                          gcgcTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAg                                                                                                                           <  6:278776/148‑1 (MQ=255)
                          gcgcTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAg                                                                                                                           <  6:273785/148‑1 (MQ=255)
                          gcgcTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGa                                                                                                                          <  4:221321/149‑1 (MQ=255)
                                aTACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGcc                                                                                                                    <  8:100595/149‑1 (MQ=255)
                                  aCCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGcccg                                                                                                                  >  2:297760/1‑148 (MQ=255)
                                  aCCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCca                                                                                                                  >  1:24489/1‑149 (MQ=255)
                                  aCCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCca                                                                                                                  >  2:17924/1‑149 (MQ=255)
                                              aaGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTgg                                                                                                                      <  8:145504/133‑1 (MQ=255)
                                              aaGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTgg                                                                                                                      >  7:145504/1‑133 (MQ=255)
                                              aaGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGc                                                                                                                     <  2:85378/134‑1 (MQ=255)
                                              aaGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGc                                                                                                                     >  1:85378/1‑134 (MQ=255)
                                              aaGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCACAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATg                                                                                                       >  3:297849/1‑148 (MQ=255)
                                               agagCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCTACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGaa                                                                                                     <  1:86557/149‑1 (MQ=255)
                                                 agCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGc                                                                                                   <  3:177353/149‑1 (MQ=255)
                                                  gCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCg                                                                                                  <  3:304581/149‑1 (MQ=255)
                                                  gCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCg                                                                                                  <  1:17924/149‑1 (MQ=255)
                                                  gCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCg                                                                                                  >  2:67595/1‑149 (MQ=255)
                                                             aGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCa                                                                                       <  1:67595/149‑1 (MQ=255)
                                                              ggTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAg                                                                                      <  1:175458/149‑1 (MQ=255)
                                                                ttGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGa                                                                                    >  5:47075/1‑149 (MQ=255)
                                                                ttGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGa                                                                                    >  1:182600/1‑149 (MQ=255)
                                                                   ggCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGaca                                                                                  >  1:21902/1‑148 (MQ=255)
                                                                    gCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGa                                                                                    >  5:307232/1‑145 (MQ=255)
                                                                    gCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGa                                                                                    <  6:307232/145‑1 (MQ=255)
                                                                     cATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACaa                                                                               <  3:193838/149‑1 (MQ=255)
                                                                       ttCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATa                                                                             <  2:182600/149‑1 (MQ=255)
                                                                          gcctgcctGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACaa                                                                               >  3:118161/1‑144 (MQ=255)
                                                                          gcctgcctGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACaa                                                                               <  4:118161/144‑1 (MQ=255)
                                                                          gcctgcctGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTc                                                                          >  4:138405/1‑149 (MQ=255)
                                                                          gcctgcctGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTc                                                                          >  1:200078/1‑149 (MQ=255)
                                                                           cctgcctgGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGACCAGGACACAATAGTCt                                                                         >  5:90641/1‑149 (MQ=255)
                                                                               cctggcctgTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTc                                                                       <  2:24489/147‑1 (MQ=255)
                                                                               cctggcctgTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCg                                                                      <  4:81916/148‑1 (MQ=255)
                                                                                 tggcctgTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCa                                                                   >  7:211578/1‑149 (MQ=255)
                                                                                   gcctgTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATg                                                                 <  6:47075/149‑1 (MQ=255)
                                                                                     ctgTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGcc                                                               >  7:20737/1‑149 (MQ=255)
                                                                                        tCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGt                                                             <  6:196043/148‑1 (MQ=255)
                                                                                          gCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCaa                                                                                                         <  1:134335/102‑1 (MQ=255)
                                                                                          gCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCaa                                                                                                         >  2:134335/1‑102 (MQ=255)
                                                                                           cAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGAcc                                                         <  6:168145/149‑1 (MQ=255)
                                                                                           cAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGAcc                                                         >  2:282414/1‑149 (MQ=255)
                                                                                            aGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGACCAGGACACAATAGTCTTCGCCATGCCGTGACCg                                                        <  6:90641/149‑1 (MQ=255)
                                                                                               ccTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATa                                                                             <  7:277368/125‑1 (MQ=255)
                                                                                               ccTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATa                                                                             >  8:277368/1‑125 (MQ=255)
                                                                                                cTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAAc                                                    >  7:31977/1‑149 (MQ=255)
                                                                                                 ttGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACt                                                   >  2:280638/1‑149 (MQ=255)
                                                                                                 ttGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACt                                                   >  5:214937/1‑149 (MQ=255)
                                                                                                   gCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACt                                                   <  3:138405/147‑1 (MQ=255)
                                                                                                     aCGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCgg                                                                                                 <  6:198310/99‑1 (MQ=255)
                                                                                                     aCGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCt                                               >  6:144675/1‑149 (MQ=255)
                                                                                                      cGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTTGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTc                                              <  1:280638/149‑1 (MQ=255)
                                                                                                           cTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGATCGGAACTTGCTCAgcgc                                         <  2:21902/149‑1 (MQ=255)
                                                                                                           cTTGGCAGAGCGGCATCAGCGACCGCCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAgcgc                                         >  3:84853/1‑149 (MQ=255)
                                                                                                              ggCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCt                                       <  2:189475/148‑1 (MQ=255)
                                                                                                               gCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTc                                      <  5:144675/148‑1 (MQ=255)
                                                                                                                 agagCTGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGgt                                   <  2:186797/149‑1 (MQ=255)
                                                                                                                         aTCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGcc                           <  8:211578/149‑1 (MQ=255)
                                                                                                                              cGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCt                                                                         <  8:133529/98‑1 (MQ=255)
                                                                                                                              cGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCt                                                                         >  7:133529/1‑98 (MQ=255)
                                                                                                                                aCCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGTACGGGCCAGGCCACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCaa                     >  4:13229/1‑148 (MQ=255)
                                                                                                                                    tgagcGCCAGGAACTCATGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATc                 >  1:315684/5‑148 (MQ=255)
                                                                                                                                    tCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATc                 >  4:192831/1‑148 (MQ=255)
                                                                                                                                    tCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATc                 >  1:317983/1‑148 (MQ=255)
                                                                                                                                      agcGCCAGGAACTCATGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCaa              <  2:315684/147‑1 (MQ=255)
                                                                                                                                       ccGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAAc             >  5:302896/1‑149 (MQ=255)
                                                                                                                                       ccGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAAc             >  4:195663/1‑149 (MQ=255)
                                                                                                                                          ccATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAACTGc          >  8:45430/1‑149 (MQ=255)
                                                                                                                                              gAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAACTGCCAg       >  1:295790/1‑148 (MQ=255)
                                                                                                                                                  tCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAACTGCAAGAACGc  >  2:53880/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCAGCTTGTGCAAGAACATCAGCTTGGCGCTGATACCGCCGATGACAAGAGCGTTCCATTGAGGTTGGGCATTCGCCTGCCTGGCCTGTCGCAGGCCTTGCACGGTGCTTGGCAGAGCGGCATCAGCGACCGTCACCGCCATGAGTTCCTGGCGCTCAGGTTGGCTCAGCGCAGATTGGCCCACAGCAGCAATGAAGCGGGCACGGGCCAGGACACAATAGTCTTCGCCATGCCGTGACCGGAACTTGCTCAGCGCCTCGGTGTAGCGCCCTGCAAAATCCAACTGCCAGAACGC  >  NC_002947/5778896‑5779190

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