Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,829,359 C→G 10.0% G77A (GGC→GCC)  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,3590CG10.0% 178.7 / 9.0 60G77A (GGC→GCC) PP_3382gluconate 2‑dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base G (3/3);  total (30/30)
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.41e-01

ATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCG  >  NC_002947/3829212‑3829501
                                                                                                                                                    |                                                                                                                                               
aTGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGc                                                                                                                                                 <  2:67948/148‑1 (MQ=255)
 tGAAATAGGCGTAGAGCGCCTGCATCTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCa                                                                                                                                               <  4:83123/149‑1 (MQ=255)
    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                                                                                                  >  4:38490/1‑79 (MQ=255)
                                                                                                                   ggTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCgg                                                                                                  <  3:38490/79‑1 (MQ=255)
                                                                                                                    gTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTCTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCgg                           >  8:157655/1‑149 (MQ=255)
                                                                                                                       gAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGgcgc                        >  6:72077/1‑149 (MQ=255)
                                                                                                                        aaGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCg                                                                                       <  5:111401/85‑1 (MQ=255)
                                                                                                                        aaGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCg                                                                                       >  6:111401/1‑85 (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                 <  5:72077/148‑1 (MQ=255)
                                                                                                                                cGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                 <  1:144929/148‑1 (MQ=255)
                                                                                                                                cGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAg                 <  2:224380/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)
                                                                                                                                                    |                                                                                                                                               
ATGAAATAGGCGTAGAGCGCCTGCATGTCCTGCTCGCTGACCCGCGCATACGACGGGTATGGCATGGCCGGATAGAGCGTACTGCCGTCCTTGCCGATCCCCCTGCGTACGGCCTGGTCGAAGTCTTCGAAGCTGTACTGGCCAATGCCACTGGCCGCCGGGGTGATGTTGGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCG  >  NC_002947/3829212‑3829501

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: