Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,809,156 G→C 8.1% L237L (CTG→CTC) ‡ dsbD‑II → thiol/disulfide interchange protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,809,1560GC 8.1% 188.4 / 6.2 62L237L (CTG→CTC) ‡dsbD‑IIthiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (2/3);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

TGCTGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTC  >  NC_002947/4809008‑4809300
                                                                                                                                                     |                                                                                                                                                
tGCTGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGctgctg                                                                                                                                                  >  4:312214/1‑149 (MQ=255)
   tGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTggg                                                                                                                                                >  3:99097/1‑148 (MQ=255)
     ggCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCa                                                                                                                                              <  1:9624/148‑1 (MQ=255)
          tGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCCCAGCAGCCAgg                                                                                                                                        >  8:78956/1‑149 (MQ=255)
          tGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCCCAGCAGCCAgg                                                                                                                                        <  7:78956/149‑1 (MQ=255)
              gcCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTAGCTGCTGGGCAGCCTGGCggg                                                                                                                                    >  2:292216/1‑149 (MQ=255)
               ccAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCgggg                                                                                                                                   <  8:293963/149‑1 (MQ=255)
                cAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCAGGCTGCCCAGCAGCCAgggttgct                                                                                                                                  <  6:133586/149‑7 (MQ=255)
                 aGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCt                                                                                                                                 <  2:225376/149‑1 (MQ=255)
                     cccGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTc                                                                                                                             <  7:57930/149‑1 (MQ=255)
                        ggCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtg                                                                                                                          <  3:282587/149‑1 (MQ=255)
                            cgGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtgctgc                                                                                                                      <  3:177341/149‑1 (MQ=255)
                             ggctggctGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTgctgct                                                                                                                     >  4:224228/1‑149 (MQ=255)
                               ctggctgCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTgg                                                                                                                   <  4:99097/149‑1 (MQ=255)
                                         gCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCg                                                                                                           <  1:81494/147‑1 (MQ=255)
                                             gggTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCCTGGCTGCAGCAACCCAGGCTGCCCAGCAGCCAgggt                                                                                                                                      <  1:219299/117‑3 (MQ=255)
                                             gggTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCAGGCTGCCCAGCAGCCAgggt                                                                                                                                      >  2:219299/1‑114 (MQ=255)
                                             gggTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTCCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTc                                                                                                     <  7:173536/149‑1 (MQ=255)
                                               gTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCgg                                                                                                   >  7:298304/1‑149 (MQ=255)
                                                    tgtgCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCt                                                                                              >  2:269303/1‑149 (MQ=255)
                                                        cTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCg                                                                                           >  3:9356/1‑148 (MQ=255)
                                                         tGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGaa                                                                                         <  4:9356/149‑1 (MQ=255)
                                                           aGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGaact                                                                                       <  4:285825/149‑1 (MQ=255)
                                                             cATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGaactg                                                                                      >  8:128635/1‑148 (MQ=255)
                                                                       gtgtACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctgc                                                                           <  3:224228/149‑1 (MQ=255)
                                                                       gtgtACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctgc                                                                           >  8:18411/1‑149 (MQ=255)
                                                                          tACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtgct                                                                                                                        <  6:44044/101‑1 (MQ=255)
                                                                          tACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtgct                                                                                                                        >  5:44044/1‑101 (MQ=255)
                                                                             gCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctg                                                                            <  4:215510/142‑1 (MQ=255)
                                                                             gCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctg                                                                            >  3:215510/1‑142 (MQ=255)
                                                                                     gggggTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTAGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGAcc                                                             <  1:292216/149‑1 (MQ=255)
                                                                                            gtTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGa                                                      >  3:74041/1‑149 (MQ=255)
                                                                                                ctgcCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGAcc                                                             >  8:105926/1‑138 (MQ=255)
                                                                                                ctgcCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGAcc                                                             <  7:105926/138‑1 (MQ=255)
                                                                                                 tgcCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACcgcg                                                 >  1:181589/1‑149 (MQ=255)
                                                                                                 tgcCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACcgcg                                                 >  1:313434/1‑149 (MQ=255)
                                                                                                    ccTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGcccg                                                                     <  5:221142/126‑2 (MQ=255)
                                                                                                    ccTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGcccg                                                                     >  6:221142/1‑125 (MQ=255)
                                                                                                         tgGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGc                                           >  5:7008/1‑147 (MQ=255)
                                                                                                          gggTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAg                                        <  8:128033/149‑1 (MQ=255)
                                                                                                           ggTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAg                                        <  2:268088/148‑1 (MQ=255)
                                                                                                               ccAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTAcc                                   >  1:139663/1‑149 (MQ=255)
                                                                                                                cAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGc                                                           <  4:180170/124‑1 (MQ=255)
                                                                                                                cAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGc                                                           >  3:180170/1‑124 (MQ=255)
                                                                                                                  gCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGtt                                                                                                                              <  8:295118/55‑1 (MQ=255)
                                                                                                                  gCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGtt                                                                                                                              >  7:295118/1‑55 (MQ=255)
                                                                                                                    cTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTgg                              >  8:68509/1‑149 (MQ=255)
                                                                                                                      gCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTggcg                            >  3:249131/1‑149 (MQ=255)
                                                                                                                       cAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTggcg                            <  2:139663/148‑1 (MQ=255)
                                                                                                                       cAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTggcg                            <  2:181589/148‑1 (MQ=255)
                                                                                                                          gCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCaa                        >  1:193488/1‑149 (MQ=255)
                                                                                                                            ccTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAAc                       >  4:296573/1‑149 (MQ=255)
                                                                                                                             cTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGGTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTGTTGGAACTGCAACTGCCTGCCGCCCTGCGCGCCCGCCTGGACCGCGCCCGGCAAGGTACCCGTGGCGGCCAcc                      >  6:23845/1‑149 (MQ=255)
                                                                                                                                gCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTg                    <  7:68509/148‑1 (MQ=255)
                                                                                                                                 cTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTa                  <  5:7776/149‑1 (MQ=255)
                                                                                                                                 cTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTa                  <  7:128635/149‑1 (MQ=255)
                                                                                                                                          cccTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCgg          >  1:244790/1‑148 (MQ=255)
                                                                                                                                          cccTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCgg          >  2:323090/1‑148 (MQ=255)
                                                                                                                                           ccTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCa        <  7:18411/149‑1 (MQ=255)
                                                                                                                                            cTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCAt       >  1:69292/1‑149 (MQ=255)
                                                                                                                                            cTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGACGATGTACGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGCCCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCAt       >  7:239887/1‑149 (MQ=255)
                                                                                                                                               gctgctGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCAtt      <  6:7008/147‑1 (MQ=255)
                                                                                                                                               gctgctGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTg     <  2:244790/148‑1 (MQ=255)
                                                                                                                                               gctgctGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGc    <  3:296573/149‑1 (MQ=255)
                                                                                                                                                 tgctgGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTc  >  6:6820/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TGCTGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTC  >  NC_002947/4809008‑4809300

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: