Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,786,527 C→T 7.6% intergenic (+14/‑117) PP_5068 → / → trxC hypothetical protein/thioredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,786,5270CT 7.6% 181.0 / 3.9 66intergenic (+14/‑117)PP_5068/trxChypothetical protein/thioredoxin
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (3/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGGGACTACAGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTGGTAATCCCCTGCCCCCAC  >  NC_002947/5786382‑5786676
                                                                                                                                                 |                                                                                                                                                     
ggggACTACCGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCtgtg                                                                                                                                                   <  1:1091348/150‑1 (MQ=255)
ggggACTACAGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCtgtg                                                                                                                                                   <  2:787823/150‑1 (MQ=255)
    aCTACAGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTCCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAg                                                                                                                                               <  1:434610/150‑1 (MQ=255)
      tACAGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCg                                                                                                                                             >  2:598258/1‑150 (MQ=255)
      tACAGTGAAGTTCGACGGCTGAACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACTGCCGGAGTGGGGCAAGCAACAGGAAATCAGCTGCTATTCCTGATTAGCGCCGAATCATGTGGGAGCg                                                                                                                                             >  1:444285/1‑150 (MQ=255)
        cAGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCggg                                                                                                                                           >  2:360241/1‑150 (MQ=255)
         aGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGc                                                                                                                                          >  1:923527/1‑150 (MQ=255)
                                  tCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAAc                                                                                                                      <  2:716517/145‑1 (MQ=255)
                                  tCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAAc                                                                                                                      >  1:716517/1‑145 (MQ=255)
                                   ccTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATg                                                                                                                <  1:790651/150‑1 (MQ=255)
                                      ggCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGAt                                                                                                             <  1:1088480/150‑1 (MQ=255)
                                            cccaggcgctcggcgggaaggcggtggaggccgttgcagtgggggcaggggattaccagtgggtcggtcatttcatctccggtgttgggttcacggtttcCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCAcc                                                                                                       <  1:517677‑M1/50‑1 (MQ=255)
                                               tttCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGc                                                                                                    >  1:911740/1‑150 (MQ=255)
                                                      gAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCCTCGCCg                                                                                             <  2:923527/150‑1 (MQ=255)
                                                            ctcGGGATGCCAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTc                                                                                       <  2:552169/148‑1 (MQ=255)
                                                               gggATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTcgcg                                                                                    <  1:708274/150‑1 (MQ=255)
                                                                gcggctgtggccgttgcagtgggggcaggggattcccagtgggtcggttagttcatctccggtgttgggttcacggtttcCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCGACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCgg                                                                                   >  2:517677‑M1/81‑150 (MQ=255)
                                                                 gATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCCACGTGCTGGATTGCACCGGCTTCGCCGGTGTTCGTggg                                                                                  <  2:620569/150‑1 (MQ=255)
                                                                 cggttgaggccgttgcagtgggggcaggggattaccagtgggtcggtcatttcatctccggtgttgggttcacggtttcCTTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCggg                                                                                  >  2:1089444‑M1/80‑150 (MQ=255)
                                                                    ttgaggccgttgcagtgggggcaggggattaccagtgggtcggtcatttcatctccggtgttgggttcacggtttcCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCggg                                                                                  <  1:554792‑M1/71‑1 (MQ=255)
                                                                    ttgaggccgttgcagtgggggcaggggattaccagtgggtcggtcatttcatctccggtgttgggttcacggtttcCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCggg                                                                                  >  2:554792‑M1/77‑147 (MQ=255)
                                                                    gCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCAc                                                                               >  1:531953/1‑150 (MQ=255)
                                                                      aaGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGc                                                                             >  2:676328/1‑150 (MQ=255)
                                                                       aGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGcc                                                                            <  1:598258/150‑1 (MQ=255)
                                                                          cTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTGCCGGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGCTTGCACCGGCTTCGCCGGTGGTCGCGGGCACGCCCGc                                                                         >  2:169765/1‑150 (MQ=255)
                                                                           tACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                              <  2:635802/96‑1 (MQ=255)
                                                                           tACGCCGAGCGACCGCCGGAGTGGGACAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                              >  1:635802/1‑96 (MQ=255)
                                                                                 gAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGAtt                                                                                                            <  1:483926/108‑1 (MQ=255)
                                                                                 gAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGAtt                                                                                                            >  2:483926/1‑108 (MQ=255)
                                                                                 gAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCcaca                                                                  >  1:1092383/1‑150 (MQ=255)
                                                                                  gattcggcgctaatcaggaagagcagctgatttccaggtgctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGt                                                                                                                                        <  2:555656‑M2/37‑1 (MQ=255)
                                                                                  gattcggcgctaatcaggaagagcagctgatttccaggtgctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGt                                                                                                                                        >  1:555656‑M2/43‑79 (MQ=255)
                                                                                         aaggggattaccagtgggtcggtcatttcatctccggtgttgggttcacggtttcCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAAcc                                                          >  1:76191‑M1/56‑150 (MQ=255)
                                                                                          ccgGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCg                                                         <  2:531953/150‑1 (MQ=255)
                                                                                             gAGTGGCGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTCGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGa                                                      <  1:714026/150‑1 (MQ=255)
                                                                                              aGTTGGGCAAGCACCTGGAAATCAGCTGCTCGTCCTGATTAGCCCCGAATCCTGTGGGAGCGCGCTTGCCCGCGAAGCAGGCAACGTGATGCATTGCACCGGCTTCAGCGATGTTCGCGGGCCCTCCCGCTCCCACAAGTAAAcctgtaa                                                     <  2:444285/150‑6 (MQ=25)
                                                                                              aGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGACAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGaa                                                     <  2:489099/150‑1 (MQ=255)
                                                                                                   ggCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCaa                                                >  1:724233/1‑150 (MQ=255)
                                                                                                   ggCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCaa                                                <  1:360241/150‑1 (MQ=255)
                                                                                                     cAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAaca                                              <  2:877959/150‑1 (MQ=255)
                                                                                                      gagcagctgatttccaggtgctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCaa                                                                                                                       <  1:868748‑M2/54‑1 (MQ=255)
                                                                                                      gagcagctgatttccaggtgctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCAAAGCAGGCaa                                                                                                                       >  2:868748‑M2/23‑76 (MQ=255)
                                                                                                        GCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCa                                                                                <  1:277268‑M1/113‑1 (MQ=255)
                                                                                                        GCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCa                                                                                >  2:277268‑M1/1‑113 (MQ=255)
                                                                                                        GCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCg                                           >  2:724233‑M1/1‑150 (MQ=255)
                                                                                                             tGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCt                                                                        >  1:879537/1‑116 (MQ=255)
                                                                                                             tGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCt                                                                        <  2:879537/116‑1 (MQ=255)
                                                                                                             tGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAAcc                                                          >  2:1008504/1‑130 (MQ=255)
                                                                                                             tGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAAcc                                                          <  1:1008504/130‑1 (MQ=255)
                                                                                                             ggtcatttcatctccgttgttgggttcacggtttcCTTGTGGGAGCGGGCGTTCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCggg                                                                                  <  2:356382‑M1/71‑1 (MQ=255)
                                                                                                             ggtcatttcatctccggtgttgggttcacggtttcCTTGTGGGAGCGGGCGTTCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCggg                                                                                  >  1:356382‑M1/36‑106 (MQ=255)
                                                                                                             ggtcatttcatctccggtgttgggttcacggtttcCTTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCgg                                                                                   <  2:807791‑M1/70‑1 (MQ=255)
                                                                                                             ggtcatttcatctccggtgttgggttcacggtttcCTTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCgg                                                                                   >  1:807791‑M1/36‑105 (MQ=255)
                                                                                                              ggAAATCAGCTGATCTTCCTGATTAGCGCCGAATCCTGTGGGCGCGGGTGGGCGCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACGGGAGATg                                     >  2:471671/1‑150 (MQ=38)
                                                                                                                      gctgctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCaa                                                >  1:464946/1‑131 (MQ=255)
                                                                                                                      gctgctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCaa                                                <  2:464946/131‑1 (MQ=255)
                                                                                                                          ctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATg                                                                                                                <  1:1016564‑M2/61‑1 (MQ=255)
                                                                                                                          ctCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATg                                                                                                                >  2:1016564‑M2/3‑63 (MQ=255)
                                                                                                                             ttCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCAc                      >  1:58708/1‑150 (MQ=255)
                                                                                                                             ttCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCAc                      <  2:911740/150‑1 (MQ=255)
                                                                                                                               ccTGATTAGCGCCGAATCCTGTGGGAGCGGGTGTGCTCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTg                    <  1:856484/150‑1 (MQ=255)
                                                                                                                               ccTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTg                    <  2:1092383/150‑1 (MQ=255)
                                                                                                                                 tGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCCCGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTGGt                  >  1:1142594/1‑150 (MQ=255)
                                                                                                                                          ccGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTGGTAATCCCCTg         >  1:93324/1‑150 (MQ=255)
                                                                                                                                               tCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAAt                                 >  1:136749/1‑121 (MQ=255)
                                                                                                                                               tCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAAt                                 <  2:136749/121‑1 (MQ=255)
                                                                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTGGTAATCCCCTGCCCCCAc  >  1:787201/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGGGACTACAGTGAAGTTCGACGGCTGCACCAAGTCCTGGCCAGGCCTTTCGAGGAACAGCCCGGGATGCAAGCCTACGCCGAGCGACCGCCGGAGTGGGGCAAGCACCTGGAAATCAGCTGCTCTTCCTGATTAGCGCCGAATCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGGCAACGTGATGGATTGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAGGAAACCGTGAACCCAACACCGGAGATGAAATGACCGACCCACTGGTAATCCCCTGCCCCCAC  >  NC_002947/5786382‑5786676

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: