Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,597,212 G→A 5.9% A866V (GCG→GTG)  PP_4071 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,597,2120GA 5.9% 207.1 / 2.3 68A866V (GCG→GTG) PP_4071membrane protein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (2/2);  total (32/36)
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.55e-01

ACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCG  >  NC_002947/4597065‑4597356
                                                                                                                                                   |                                                                                                                                                
aCTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGc                                                                                                                                                 >  3:349324/1‑149 (MQ=255)
       gcagcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTg                                                                                                                                           >  5:200061/1‑148 (MQ=255)
       gcagcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTg                                                                                                                                           >  2:148238/1‑148 (MQ=255)
         agcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAAt                                                                                                                                        <  7:350375/149‑1 (MQ=255)
         agcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAAt                                                                                                                                        <  1:189206/149‑1 (MQ=255)
                ccTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGTCGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAgc                                                                                                                                 <  1:345791/149‑1 (MQ=255)
                ccTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAgc                                                                                                                                 >  1:45104/1‑149 (MQ=255)
                      cTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCtgt                                                                                                                           >  7:212775/1‑149 (MQ=255)
                      cTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCtgt                                                                                                                           <  8:354947/149‑1 (MQ=255)
                         ggTCTTGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTGATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCtgttg                                                                                                                         <  4:127566/148‑1 (MQ=255)
                          gTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTgttgtt                                                                                                                       <  6:1907/149‑1 (MQ=255)
                           tCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGttgttg                                                                                                                      <  6:200061/149‑1 (MQ=255)
                               ggCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCt                                                                                                                    <  5:126960/147‑1 (MQ=255)
                               ggCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCa                                                                                                                  <  3:88723/149‑1 (MQ=255)
                               ggCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTACGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCa                                                                                                                  <  8:360422/149‑1 (MQ=255)
                                  tACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGc                                                                                                               >  2:285955/1‑149 (MQ=255)
                                    cccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCtt                                                                                                             >  4:118311/1‑149 (MQ=255)
                                    cccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACACCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCtt                                                                                                             >  3:52661/1‑149 (MQ=255)
                                    cccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACACCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCtt                                                                                                             >  3:47631/1‑149 (MQ=255)
                                     ccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTg                                                                                                            >  2:49149/1‑149 (MQ=255)
                                         tACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCt                                                                                                          <  3:106912/147‑1 (MQ=255)
                                         tACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGt                                                                                                        <  2:45104/149‑1 (MQ=255)
                                         tACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGt                                                                                                        <  8:212775/149‑1 (MQ=255)
                                           ccGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGtt                                                                                                       <  7:247346/148‑1 (MQ=255)
                                           ccGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTg                                                                                                      <  4:349324/149‑1 (MQ=255)
                                             gCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGAt                                                                                                    <  3:118311/149‑1 (MQ=255)
                                              ccTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATc                                                                                                   >  3:175631/1‑149 (MQ=255)
                                               cTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGGTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATc                                                                                                   >  4:297841/1‑148 (MQ=255)
                                                   gCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGaaa                                                                                               <  1:148238/148‑1 (MQ=255)
                                                     aTAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACgg                                                                                            >  6:204171/1‑149 (MQ=255)
                                                     aTAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACgg                                                                                            >  7:119680/1‑149 (MQ=255)
                                                         tCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCg                                                                                         >  6:364552/1‑148 (MQ=255)
                                                          cGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGaa                                                                                       <  8:115903/149‑1 (MQ=255)
                                                            tGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAgg                                                                                     <  5:304004/149‑1 (MQ=255)
                                                                   cGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCCGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCa                                                                              >  8:83940/1‑149 (MQ=255)
                                                                            cTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCTACATGGCGCTCACGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCcgccacg                                                                     >  4:216174/1‑149 (MQ=255)
                                                                                gTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCaa                                                                 <  1:133050/149‑1 (MQ=255)
                                                                                 tGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCaa                                                                 >  8:195331/1‑148 (MQ=255)
                                                                                   ggTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACa                                                               <  5:364552/148‑1 (MQ=255)
                                                                                        cATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAAc                                                          <  1:49149/148‑1 (MQ=255)
                                                                                          tCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCCTCCGCACCGCCAACATTAACCgg                                                       >  6:27746/1‑149 (MQ=255)
                                                                                          tCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCgg                                                       <  4:374167/149‑1 (MQ=255)
                                                                                           cGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGt                                                      >  1:12813/1‑149 (MQ=255)
                                                                                               gtagCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCCGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGAt                                                   <  7:83940/148‑1 (MQ=255)
                                                                                                    tCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGc                                             >  7:171508/1‑149 (MQ=255)
                                                                                                    tCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGc                                             >  2:99994/1‑149 (MQ=255)
                                                                                                      ggTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGc                                             <  5:204171/147‑1 (MQ=255)
                                                                                                       gTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGa                                           >  3:131801/1‑148 (MQ=255)
                                                                                                         ttttCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACACCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACt                                        <  4:52661/149‑1 (MQ=255)
                                                                                                          tttCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACACCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACt                                        <  4:47631/148‑1 (MQ=255)
                                                                                                           ttCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACtg                                       <  1:99994/148‑1 (MQ=255)
                                                                                                            tCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGa                                                                                                     <  4:253058/85‑1 (MQ=255)
                                                                                                             cGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGt                                                                                                        >  8:216303/1‑81 (MQ=255)
                                                                                                             cGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGt                                                                                                        <  7:216303/81‑1 (MQ=255)
                                                                                                             cGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACtgtt                                     <  3:297841/148‑1 (MQ=255)
                                                                                                             cGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACtgtt                                     <  8:218964/148‑1 (MQ=255)
                                                                                                             cGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGAGAAGGGGCCTCATCCGCCACGCCAACATTAACCGCTGATAACAGCGAACtgttg                                    >  3:168939/1‑149 (MQ=255)
                                                                                                              gCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACtgttgt                                   >  8:156703/1‑149 (MQ=255)
                                                                                                                    aGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTAAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGaa                             >  7:84065/1‑149 (MQ=255)
                                                                                                                         gCAACATGGCGCTCAAGCCTGCAAACGCCCGCTTAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGt                         <  8:171508/148‑1 (MQ=255)
                                                                                                                           aaCATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGtctc                      >  2:168798/1‑149 (MQ=255)
                                                                                                                           aaCATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGtctc                      >  2:94078/1‑149 (MQ=255)
                                                                                                                            aCATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTcc                     <  4:175631/149‑1 (MQ=255)
                                                                                                                                   gcTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATc              >  1:53451/1‑149 (MQ=255)
                                                                                                                                       aaGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGt           <  7:156703/148‑1 (MQ=255)
                                                                                                                                               aaaCGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGt                                                      >  4:132855/1‑97 (MQ=255)
                                                                                                                                               aaaCGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGt                                                      <  3:132855/97‑1 (MQ=255)
                                                                                                                                               aaaCGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGgcg  >  4:87732/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCG  >  NC_002947/4597065‑4597356

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