Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 971,679 A→G 6.4% intergenic (+34/‑30) yajC → / → secD protein translocase subunit YajC/protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947971,6790AG 6.4% 179.0 / 2.5 63intergenic (+34/‑30)yajC/secDprotein translocase subunit YajC/protein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base G (2/2);  total (33/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 = 8.63e-01

TCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGATCAGCGGT  >  NC_002947/971533‑971825
                                                                                                                                                   |                                                                                                                                                  
tCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAAcg                                                                                                                                                  >  2:229502/1‑149 (MQ=255)
tCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAAcg                                                                                                                                                  >  3:130880/1‑149 (MQ=255)
      tttCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGcccc                                                                                                                                            <  6:229227/149‑1 (MQ=255)
          cGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGt                                                                                                                                        >  1:147287/1‑149 (MQ=255)
               tcTTCGGGGTGGTGGATGTCTCAGACAGTGTCGAGCTGTAGCTCCAGCTGGGCGACATTGCTGCGACCCTGACAAATGTGACGCTCAAGGATATCTGAGTCAGGGTTTTTATTCTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTg                                                                                                                                    <  8:110644/147‑1 (MQ=255)
                cTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGaa                                                                                                                                  >  4:135875/1‑149 (MQ=255)
                 ttCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGTTGCGCGCCCcgact                                                                                                                                        <  4:86136/142‑4 (MQ=255)
                 ttCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGTTGCGCGCCCcgact                                                                                                                                        >  3:86136/1‑139 (MQ=255)
                 ttCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAAc                                                                                                                                 >  5:131596/1‑149 (MQ=255)
                 ttCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAAc                                                                                                                                 >  3:9862/1‑149 (MQ=255)
                    gtggtgCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACggg                                                                                                                              <  6:131596/149‑1 (MQ=255)
                             gAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTCTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGt                                                                                                                        <  1:229502/147‑1 (MQ=255)
                              aaGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATg                                                                                                                    >  4:200640/1‑149 (MQ=255)
                                gtgtCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCt                                                                                                                  <  5:121474/149‑1 (MQ=255)
                                   tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTTCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                               <  8:295222/149‑1 (MQ=255)
                                   tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                               <  7:14661/149‑1 (MQ=255)
                                   tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGa                                                                                                                >  3:141956/1‑148 (MQ=255)
                                   tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGa                                                                                                                <  6:80712/148‑1 (MQ=255)
                                           tGTCTAGCTTAATTTCCAGAAGGGTGCCATTGCTGCGACCCGGCCAAAAGGTACGCTCAAGGCTATCTTAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATg                                                                                                                    <  6:131169/136‑1 (MQ=255)
                                           tGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATc                                                                                                       >  3:161977/1‑149 (MQ=255)
                                           tGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATc                                                                                                       >  3:204696/1‑149 (MQ=255)
                                           tGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCGGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATg                                                                                                                    >  5:131169/1‑136 (MQ=255)
                                             tCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACgg                                                                                                                               <  7:260817/123‑1 (MQ=255)
                                             tCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACgg                                                                                                                               >  8:260817/1‑123 (MQ=255)
                                             tCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATc                                                                                                       <  3:200640/147‑1 (MQ=255)
                                                  cTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGGCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtgtg                                                                                                <  6:109615/149‑1 (MQ=255)
                                                       gTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAat                                                                                           <  4:130880/149‑1 (MQ=255)
                                                          ccAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATg                                                                                        >  8:182889/1‑149 (MQ=255)
                                                          ccAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATg                                                                                        >  8:247081/1‑149 (MQ=255)
                                                            aGAAGGTTGCCATTGCTGCGACCCTGCAAAAAGGTACGCTCATGGCTATCTTAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCa                                                                                      <  4:141956/149‑1 (MQ=255)
                                                                               gACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACtgg                                                                   <  4:161977/149‑1 (MQ=255)
                                                                                    tGCCAAAAGGTGCGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTc                                                              >  1:206511/1‑149 (MQ=255)
                                                                                    tGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATg                                                                                        >  5:166616/1‑123 (MQ=255)
                                                                                    tGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATg                                                                                        <  6:166616/123‑1 (MQ=255)
                                                                                     gCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTgg                                                                         <  7:58293/137‑1 (MQ=255)
                                                                                     gCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTgg                                                                         >  8:58293/1‑137 (MQ=255)
                                                                                        aaaaGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTGTGTGGAAATATGCACTGCTCGGGGGGGTACGGGTGGTCGGtt                                                          >  4:17507/1‑149 (MQ=255)
                                                                                        aaaaGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtgt                                                                                                 >  4:157976/1‑110 (MQ=255)
                                                                                        aaaaGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtgt                                                                                                 <  3:157976/110‑1 (MQ=255)
                                                                                        aaaaGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAatat                                                                                         >  3:296973/1‑118 (MQ=255)
                                                                                           cggTACGCTCAAGGCTATCTGAGTTACGGTTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCa                                                       <  5:47572/148‑1 (MQ=255)
                                                                                           aGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCa                                                       <  3:135875/149‑1 (MQ=255)
                                                                                                 gCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTAtt                                                 >  1:303394/1‑149 (MQ=255)
                                                                                                  cTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGCGGCGCACCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTc                                                <  2:303394/149‑1 (MQ=255)
                                                                                                   tCAAGGCTATCTGAGTTACCGGTATTATTTTCCAGTCGGGGAGCTCAAAGCGCCCCGCGTCTAGAACGGGCTGCGTGATGATGAACAAATATACTCTGTGGAAATATGAACTGATCGTGCTGTTATTGGTGGTAGGTTTCATTTAttac                                               >  3:210268/1‑147 (MQ=255)
                                                                                                        gCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTcc                                          <  4:204696/149‑1 (MQ=255)
                                                                                                         cTATCTGAGTTACCGGTTTTATTTTCCCGGCGGGGCGCGAACCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGGCGGTTTCATTTATTCCGCTccc                                         <  4:9862/149‑1 (MQ=255)
                                                                                                            tctttgTTTCCGTTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCAGTGATCGTGCTGGTACTGGTTGTCGGTTTCATTTATTCCGATCCCaa                                       <  4:210268/143‑1 (MQ=255)
                                                                                                            tCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTAt                                                  <  8:88100/137‑1 (MQ=255)
                                                                                                            tCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTAt                                                  >  7:88100/1‑137 (MQ=255)
                                                                                                            tCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTCTGGAAATATGCACTGATCGTGCAGGTACTGGTGGTCGGTTTCATTTATTCCGCTACCaa                                       >  3:188831/1‑148 (MQ=255)
                                                                                                               gAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCAtt                                                     >  2:118251/1‑131 (MQ=255)
                                                                                                               gAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCAtt                                                     <  1:118251/131‑1 (MQ=255)
                                                                                                                  ttACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTAc                                >  3:212990/1‑149 (MQ=255)
                                                                                                                               tttCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGGTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCcgg                   >  2:40572/1‑149 (MQ=255)
                                                                                                                                      gTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa              >  6:264733/1‑148 (MQ=255)
                                                                                                                                      gTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTAc             <  7:182889/149‑1 (MQ=255)
                                                                                                                                      cgcGGGGCGCGTTGCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACctc                                   >  7:216693/3‑127 (MQ=255)
                                                                                                                                      cgcGGGGCGCGTTGCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACctc                                   <  8:216693/125‑1 (MQ=255)
                                                                                                                                           ggCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTATTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGAt         >  2:198346/1‑148 (MQ=255)
                                                                                                                                               cgcAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCg                                                                               >  4:129611/1‑74 (MQ=255)
                                                                                                                                               cgcAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCg                                                                               <  3:129611/74‑1 (MQ=255)
                                                                                                                                                 cAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGATCAGCGGt  >  6:297493/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGATCAGCGGT  >  NC_002947/971533‑971825

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: