Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,123,943 T→A 8.5% intergenic (‑138/+83) PP_mr13 ← / ← gcvT‑I cspA regulatory region/aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,123,9430TA 8.5% 180.0 / 6.7 59intergenic (‑138/+83)PP_mr13/gcvT‑IcspA regulatory region/aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base A (3/2);  total (33/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.88e-01

GCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCT  >  NC_002947/1123799‑1124084
                                                                                                                                                |                                                                                                                                             
gCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcct                                                                                                                                          >  1:1516075/1‑150 (MQ=255)
        cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                  <  2:1447939/150‑1 (MQ=255)
         aCGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATg                                                                                                                                 <  1:1246888/150‑1 (MQ=255)
           gAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGtt                                                                                                                               >  2:1531147/1‑150 (MQ=255)
                     gACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAt                                                                                                                     >  1:1582996/1‑150 (MQ=255)
                      aCCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAtt                                                                                                                    >  2:103925/1‑150 (MQ=255)
                       cccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTg                                                                                                                   >  1:153168/1‑150 (MQ=255)
                              aaaCACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGt                                                                                                            >  2:50668/1‑150 (MQ=255)
                               aaCACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCACAAGTTATTGGGCAGGtt                                                                                                           >  1:157547/1‑150 (MQ=255)
                                acacAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGCTTc                                                                                                          <  2:153168/150‑1 (MQ=255)
                                     gtgcattgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGaaaa                                                                                                     <  2:1582996/150‑1 (MQ=255)
                                      tgcagtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAAc                                                                                                    <  1:1081917/150‑1 (MQ=255)
                                               tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGACATGGCGCCAGGCCAGGCGCCATgcc                                                                                                                               >  1:239886/1‑112 (MQ=255)
                                               tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGACATGGCGCCAGGCCAGGCGCCATgcc                                                                                                                               <  2:239886/114‑3 (MQ=255)
                                                 tACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcc                                                                                                                                           >  2:901896/1‑100 (MQ=255)
                                                 tACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcc                                                                                                                                           <  1:901896/100‑1 (MQ=255)
                                                      cAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtgctg                                                                                    >  1:314164/1‑150 (MQ=255)
                                                          tctcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtg                                                                                       <  2:1281870/143‑1 (MQ=255)
                                                          tctcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtg                                                                                       >  1:1281870/1‑143 (MQ=255)
                                                          tctcGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCa                                                                                >  1:954057/1‑150 (MQ=255)
                                                                 cGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGAGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCtt                                                                         <  2:729758/150‑1 (MQ=255)
                                                                  gAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcct                                                                                                                                          <  2:346105/84‑1 (MQ=255)
                                                                  gAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctggcct                                                                                                                                          >  1:346105/1‑84 (MQ=255)
                                                                   aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                <  1:615280/125‑1 (MQ=255)
                                                                   aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                >  2:615280/1‑125 (MQ=255)
                                                                   aaaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTca                                                                       <  2:1516075/150‑1 (MQ=255)
                                                                    aaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCctggcctggc                                                                                                                                       <  1:293565/85‑1 (MQ=255)
                                                                    aaaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCctggcctggc                                                                                                                                       >  2:293565/1‑85 (MQ=255)
                                                                     aaaCAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGTATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGTCCAAAACCTTcaca                                                                     >  2:1185236/1‑150 (MQ=255)
                                                                          aGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc                                                                                                                >  1:833666/1‑102 (MQ=255)
                                                                          aGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc                                                                                                                <  2:833666/102‑1 (MQ=255)
                                                                              cTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGttt                                                            <  2:314164/150‑1 (MQ=255)
                                                                                         ggAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgta                                                 >  2:106827/1‑150 (MQ=255)
                                                                                          gAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgtag                                                <  1:1531147/150‑1 (MQ=255)
                                                                                          gAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGgtag                                                >  2:1298671/1‑150 (MQ=255)
                                                                                             aaCTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTaa                                             >  1:440247/1‑150 (MQ=255)
                                                                                                    ggTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtggtg                                      >  2:1485325/1‑150 (MQ=255)
                                                                                                      tttGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTc                                    <  1:1298671/150‑1 (MQ=255)
                                                                                                          ccGATAACCGACCGGTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtgctgc                                                                                   <  2:949789/99‑1 (MQ=255)
                                                                                                          ccGATAACCGACCGGTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtgctgc                                                                                   >  1:949789/1‑99 (MQ=255)
                                                                                                            gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa                              >  1:1053696/1‑150 (MQ=255)
                                                                                                                  cgaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                >  2:1303301/1‑78 (MQ=255)
                                                                                                                  cgaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAat                                                                                                <  1:1303301/78‑1 (MQ=255)
                                                                                                                  cgaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCtt                                                                         >  1:1017413/1‑101 (MQ=255)
                                                                                                                  cgaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCtt                                                                         <  2:1017413/101‑1 (MQ=255)
                                                                                                                   gaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtggt                                       <  1:120331/134‑1 (MQ=255)
                                                                                                                   gaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtggt                                       >  2:120331/1‑134 (MQ=255)
                                                                                                                   gaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg                                         <  2:1363403/132‑1 (MQ=255)
                                                                                                                   gaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg                                         >  1:1363403/1‑132 (MQ=255)
                                                                                                                   gaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCAttt                       >  2:560696/1‑150 (MQ=255)
                                                                                                                    accgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTc                                    >  2:1253721/1‑136 (MQ=255)
                                                                                                                    accgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTc                                    <  1:1253721/136‑1 (MQ=255)
                                                                                                                       gaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCt                   <  2:1494972/150‑1 (MQ=255)
                                                                                                                        aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTaa                                             <  1:1464751/123‑1 (MQ=255)
                                                                                                                        aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTaa                                             >  2:1464751/1‑123 (MQ=255)
                                                                                                                             tCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGttt                                                            >  2:600539/1‑103 (MQ=255)
                                                                                                                             tCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGttt                                                            <  1:600539/103‑1 (MQ=255)
                                                                                                                                gAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCaa          >  2:829974/1‑150 (MQ=255)
                                                                                                                                  aaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAtt                                                                                                                    >  2:1599443/1‑42 (MQ=255)
                                                                                                                                  aaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAtt                                                                                                                    <  1:1599443/42‑1 (MQ=255)
                                                                                                                                     aacATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCa     >  1:1414829/3‑150 (MQ=255)
                                                                                                                                      gCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTCCAACGGGATGTCCTGAAACTAATATGTCTATGCTGCCAGTAGGCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAc    <  1:1185236/150‑1 (MQ=255)
                                                                                                                                       cATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAcc   <  1:50668/150‑1 (MQ=255)
                                                                                                                                       cATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCAcc   >  2:718758/1‑150 (MQ=255)
                                                                                                                                        aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTACCCAGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct  >  1:15843/1‑150 (MQ=255)
                                                                                                                                        aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct  >  2:1621981/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCT  >  NC_002947/1123799‑1124084

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