Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,332,293 C→G 6.3% intergenic (+28/‑21) dksA → / → gluQ RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,332,2930CG 6.3% 190.7 / 2.7 62intergenic (+28/‑21)dksA/gluQRNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base G (2/2);  total (36/26)
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.56e-01

CGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCC  >  NC_002947/5332148‑5332433
                                                                                                                                                  |                                                                                                                                             
cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATgc                                                                                                                                           <  1:554215/150‑1 (MQ=255)
      gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCAACTGAAAGGAGAATGCCGAGAAAAAGG‑AAAATTCGGGCGACAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                                     >  2:273738/1‑150 (MQ=255)
           gCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCa                                                                                                                                <  1:303172/150‑1 (MQ=255)
             aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                              >  1:861972/1‑150 (MQ=255)
              tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                                             >  1:792683/1‑150 (MQ=255)
              tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                                             <  1:142897/150‑1 (MQ=255)
                 agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                                          <  1:244332/150‑1 (MQ=255)
                  gaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGTAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                                      >  2:1418811/1‑137 (MQ=255)
                  gaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                                      <  1:1418811/137‑1 (MQ=255)
                       gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGTGCGTGGCGCATCATGCGCCCCGGTTCATTTATATGCCCa                                                                                                                    >  2:1330033/1‑150 (MQ=255)
                       gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTta                                                                                                                            <  1:176272/142‑1 (MQ=255)
                       gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTta                                                                                                                            >  2:176272/1‑142 (MQ=255)
                       gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCa                                                                                                                    <  1:1532026/150‑1 (MQ=255)
                               cgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAcc                                                                                                            <  1:643661/150‑1 (MQ=255)
                                 ccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                          <  2:792683/150‑1 (MQ=255)
                                     gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                                      <  2:964404/150‑1 (MQ=255)
                                        gcccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                                           >  1:271988/1‑110 (MQ=255)
                                        ccccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                                           <  2:271988/109‑1 (MQ=255)
                                         cccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGACGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACt                                                                                                  >  1:1610255/1‑150 (MQ=255)
                                          ccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTc                                                                                                 >  2:916696/1‑150 (MQ=255)
                                           cgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCCTTTATATGCCCAGCCTGACCATGACCGACGcc                                                                                                >  1:509798/1‑150 (MQ=255)
                                             cccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTCATGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAg                                                                                              <  1:1330033/150‑1 (MQ=255)
                                              ccgaCCTCCGACCTGTCCTTCGACTCCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGc                                                                                             <  2:509798/147‑1 (MQ=255)
                                               cTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                               >  2:1293061/1‑115 (MQ=255)
                                               cTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                                               <  1:1293061/115‑1 (MQ=255)
                                                          cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGgcgc                                                                                >  1:1245952/1‑150 (MQ=255)
                                                          cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGAAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCAGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGgcgc                                                                                >  1:1043623/1‑150 (MQ=255)
                                                            gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                                            >  2:1369080/1‑137 (MQ=255)
                                                            gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                                            <  1:1369080/137‑1 (MQ=255)
                                                               cTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGc                                                                           >  1:833747/1‑150 (MQ=255)
                                                                   gACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGtt                                                                                                                                   <  1:157722/91‑1 (MQ=255)
                                                                   gACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGtt                                                                                                                                   >  2:157722/1‑91 (MQ=255)
                                                                        cAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGcccccacccc                                                                  >  2:834750/1‑150 (MQ=255)
                                                                          aGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGcccccaccccca                                                                <  2:1245952/150‑1 (MQ=255)
                                                                           ggAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAg                                                               >  2:1691855/1‑150 (MQ=255)
                                                                            gacaTTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGc                                                              >  1:2400849/4‑150 (MQ=255)
                                                                              aaTTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCgg                                                            >  2:418643/1‑150 (MQ=255)
                                                                                    cGAGAAAAAGGCAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGCGGCGCATCATGCGCCCCGTTTCCTTTATATGCCCAGCCTGACCCTGACCGACTCCCGCTACA‑TCGCGCGCTTCGCCCCCACCCCAAGCAGCTTcc                                                       >  2:280348/1‑150 (MQ=255)
                                                                                    cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                      >  2:1496531/1‑150 (MQ=255)
                                                                                    cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                                      <  1:1199635/150‑1 (MQ=255)
                                                                                    cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGc                                                                                                                  <  2:2151935/91‑1 (MQ=255)
                                                                                    cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGc                                                                                                                  >  1:2151935/1‑91 (MQ=255)
                                                                                                aaaaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCt                                           >  2:297028/1‑150 (MQ=255)
                                                                                                aaaaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACCGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCt                                           <  2:1650317/150‑1 (MQ=255)
                                                                                                   aaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGc                                        <  1:1792509/150‑1 (MQ=255)
                                                                                                    aaCGGTCGGCAAAGGCGGATCTTTCTTCCACGTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCt                                       >  1:195894/1‑150 (MQ=255)
                                                                                                        gTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTc                                   <  1:916696/150‑1 (MQ=255)
                                                                                                             cAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgccgc                              >  1:264455/1‑150 (MQ=255)
                                                                                                             cAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgccgc                              >  1:2046528/1‑150 (MQ=255)
                                                                                                                  gCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCATGCACTTCGGCTCGCTGGTCGCCGCCCTCg                         >  2:1671964/1‑150 (MQ=255)
                                                                                                                    tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCTCAGCTTAAAAAATAAGCCCAGCCGGGCCATGACCTACTACAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc                       >  1:1855434/1‑150 (MQ=255)
                                                                                                                      aTCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc                     >  1:274991/1‑150 (MQ=255)
                                                                                                                             ttcCCCCTGCACGGGACGCTTGCTGCGCACATTGTTATTGATATGCCCACCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc              <  2:264455/150‑1 (MQ=255)
                                                                                                                             ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATTCACAAGC‑GCTCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc              >  1:1955778/1‑150 (MQ=255)
                                                                                                                             ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc              >  2:995191/1‑150 (MQ=255)
                                                                                                                             ttcCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTc              <  2:833747/150‑1 (MQ=255)
                                                                                                                               ccACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGa            >  2:565316/1‑150 (MQ=255)
                                                                                                                                   cTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGcc        >  1:1718206/1‑150 (MQ=255)
                                                                                                                                       aCGGGGTGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcgcg    >  1:314965/1‑150 (MQ=255)
                                                                                                                                       aCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCcgcg    >  1:182298/1‑150 (MQ=255)
                                                                                                                                         ggggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGCGATGATCGACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGcc  <  1:1233600/150‑1 (MQ=255)
                                                                                                                                         ggggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGCCTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCAATGCCCGCGcc  <  2:1610255/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACA‑TCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCC  >  NC_002947/5332148‑5332433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑