Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,353,367 C→T 7.4% E250E (GAG→GAA)  pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,353,3670CT 7.4% 198.4 / 5.0 67E250E (GAG→GAA) pnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (3/2);  total (35/32)
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.62e-01

GCCGAAGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACC  >  NC_002947/5353220‑5353507
                                                                                                                                                    |                                                                                                                                            
gCCGAAGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACt                                                                                                                                             >  6:132900/1‑149 (MQ=255)
       aTCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACCGGGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGgcgc                                                                                                                                       >  1:52046/1‑149 (MQ=255)
          tctTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGa                                                                                                                                    >  7:10271/1‑148 (MQ=255)
          tctTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGa                                                                                                                                    >  2:198587/1‑148 (MQ=255)
            ttttACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaagc                                                                                                                                    <  8:163627/146‑5 (MQ=255)
            ttttACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaagc                                                                                                                                    >  7:163627/1‑142 (MQ=255)
                 cTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCAtt                                                                                                                             >  2:35970/1‑148 (MQ=255)
                 cTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTg                                                                                                                            <  1:198587/149‑1 (MQ=255)
                  ttCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaagcgg                                                                                                                                  <  6:270810/142‑7 (MQ=255)
                  ttCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTg                                                                                                                            >  7:143468/1‑148 (MQ=255)
                  ttCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGGGATGGTGTAACCCTGCGAAACCGCTTCGCCGAGTTCGgcgaagcgg                                                                                                                                  >  5:270810/1‑136 (MQ=255)
                    cGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGaa                                                                                                                          <  8:10271/148‑1 (MQ=255)
                    cGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGaa                                                                                                                          >  7:222378/1‑148 (MQ=255)
                     ggAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACa                                                                                                                        <  4:234883/149‑1 (MQ=255)
                           cGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCgg                                                                                                                  >  2:4540/1‑149 (MQ=255)
                             aaGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGgtg                                                                                                                <  6:12866/149‑1 (MQ=255)
                               ggCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGGAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGtt                                                                                                              <  2:35834/149‑1 (MQ=255)
                                  tGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCTAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTgg                                                                                                            >  1:267189/1‑148 (MQ=255)
                                     ccTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAAc                                                                                                        <  2:52046/149‑1 (MQ=255)
                                          ttcACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAg                                                                                                   >  4:182691/1‑149 (MQ=255)
                                               c‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTc                                                                                              <  1:35970/149‑1 (MQ=255)
                                                    gaACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGt                                                                                          >  2:44175/1‑149 (MQ=255)
                                                        ttGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACgg                                                                                                      >  7:180683/1‑133 (MQ=255)
                                                        ttGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACgg                                                                                                      <  8:180683/133‑1 (MQ=255)
                                                            cAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTTGCAACGGCAGGTTTCCAGTCCCAAGTc                                                                                  <  1:4540/149‑1 (MQ=255)
                                                             aaCCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCg                                                                                 >  1:139401/1‑149 (MQ=255)
                                                             aaCCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCg                                                                                 >  5:143108/1‑149 (MQ=255)
                                                             aaCCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCg                                                                                 >  4:152479/1‑149 (MQ=255)
                                                              aCCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCgg                                                                                >  2:176124/1‑149 (MQ=255)
                                                                   cTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGc                                                                           <  8:222378/149‑1 (MQ=255)
                                                                   cTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGc                                                                           <  3:182691/149‑1 (MQ=255)
                                                                       tCACGCAGCTCACCGAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACTGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGc                                                                       <  4:245974/149‑1 (MQ=255)
                                                                       tCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCGAGc                                                                       <  5:132900/149‑1 (MQ=255)
                                                                           gCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCg                                                                   <  3:233276/149‑1 (MQ=255)
                                                                            cAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCgg                                                                  <  8:143468/149‑1 (MQ=255)
                                                                                  aCCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagc                                                              >  8:105549/1‑147 (MQ=255)
                                                                                  aCCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagc                                                              >  7:211738/1‑147 (MQ=255)
                                                                                    ccAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCtt                                                                     >  2:43489/1‑138 (MQ=255)
                                                                                    ccAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCtt                                                                     <  1:43489/138‑1 (MQ=255)
                                                                                    ccAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCAGTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGcagccagc                                                          >  5:220268/1‑149 (MQ=255)
                                                                                      aGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGctc                                                        <  7:74688/149‑1 (MQ=255)
                                                                                          gcgcgTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCttt                                                     >  2:56878/1‑148 (MQ=255)
                                                                                          gcgcgTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTg                                                    <  4:232594/149‑1 (MQ=255)
                                                                                           cgcgTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGa                                                   <  6:143108/149‑1 (MQ=255)
                                                                                            gcgTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGa                                                   >  8:113491/1‑148 (MQ=255)
                                                                                                 gCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCt                                             <  1:176124/149‑1 (MQ=255)
                                                                                                    gTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACATACTggc                                          >  4:14634/1‑148 (MQ=255)
                                                                                                      cGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTc                                                                                  >  6:219048/1‑107 (MQ=255)
                                                                                                      cGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTc                                                                                  <  5:219048/107‑1 (MQ=255)
                                                                                                               cGGTGATGGTGTAACCCTGCGAAACCGCTTCGACGTACTCGGCGCGGATAGCATTGATCAGCGCGGTGTTGGCAACGGCAGGTTTCCAGTCCGAAGGCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTgg                               <  2:267189/149‑1 (MQ=255)
                                                                                                                 gTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGTACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGaa                             <  3:152479/149‑1 (MQ=255)
                                                                                                                  tGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAAt                            >  5:3726/1‑149 (MQ=255)
                                                                                                                         gtAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTg                                >  8:184905/1‑138 (MQ=255)
                                                                                                                         gtAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTg                                <  7:184905/138‑1 (MQ=255)
                                                                                                                          tAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTg                    >  7:139895/1‑149 (MQ=255)
                                                                                                                           aaCCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTg                    <  8:139895/148‑1 (MQ=255)
                                                                                                                              ccTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGa                                          <  4:73856/123‑1 (MQ=255)
                                                                                                                              ccTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGa                                          >  3:73856/1‑123 (MQ=255)
                                                                                                                                  cGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGaa              <  7:105549/147‑1 (MQ=255)
                                                                                                                                   gAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTCGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTATTTGGCATCCTGGATAACAGCCTGGAATTAGTCGTGGGCGAACAg           >  8:53374/1‑149 (MQ=255)
                                                                                                                                       ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTg                                >  2:173431/1‑124 (MQ=255)
                                                                                                                                       ccGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTg                                <  1:173431/124‑1 (MQ=255)
                                                                                                                                         tccgcGCCGAGTTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGATCAGTACgg      >  2:48105/5‑148 (MQ=255)
                                                                                                                                           ttCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCAc   >  6:30574/1‑149 (MQ=255)
                                                                                                                                             cGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCAcc  <  8:211738/148‑1 (MQ=255)
                                                                                                                                                cGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACg       >  4:71312/1‑140 (MQ=255)
                                                                                                                                                cGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACg       <  3:71312/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GCCGAAGATCTCTTTTACTTCGGAAGCCGAAGGCTGGCCTTCTTCACC‑GGAGAACTTGGCAACCGCCTGGTCACGCAGCTCACCCAGACGCGCGTAGCGGTCGGCTTTGACGGTGATGGTGTAACCCTGCGAAACCGCTTCGCCGAACTCGGCGCGGATAGCATTGAACAGCTCGGTGTTGGCAACGGCAGGTTTCCAGTCCCAAGTCGGCTTGCCAGCTTCGGCAGCCAGCTCTTTGACAGCCTGGATAACAGCCTGGAATTCGTCGTGGGCGAACAGTACGGCACC  >  NC_002947/5353220‑5353507

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