Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,417,744 T→G 6.1% E99D (GAA→GAC) ‡ PP_3030 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,417,7440TG 6.1% 207.1 / 4.7 66E99D (GAA→GAC) ‡PP_3030hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base G (2/2);  total (33/33)
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.44e-01

GGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGG  >  NC_002947/3417602‑3417887
                                                                                                                                                |                                                                                                                                                 
ggTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCgcag                                                                                                                                             <  3:185583/149‑1 (MQ=255)
ggTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCgca                                                                                                                                              <  4:169416/148‑1 (MQ=255)
  tGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCgcagca                                                                                                                                           <  8:260996/149‑1 (MQ=255)
   gCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCGAAGcagcag                                                                                                                                          <  5:198215/149‑1 (MQ=255)
        taCCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCGAAGcagcaggca                                                                                                                                       <  1:79650/148‑4 (MQ=255)
         aCCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtt                                                                                                                                    >  7:56972/1‑149 (MQ=255)
          ccccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtt                                                                                                                                    <  3:91677/148‑1 (MQ=255)
           cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtttt                                                                                                                                  >  4:211914/1‑149 (MQ=255)
           cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtttt                                                                                                                                  >  5:12312/1‑149 (MQ=255)
           cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtttt                                                                                                                                  >  3:293723/1‑149 (MQ=255)
           cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACtttt                                                                                                                                  >  2:22245/1‑149 (MQ=255)
            ccACGGCAGGCTGCGCGTGCGGGCGACTTATAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGGAGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACttttt                                                                                                                                 >  7:207/1‑149 (MQ=255)
                     ggtGCGCGTGCGGGCGAGTTTTAGGGATTAGGGTTGCATGAGGGGGCGGTGCCTCGCGCAGGTTGTTGTGGGGTATGTGTGCCTGGGCGGGGCTGCTGATGGGGGGAATTATGCCTGCTGGCTTCGCAGCAGAGACTTTTTGTTTTTTg                                                                                                                         <  1:210046/147‑1 (MQ=255)
                       tGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTgcgc                                                                                                                      >  5:162965/1‑149 (MQ=255)
                           cgTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCg                                                                                                                  <  1:73311/149‑1 (MQ=255)
                             tGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGac                                                                                                               <  4:293723/149‑1 (MQ=255)
                                          tctcGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCTGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTc                                                                                                   <  1:40236/145‑1 (MQ=255)
                                              gATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCTGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGa                                                                                             <  3:260853/149‑1 (MQ=255)
                                              gATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTc                                                                                                   >  2:40236/1‑144 (MQ=255)
                                                   gCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCg                                                                                          >  8:215324/1‑147 (MQ=255)
                                                      ttGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTGTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAATCg                                                                                     >  8:146646/1‑149 (MQ=255)
                                                        gCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCTCCGCGAAgtcc                                                                                    <  6:12312/148‑4 (MQ=255)
                                                          aTCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAg                                                                                 >  2:225374/1‑149 (MQ=255)
                                                          aTCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAg                                                                                 >  2:266358/1‑149 (MQ=255)
                                                            cAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTg                                                                               >  5:276367/1‑149 (MQ=255)
                                                             agagggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTgg                                                                              >  3:349824/4‑149 (MQ=255)
                                                               cgggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGt                                                                                >  3:52791/1‑145 (MQ=255)
                                                               cgggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGt                                                                                <  4:52791/145‑1 (MQ=255)
                                                                   cggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCAc                                                                        <  3:211914/149‑1 (MQ=255)
                                                                     ggtccTCGCGCAGCTTGTTGTGGGGTCTCAGTGGCTCGGGGGGGCTGCTGATGGGGGGAATATTGCCTGCTGCTTCGCAGCAGAGAGTTTTTGTTTTTTGCGCTTCGACTACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACtt                                                                      <  8:207/146‑1 (MQ=255)
                                                                     gTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACtt                                                                      <  6:276367/149‑1 (MQ=255)
                                                                      tGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCAGCGGAGCGAAGCGGGAGTGGCCTCACtt                                                                      <  1:22245/149‑1 (MQ=255)
                                                                          tCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTc                                                                 >  3:341027/1‑149 (MQ=255)
                                                                          tCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTc                                                                 >  5:131764/1‑149 (MQ=255)
                                                                            gcgcAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTc                                                                 <  8:56972/147‑1 (MQ=255)
                                                                            gcgcAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCt                                                               >  6:166640/1‑149 (MQ=255)
                                                                               cAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCtgtg                                                            >  1:100705/1‑149 (MQ=255)
                                                                                  cTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTg                                                         <  4:349824/149‑1 (MQ=255)
                                                                                     gttgtGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTc                                                       <  5:131758/148‑1 (MQ=255)
                                                                                        gtgGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTc                                                                                                                   >  6:98353/1‑87 (MQ=255)
                                                                                        gtgGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTc                                                                                                                   <  5:98353/87‑1 (MQ=255)
                                                                                         tgGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTg                                                  <  5:166640/149‑1 (MQ=255)
                                                                                           gggTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGcc                                                >  8:115924/1‑149 (MQ=255)
                                                                                             gTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCg                                              >  4:150762/1‑149 (MQ=255)
                                                                                              tATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCg                                              >  1:310508/1‑148 (MQ=255)
                                                                                              tATCAGTCCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGc                                             <  7:115924/149‑1 (MQ=255)
                                                                                                   gTGCCTCAGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTccc                                        >  1:1967/1‑149 (MQ=255)
                                                                                                        tCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGCCTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCgg                                   >  6:99683/1‑149 (MQ=255)
                                                                                                          ggcggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCgg                                   <  5:222651/147‑1 (MQ=255)
                                                                                                          ggcggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCgg                                   <  7:215324/147‑1 (MQ=255)
                                                                                                          agcggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCggg                                  <  2:1967/147‑1 (MQ=255)
                                                                                                           gcggGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGt                                <  3:150762/149‑1 (MQ=255)
                                                                                                              ggggCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCTGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGt                                >  6:71909/1‑147 (MQ=255)
                                                                                                               gggCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAg                             <  2:116474/149‑1 (MQ=255)
                                                                                                                              ggAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGc                                             <  3:190408/118‑1 (MQ=255)
                                                                                                                              ggAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGc                                             >  4:190408/1‑118 (MQ=255)
                                                                                                                                 ataatGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCtgtg                                                            <  2:222854/100‑1 (MQ=255)
                                                                                                                                    gtctctgctgCGAAGCAGCAGGGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCat        >  6:222070/5‑149 (MQ=255)
                                                                                                                                    atGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCat        >  6:40508/1‑149 (MQ=255)
                                                                                                                                    atGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCa         <  6:131764/148‑1 (MQ=255)
                                                                                                                                       tctgctgCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCatgat     >  4:40383/2‑149 (MQ=255)
                                                                                                                                       ccTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGAACAACACGGACTTCGCGGACCGAAGCGGGAGTGGCATCACTTGAGTCATGTGATTTCGGGTGCCCGCTTCCATCGGGGGTCAGCAATTgtgg                    >  1:68091/1‑135 (MQ=255)
                                                                                                                                       ccTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                  >  2:44969/1‑136 (MQ=255)
                                                                                                                                       ccTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                  <  1:44969/136‑1 (MQ=255)
                                                                                                                                       ccTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtgg                    <  2:68091/134‑1 (MQ=255)
                                                                                                                                        ctgctgCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCatgatg    >  4:17906/1‑149 (MQ=255)
                                                                                                                                           ctgctTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATggg  <  6:162965/148‑1 (MQ=255)
                                                                                                                                                ttCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATggg  >  6:13331/1‑144 (MQ=255)
                                                                                                                                                ttCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATggg  <  5:13331/144‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCGGGGCTGCTGATGGGGGGAATAATGCCTGCTGCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGACAACACGGACTTCGCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGG  >  NC_002947/3417602‑3417887

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