Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,110,085 C→G 6.4% N621K (AAC→AAG PP_3616 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,110,0850CG 6.4% 179.5 / 4.6 62N621K (AAC→AAGPP_3616membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/33);  new base G (2/2);  total (27/35)
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.65e-01

TTACGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGG  >  NC_002947/4109943‑4110234
                                                                                                                                              |                                                                                                                                                     
ttACGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTg                                                                                                                                                <  1:1464283/150‑1 (MQ=32)
   cGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCAGCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGt                                                                                                                                             <  2:889690/150‑1 (MQ=17)
         aCCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCa                                                                                                                                       >  2:716503/1‑150 (MQ=31)
                   ccAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTAACCTCcggcgg                                                                                                                             >  1:230347/1‑150 (MQ=17)
                    cAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggt                                                                                                                            <  1:1311538/150‑1 (MQ=32)
                     aggaggTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtg                                                                                                                           <  1:716503/150‑1 (MQ=32)
                         ggTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtgg                                                                                                                       <  1:309653/150‑1 (MQ=32)
                          gTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggt                                                                                                                      <  1:653167/150‑1 (MQ=32)
                           tGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGAGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtg                                                                                                                     <  2:490092/150‑1 (MQ=17)
                             gatgttgaTTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtggt                                                                                                                   <  2:588995/143‑1 (MQ=17)
                             cAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtggt                                                                                                                   <  1:490241/150‑1 (MQ=32)
                                cGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCAATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTa                                                                                                                >  1:148504/1‑150 (MQ=17)
                                 ggTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTag                                                                                                               <  2:1509598/150‑2 (MQ=17)
                                  gTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtg                                                                                                                     >  1:907485/1‑143 (MQ=32)
                                  gTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtg                                                                                                                     <  2:907485/143‑1 (MQ=32)
                                       tGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         <  2:230347/150‑1 (MQ=32)
                                       tGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         <  1:1483740/150‑1 (MQ=32)
                                       tGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         <  1:516868/150‑1 (MQ=32)
                                        gAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCggg                                                                                                        >  1:771376/1‑150 (MQ=32)
                                               ggCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCtg                                                                                                 <  1:26214/150‑1 (MQ=32)
                                                           gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggt                                                                                                                      >  2:6172/1‑117 (MQ=33)
                                                           gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggt                                                                                                                      <  1:6172/117‑1 (MQ=33)
                                                           gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATc                                                                                     <  2:941859/150‑1 (MQ=32)
                                                                agagATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATc                                                                                <  1:720723/150‑1 (MQ=32)
                                                                   gaTGAAAGCGTTGCGCATCTCGGGTAGCCGCATAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCt                                                                             >  1:279165/1‑150 (MQ=17)
                                                                    aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTaaa                                                                                                              >  1:120067/1‑116 (MQ=33)
                                                                    aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTaaa                                                                                                              <  2:120067/116‑1 (MQ=33)
                                                                    aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTc                                                                            >  1:18844/1‑150 (MQ=32)
                                                                             ttGCGCATCTCGGGTAGCCGCATAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGcc                                                                   <  2:279165/150‑1 (MQ=17)
                                                                               gcgcATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTa                                                                 <  2:18844/150‑1 (MQ=32)
                                                                                   aTCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaac                                                             >  1:1040016/1‑150 (MQ=32)
                                                                                    tctcGGATAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaaca                                                            >  2:602404/1‑150 (MQ=17)
                                                                                         ggTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         <  1:284784/100‑1 (MQ=34)
                                                                                         ggTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         >  2:284784/1‑100 (MQ=34)
                                                                                                gCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                                                >  1:1268128/1‑150 (MQ=35)
                                                                                                gCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                                                <  2:1040016/150‑1 (MQ=35)
                                                                                                  acaGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGc                                              >  1:1458902/1‑150 (MQ=35)
                                                                                                         ccATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGTTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         <  1:1330393/84‑1 (MQ=25)
                                                                                                         ccATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                         >  2:1330393/1‑84 (MQ=35)
                                                                                                         ccATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCAc                                       >  1:1612695/1‑150 (MQ=35)
                                                                                                          cATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCAcc                                      >  1:672044/1‑150 (MQ=35)
                                                                                                           aTCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCg                                     <  2:1458902/150‑1 (MQ=35)
                                                                                                                  tGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGc                                                          >  2:156551/1‑122 (MQ=32)
                                                                                                                  tGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGc                                                          <  1:156551/122‑1 (MQ=33)
                                                                                                                    accaccACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACGGCCTGCAGAAGAGCGATGCACCGGAATACAcc                            <  2:1268128/150‑1 (MQ=25)
                                                                                                                     ccaccaCTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACAccc                           >  2:59465/1‑150 (MQ=35)
                                                                                                                     ccaccaCTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCGCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGATCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACAccc                           <  1:1446760/150‑1 (MQ=21)
                                                                                                                      caccgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGc                                                      <  1:1455195‑M1/105‑1 (MQ=255)
                                                                                                                      caccgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGc                                                      >  2:1455195‑M1/18‑122 (MQ=255)
                                                                                                                       accgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGa                                                                        >  2:1609045‑M1/17‑103 (MQ=255)
                                                                                                                       accgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGa                                                                        <  1:1609045‑M1/87‑1 (MQ=255)
                                                                                                                       accacTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTg                         >  2:1577703/1‑150 (MQ=35)
                                                                                                                               gcgTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGc                 >  1:1386341/1‑150 (MQ=35)
                                                                                                                                 gTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCtt               >  1:1058588/1‑150 (MQ=35)
                                                                                                                                     gAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTac                                                                >  1:1465648/1‑97 (MQ=33)
                                                                                                                                     gAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTac                                                                <  2:1465648/97‑1 (MQ=34)
                                                                                                                                        ACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCtg                                                                                                 <  2:272975‑M1/61‑1 (MQ=255)
                                                                                                                                        ACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCtg                                                                                                 >  1:272975‑M1/1‑61 (MQ=255)
                                                                                                                                          CCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATc                                                                                <  1:891902‑M1/76‑1 (MQ=255)
                                                                                                                                          CCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATc                                                                                >  2:891902‑M1/1‑76 (MQ=255)
                                                                                                                                           CAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGtt     <  1:1450306‑M1/150‑1 (MQ=255)
                                                                                                                                             ACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCg   <  2:487406‑M1/150‑1 (MQ=255)
                                                                                                                                              cTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCTTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCgg  <  1:34593/150‑1 (MQ=38)
                                                                                                                                              |                                                                                                                                                     
TTACGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGG  >  NC_002947/4109943‑4110234

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