Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,809,157 G→A 7.9% G238S (GGC→AGC)  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,1570GA 7.9% 183.0 / 4.4 63G238S (GGC→AGC) dsbD‑IIthiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/3);  total (30/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 = 9.34e-01

TGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTCGGGGC  >  NC_002947/4809011‑4809305
                                                                                                                                                   |                                                                                                                                                    
tGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTggg                                                                                                                                                     >  3:99097/1‑148 (MQ=255)
  ggCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCa                                                                                                                                                   <  1:9624/148‑1 (MQ=255)
       tGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCCCAGCAGCCAgg                                                                                                                                             >  8:78956/1‑149 (MQ=255)
       tGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCCCAGCAGCCAgg                                                                                                                                             <  7:78956/149‑1 (MQ=255)
           gcCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTAGCTGCTGGGCAGCCTGGCggg                                                                                                                                         >  2:292216/1‑149 (MQ=255)
            ccAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCgggg                                                                                                                                        <  8:293963/149‑1 (MQ=255)
             cAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCAGGCTGCCCAGCAGCCAgggttgct                                                                                                                                       <  6:133586/149‑7 (MQ=255)
              aGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCt                                                                                                                                      <  2:225376/149‑1 (MQ=255)
                  cccGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTc                                                                                                                                  <  7:57930/149‑1 (MQ=255)
                     ggCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtg                                                                                                                               <  3:282587/149‑1 (MQ=255)
                         cgGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtgctgc                                                                                                                           <  3:177341/149‑1 (MQ=255)
                          ggctggctGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTgctgct                                                                                                                          >  4:224228/1‑149 (MQ=255)
                            ctggctgCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTgg                                                                                                                        <  4:99097/149‑1 (MQ=255)
                                      gCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCg                                                                                                                <  1:81494/147‑1 (MQ=255)
                                          gggTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGGCCCTGGCTGCAGCAACCCAGGCTGCCCAGCAGCCAgggt                                                                                                                                           <  1:219299/117‑3 (MQ=255)
                                          gggTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCAGGCTGCCCAGCAGCCAgggt                                                                                                                                           >  2:219299/1‑114 (MQ=255)
                                          gggTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTCCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTc                                                                                                          <  7:173536/149‑1 (MQ=255)
                                            gTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCgg                                                                                                        >  7:298304/1‑149 (MQ=255)
                                                 tgtgCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCt                                                                                                   >  2:269303/1‑149 (MQ=255)
                                                     cTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCg                                                                                                >  3:9356/1‑148 (MQ=255)
                                                      tGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGaa                                                                                              <  4:9356/149‑1 (MQ=255)
                                                        aGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGaact                                                                                            <  4:285825/149‑1 (MQ=255)
                                                          cATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGaactg                                                                                           >  8:128635/1‑148 (MQ=255)
                                                                    gtgtACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctgc                                                                                >  8:18411/1‑149 (MQ=255)
                                                                    gtgtACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctgc                                                                                <  3:224228/149‑1 (MQ=255)
                                                                       tACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtgct                                                                                                                             >  5:44044/1‑101 (MQ=255)
                                                                       tACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGtgct                                                                                                                             <  6:44044/101‑1 (MQ=255)
                                                                          gCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctg                                                                                 <  4:215510/142‑1 (MQ=255)
                                                                          gCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAActgcctg                                                                                 >  3:215510/1‑142 (MQ=255)
                                                                                  gggggTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTAGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGAcc                                                                  <  1:292216/149‑1 (MQ=255)
                                                                                         gtTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGa                                                           >  3:74041/1‑149 (MQ=255)
                                                                                             ctgcCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGAcc                                                                  >  8:105926/1‑138 (MQ=255)
                                                                                             ctgcCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGAcc                                                                  <  7:105926/138‑1 (MQ=255)
                                                                                              tgcCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACcgcg                                                      >  1:181589/1‑149 (MQ=255)
                                                                                              tgcCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACcgcg                                                      >  1:313434/1‑149 (MQ=255)
                                                                                                 ccTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGcccg                                                                          <  5:221142/126‑2 (MQ=255)
                                                                                                 ccTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGcccg                                                                          >  6:221142/1‑125 (MQ=255)
                                                                                                      tgGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGc                                                >  5:7008/1‑147 (MQ=255)
                                                                                                       gggTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAg                                             <  8:128033/149‑1 (MQ=255)
                                                                                                        ggTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAg                                             <  2:268088/148‑1 (MQ=255)
                                                                                                            ccAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTAcc                                        >  1:139663/1‑149 (MQ=255)
                                                                                                             cAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGc                                                                >  3:180170/1‑124 (MQ=255)
                                                                                                             cAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGc                                                                <  4:180170/124‑1 (MQ=255)
                                                                                                               gCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGtt                                                                                                                                   >  7:295118/1‑55 (MQ=255)
                                                                                                               gCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGtt                                                                                                                                   <  8:295118/55‑1 (MQ=255)
                                                                                                                 cTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTgg                                   >  8:68509/1‑149 (MQ=255)
                                                                                                                   gCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTggcg                                 >  3:249131/1‑149 (MQ=255)
                                                                                                                    cAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTggcg                                 <  2:181589/148‑1 (MQ=255)
                                                                                                                    cAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTggcg                                 <  2:139663/148‑1 (MQ=255)
                                                                                                                       g‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCaa                             >  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                       <  7:128635/149‑1 (MQ=255)
                                                                                                                              cTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTa                       <  5:7776/149‑1 (MQ=255)
                                                                                                                                       cccTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCgg               >  2:323090/1‑148 (MQ=255)
                                                                                                                                       cccTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCgg               >  1:244790/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)
                                                                                                                                                   ggCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTCGGGGc  <  2:69292/149‑1 (MQ=255)
                                                                                                                                                   ggCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTCGGGGc  <  2:193488/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGGCAATGGCGCCAGTGCCCGGCGCGGCTGGCTGCTGGCCGGGGTCTATGTGCTGAGCATGGCCCTGGTGTACGCCGGGCTGGGGGTGGTTGCTGCCCTGCTGGGTGCCAGCCTGCAGG‑CCTGGCTGCAGCAACCCTGGCTGCTGGGCAGCCTGGCGGGGCTGTTCGTGCTGCTGGCCCTGCCGATGTTCGGCGTCTTCGAACTGCAACTGCCTGCCGCCCTGCGCGACCGCCTGGACCGCGCCGGGCAAGGTACCCGTGGCGGCAACCTGTATGGTGCGGCATTGCTCGGGGC  >  NC_002947/4809011‑4809305

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