Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,004,815 C→T 6.8% R466C (CGC→TGC) ‡ PP_2627 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,004,8150CT 6.8% 212.8 / 4.2 73R466C (CGC→TGC) ‡PP_2627membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base T (3/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

AGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGTC  >  NC_002947/3004672‑3004958
                                                                                                                                               |                                                                                                                                                 
aGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCATAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCa                                                                                                                                               <  8:152575/148‑1 (MQ=255)
       cacGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCtgt                                                                                                                                       >  8:175614/1‑149 (MQ=255)
         cgacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgg                                                                                                                                     <  7:175614/149‑1 (MQ=255)
         cgacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgg                                                                                                                                     <  2:72145/149‑1 (MQ=255)
         cgacgaAGCTGCCACTGCGTGAAGGTAAGGCACAGGTGAAGTCGTTTGTCTATTTCATGCTGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACTATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgg                                                                                                                                     <  2:276878/149‑1 (MQ=255)
          gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTg                                                                                                                                           >  1:195476/1‑142 (MQ=255)
          gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTg                                                                                                                                           <  2:195476/142‑1 (MQ=255)
          gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtggg                                                                                                                                    >  1:379661/1‑149 (MQ=255)
            cgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtg                                                                                                                                  <  2:379661/149‑1 (MQ=255)
             gaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtgg                                                                                                                                 <  6:71926/149‑1 (MQ=255)
             gaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtgg                                                                                                                                 <  3:330075/149‑1 (MQ=255)
              aaGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGtgggtggg                                                                                                                                >  4:184081/1‑149 (MQ=255)
                   gCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTg                                                                                                                           >  3:80979/1‑149 (MQ=255)
                   gCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTg                                                                                                                           >  1:217849/1‑149 (MQ=255)
                      ccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTTTTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGgtg                                                                                                                        >  6:227325/1‑149 (MQ=255)
                      ccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGgtg                                                                                                                        <  4:80979/149‑1 (MQ=255)
                      ccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGAATTCAGGCGGGTTATCTATCAGGAAGCCGGGATGGCGGCTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGATCTCgtg                                                                                                                        >  5:245356/1‑149 (MQ=255)
                                ggTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCt                                                                                                              >  6:91999/1‑149 (MQ=255)
                                 gTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTg                                                                                                             <  3:229213/149‑1 (MQ=255)
                                 gTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTg                                                                                                             >  3:86518/1‑149 (MQ=255)
                                   aaGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑c                                                                                                         >  2:336762/1‑149 (MQ=255)
                                    agtcacaGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGACACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑Cg                                                                                                        >  6:292060/4‑149 (MQ=255)
                                        acaGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCa                                                                                                    <  2:217849/149‑1 (MQ=255)
                                        acaGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCa                                                                                                    <  8:207999/149‑1 (MQ=255)
                                            gggAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAg                                                                                                                     >  7:131815/1‑130 (MQ=255)
                                            gggAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAg                                                                                                                     <  8:131815/130‑1 (MQ=255)
                                               aaGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGc                                                                                             >  3:295213/1‑149 (MQ=255)
                                                     ttGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtg                                                                                                                                  <  6:159871/108‑1 (MQ=255)
                                                             tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTgg                                                                                >  6:302237/1‑148 (MQ=255)
                                                             tatCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTgg                                                                                >  2:63510/3‑148 (MQ=255)
                                                                 aTCCAGCAAGCATTCAGGCGTGTTTTTTATCTGGAAGGTGGGCTGGGGGGTGACAATGTTCGGGTTGGGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGgtg                                                                                                                        <  4:144216/106‑1 (MQ=255)
                                                                 aTCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGgtg                                                                                                                        >  3:144216/1‑106 (MQ=255)
                                                                 aTCCAGCAAGCATTCAGGCGGGCTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCg                                                                           >  5:349163/1‑149 (MQ=255)
                                                                  tCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGt                                                                          >  6:323843/1‑149 (MQ=255)
                                                                    cAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGtt                                                                         <  1:336762/148‑1 (MQ=255)
                                                                          gcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGa                                                                  >  4:367519/1‑149 (MQ=255)
                                                                            aTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGc                                                                     <  3:245641/144‑1 (MQ=255)
                                                                            aTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGc                                                                     >  4:245641/1‑144 (MQ=255)
                                                                            aTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGAc                                                                 <  7:81814/148‑1 (MQ=255)
                                                                            aTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGAc                                                                 >  8:81814/1‑148 (MQ=255)
                                                                                 ggCGGGGTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCaa                                                            <  3:367519/148‑1 (MQ=255)
                                                                                  gCGGGCTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATg                                                          <  6:349163/149‑1 (MQ=255)
                                                                                      gTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGATGCCAGCTGGCAGCGTTACGCCGACATCAATGgcgc                                                      <  3:184081/149‑1 (MQ=255)
                                                                                       ttATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCg                                                                           >  6:122350/1‑127 (MQ=255)
                                                                                            tatcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGc                                                >  3:295429/1‑149 (MQ=255)
                                                                                             atcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgc                                                        >  7:207444/1‑140 (MQ=255)
                                                                                             atcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgc                                                        <  8:207444/140‑1 (MQ=255)
                                                                                             atcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCt                                               >  8:108395/1‑149 (MQ=255)
                                                                                              tcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTg                                              >  8:104133/1‑149 (MQ=255)
                                                                                              tcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTg                                              <  1:63510/149‑1 (MQ=255)
                                                                                              tcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTAAGCCGACATCAATGGCGCCAAGGCTg                                              <  2:31895/149‑1 (MQ=255)
                                                                                                  gAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACAt                                                               >  8:219873/1‑128 (MQ=255)
                                                                                                  gAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                           >  4:359839/1‑132 (MQ=255)
                                                                                                  gAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                           <  3:359839/132‑1 (MQ=255)
                                                                                                   aaGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAgt                                         <  5:227325/149‑1 (MQ=255)
                                                                                                              gCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGcc                                <  4:295429/147‑1 (MQ=255)
                                                                                                                     aCAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGTAGGTTCAGTATCGCCGTTGCCAGCTGGAAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGc                          >  1:112622/1‑148 (MQ=255)
                                                                                                                     aCAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                       <  5:302237/149‑1 (MQ=255)
                                                                                                                     aCAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCGGGTGTATGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                       <  6:245356/149‑1 (MQ=255)
                                                                                                                       aaTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                       >  1:142943/1‑147 (MQ=255)
                                                                                                                       aaTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGa                     <  2:142943/149‑1 (MQ=255)
                                                                                                                              gggTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGTCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGc              <  5:323843/149‑1 (MQ=255)
                                                                                                                              gggTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                           >  7:81784/1‑104 (MQ=255)
                                                                                                                              gggTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                           <  8:81784/104‑1 (MQ=255)
                                                                                                                                     gcgcAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCAc        >  7:15471/1‑148 (MQ=255)
                                                                                                                                     gcgcAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCAc        >  2:170067/1‑148 (MQ=255)
                                                                                                                                     ccacagcaAGTGGCGCTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCACTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCCGTCAc        >  5:19292/4‑148 (MQ=255)
                                                                                                                                     ccacagcaAGTGGCGCTTGCTGCGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTc                                   >  4:367707/4‑121 (MQ=255)
                                                                                                                                     ccacagcaAGTGGCGCTTGCTGCGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTc                                   <  3:367707/118‑1 (MQ=255)
                                                                                                                                       gcagcaAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACg       >  1:331681/1‑147 (MQ=255)
                                                                                                                                       acagcaAGTGGCGCTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgc                                                        <  6:290791/97‑1 (MQ=255)
                                                                                                                                       acagcaAGTGGCGCTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgc                                                        >  5:290791/2‑98 (MQ=255)
                                                                                                                                          gcaAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGTc  >  7:297136/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
AGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGG‑‑CGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGTC  >  NC_002947/3004672‑3004958

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