Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,829,362 C→T 10.0% S76N (AGT→AAT)  PP_3382 ← gluconate 2‑dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,829,3620CT10.0% 165.7 / 7.2 60S76N (AGT→AAT) PP_3382gluconate 2‑dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base T (3/3);  total (31/29)
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.42e-01

AATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGC  >  NC_002947/3829216‑3829510
                                                                                                                                                   |                                                                                                                                                     
aaTAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTg                                                                                                                                                     >  5:151520/1‑149 (MQ=255)
          gagCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAGTGGCATTGgccagtaca                                                                                                                                            <  5:159338/148‑7 (MQ=255)
          gaccgcgTGGATGTCCTGCTCGCGGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGGCTGGGCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCgggg                                                                                                                                          <  5:168180/146‑1 (MQ=255)
                  gCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTg                                                                                                                                                     >  6:332169/1‑131 (MQ=255)
                  gCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTg                                                                                                                                                     <  5:332169/131‑1 (MQ=255)
                   cATGTCCTGCTCGCTGACCCGCCCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtg                                                                                                                                  <  2:150544/148‑1 (MQ=255)
                     tGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtg                                                                                                                               >  2:52049/1‑149 (MQ=255)
                       tCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTtggtgg                                                                                                                              <  7:182947/148‑1 (MQ=255)
                       tCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGa                                                                                                                             <  2:53215/149‑1 (MQ=255)
                         cTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAg                                                                                                                            >  6:289829/1‑148 (MQ=255)
                         cTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAg                                                                                                                            >  4:46961/1‑148 (MQ=255)
                         cTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGt                                                                                                                           <  3:283881/149‑1 (MQ=255)
                                       gcgcATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATc                                                                                                             <  8:48010/149‑1 (MQ=255)
                                          cATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGc                                                                                                          <  7:285352/149‑1 (MQ=255)
                                               gacgGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTc                                                                                                     <  8:54096/149‑1 (MQ=255)
                                               gacgGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTc                                                                                                     <  3:195671/149‑1 (MQ=255)
                                                acgGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTGGGTGGAGTACACCGTGCCAGTCGGCGtagg                                                                                                     >  8:118910/1‑145 (MQ=255)
                                                          catggcCGGATAGAGCTTACTGCCGTGCTTGGCGTTCCCCGTGCGTAGGGCCTGGTCTTAGTCTTCGAAGGTGTTCTGGCCAATGCCACTGGCCGGCGGGGTGATGTTGTTGGAg                                                                                                                            <  7:142803/115‑1 (MQ=255)
                                                           atggcCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAAccc                                                                                         >  1:132703/1‑149 (MQ=255)
                                                                cGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATc                                                                                                             >  7:252723/1‑124 (MQ=255)
                                                                cGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATc                                                                                                             <  8:252723/124‑1 (MQ=255)
                                                                 ggATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGc                                                                                   <  3:46961/149‑1 (MQ=255)
                                                                    tAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCaaa                                                                                <  8:45708/149‑1 (MQ=255)
                                                                    tAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCaaa                                                                                <  6:151520/149‑1 (MQ=255)
                                                                       agCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAggg                                                                             >  6:330505/1‑149 (MQ=255)
                                                                             cTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTgc                                                                        <  5:289829/148‑1 (MQ=255)
                                                                                ccGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAGAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTgccgc                                                                     >  1:224380/1‑149 (MQ=255)
                                                                                ccGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTgccgcc                                                                    >  1:303370/1‑149 (MQ=255)
                                                                                               ccccTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCTGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATgg                                                       >  2:280515/1‑147 (MQ=255)
                                                                                                  cTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACAGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCCAAGCGTTGGCCGCCTTTGGGGGTATGGCAGGc                                                  >  8:274246/1‑149 (MQ=255)
                                                                                                      gTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGc                                                      <  7:299954/141‑1 (MQ=255)
                                                                                                      gTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGc                                                      >  8:299954/1‑141 (MQ=255)
                                                                                                        aCGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACg                                              >  2:144929/1‑147 (MQ=255)
                                                                                                           gCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGTTATGGCAGGCCACGCAGTc                                         <  5:330505/149‑1 (MQ=255)
                                                                                                            ccTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGt                                          >  1:234482/1‑148 (MQ=255)
                                                                                                               ggTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCgg                                                                                                           <  3:38490/79‑1 (MQ=255)
                                                                                                               ggTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCgg                                                                                                           >  4:38490/1‑79 (MQ=255)
                                                                                                                gTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTCTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCgg                                    >  8:157655/1‑149 (MQ=255)
                                                                                                                   gAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGgcgc                                 >  6:72077/1‑149 (MQ=255)
                                                                                                                    aaGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCg                                                                                                >  6:111401/1‑85 (MQ=255)
                                                                                                                    aaGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCg                                                                                                <  5:111401/85‑1 (MQ=255)
                                                                                                                      aGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTTGTGGAGTACACCTTGCCATTCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCgg                               <  7:118910/149‑1 (MQ=255)
                                                                                                                       gTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCggg                              >  1:222803/1‑149 (MQ=255)
                                                                                                                          ttCGAAGCTGTACTGGCCAATGCCACTGGCCGGCGGGGTGATGTTGGTGGAGTCCAAGCTGCCAATCGGCGTTTCCATCGGAAACCCGCCGGAAAAGGGTTTGCCGGCTTTGGCGGTATGGCAGGCCACGAAGGCCCCGGCGCGGGc                             >  4:221913/1‑147 (MQ=255)
                                                                                                                          ttCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGc                             >  8:133250/1‑147 (MQ=255)
                                                                                                                            cGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                          <  1:144929/148‑1 (MQ=255)
                                                                                                                            cGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                          <  2:224380/148‑1 (MQ=255)
                                                                                                                            cGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                          <  5:72077/148‑1 (MQ=255)
                                                                                                                            cGAAGCTGGACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCa                           <  5:309398/147‑1 (MQ=255)
                                                                                                                                  cccggcggCCAGTGGCATTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                          >  6:81793/7‑142 (MQ=255)
                                                                                                                                  cccggcggCCAGTGGCATTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                          <  5:81793/136‑1 (MQ=255)
                                                                                                                                  cccggcggCCAGTGGCATTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCg                   >  1:101732/7‑149 (MQ=255)
                                                                                                                                   ccggcggCCAGTGGCATTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCGTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGc                  >  5:318286/6‑149 (MQ=255)
                                                                                                                                       tGGTCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAg              >  6:303823/1‑149 (MQ=255)
                                                                                                                                       tGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAg              >  3:164437/1‑149 (MQ=255)
                                                                                                                                         gCCAGTGGCATTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGc            <  2:101732/149‑1 (MQ=255)
                                                                                                                                         gCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGc            >  4:183705/1‑149 (MQ=255)
                                                                                                                                         gCCAATGCCACTGGCCGCCGGGGTGATGGTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAgg             <  2:303370/148‑1 (MQ=255)
                                                                                                                                          ccAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCg           >  3:133397/1‑149 (MQ=255)
                                                                                                                                                 cACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTgc    >  2:187775/1‑149 (MQ=255)
                                                                                                                                                   cTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTgcgc  <  2:234482/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGC  >  NC_002947/3829216‑3829510

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: