Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,637,326 A→T 9.1% intergenic (+35/+31) pdxJ → / ← PP_1437 pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,637,3260AT 9.1% 175.6 / 4.9 63intergenic (+35/+31)pdxJ/PP_1437pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (3/3);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAA  >  NC_002947/1637181‑1637461
                                                                                                                                                 |                                                                                                                                       
tgTATCAATGAGTTGCACATTGGCCATGCGCTGGTGGCGCATGCCGTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGTCAAGcc                                                                                                                                      <  8:279104/148‑1 (MQ=255)
  tATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcg                                                                                                                                    <  7:94088/149‑1 (MQ=255)
     cAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgcg                                                                                                                                  >  6:330579/1‑148 (MQ=255)
        tGAGTTGAACATTGGCGATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTCTGGCCGAGATGTAGTCGCTGATCGTGGCGGCTTCGCTCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCTGGCAAGCCCGCGAaga                                                                                                                              <  5:221858/149‑1 (MQ=255)
          aGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCTTGCCCGCGAaaag                                                                                                                             <  1:213546/148‑1 (MQ=255)
          aGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAagag                                                                                                                             <  6:136933/148‑1 (MQ=255)
              gAACTTTGGCCTTGCGCGGGTGGCGCATGCCTGGTTTGGGGGCTTCACGGCTGCTGGGGCCGATATCCTGGCGGTTAGCGGGGCGGCTTCGCGCTTATTGTTGGGCGCCTGTACCGGCCTTTTCGCGGGCTTGCCCGCTAAAAGGccg                                                                                                                         <  2:173367/148‑2 (MQ=255)
              gAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                         <  5:285730/148‑1 (MQ=255)
              gAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGACTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAg                                                                                                                        >  1:88386/1‑149 (MQ=255)
               aaCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                           >  5:160588/1‑145 (MQ=255)
               aaCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                           <  6:160588/145‑1 (MQ=255)
               aaCATTGGCCATGCGCTGGTGGCGCATGCATTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAg                                                                                                                        >  5:371677/1‑148 (MQ=255)
                   ttGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAg                                                                                                                   <  4:215736/149‑1 (MQ=255)
                       ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcgg                                                                                                               >  3:258235/1‑149 (MQ=255)
                        cATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGGTAAGAGGCCAGCACAggcggc                                                                                                              <  2:190592/149‑1 (MQ=255)
                                     cgcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                                 <  1:41710/149‑1 (MQ=255)
                                                    tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAgg                                                                                                                  >  3:265549/1‑117 (MQ=255)
                                                    tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAgg                                                                                                                  <  4:265549/117‑1 (MQ=255)
                                                       gCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAc                                                                               >  3:123771/1‑149 (MQ=255)
                                                       gCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAc                                                                               >  4:229834/1‑149 (MQ=255)
                                                        cTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACAcc                                                                                        >  2:345583/1‑139 (MQ=255)
                                                        cTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACAcc                                                                                        <  1:345583/139‑1 (MQ=255)
                                                         ttCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCg                                                                             >  3:324772/1‑149 (MQ=255)
                                                                        gCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                               <  2:347674/148‑1 (MQ=255)
                                                                         ccGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGAACACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                             >  6:82503/1‑149 (MQ=255)
                                                                          cGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGc                                                            <  4:123771/149‑1 (MQ=255)
                                                                           gagaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGt                                                                          >  7:97663/1‑134 (MQ=255)
                                                                           gagaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGt                                                                          <  8:97663/134‑1 (MQ=255)
                                                                           gagaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCa                                                           <  4:324772/149‑1 (MQ=255)
                                                                            agaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCGTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAg                                                          <  6:371677/149‑1 (MQ=255)
                                                                               tGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctg                                                       <  4:258235/149‑1 (MQ=255)
                                                                                gAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                      >  1:249611/1‑149 (MQ=255)
                                                                                  aGGCGCTGAACGTGGCGGCTTCGCGCTTATTGTAGGTCGCCTGTACAGGCCTTTTCGCGGGCAAGCACGCGAAGAGGACAGCACAGGCGTCCGATCACTCGGCCGGAAACACCAGGGTAATACGGGTCGGCCCCAGGGGGCCgctgatt                                                    >  4:265313/1‑146 (MQ=255)
                                                                                   ggCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAgg                                                                                               <  7:315833/105‑1 (MQ=255)
                                                                                   ggCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAgg                                                                                               >  8:315833/1‑105 (MQ=255)
                                                                                      gcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                >  8:100192/1‑149 (MQ=255)
                                                                                          aTCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgct                                                     <  3:349164/140‑1 (MQ=255)
                                                                                          aTCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgct                                                     >  4:349164/1‑140 (MQ=255)
                                                                                          aTCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                            <  5:330579/149‑1 (MQ=255)
                                                                                              tGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATGACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                        >  5:166711/1‑149 (MQ=255)
                                                                                                  ggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACccgcccgc                                    <  5:82503/149‑1 (MQ=255)
                                                                                                         cgcTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCaa                             <  3:86779/149‑1 (MQ=255)
                                                                                                           cTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                             <  7:217029/115‑1 (MQ=255)
                                                                                                           cTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                             >  8:217029/1‑115 (MQ=255)
                                                                                                           cTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                           >  7:203568/1‑149 (MQ=255)
                                                                                                              aTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCa                        >  3:92653/1‑149 (MQ=255)
                                                                                                              aTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCa                        <  2:249611/149‑1 (MQ=255)
                                                                                                              aTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCa                        >  3:97825/1‑149 (MQ=255)
                                                                                                                      tCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGAc                >  4:13751/1‑149 (MQ=255)
                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa            <  2:109313/149‑1 (MQ=255)
                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa            <  2:88386/149‑1 (MQ=255)
                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa            >  5:102022/1‑149 (MQ=255)
                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa            >  1:335367/1‑149 (MQ=255)
                                                                                                                               ccGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc       <  6:166711/149‑1 (MQ=255)
                                                                                                                                tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCg      >  7:288673/2‑149 (MQ=39)
                                                                                                                                cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGa         >  1:274470/1‑146 (MQ=255)
                                                                                                                                cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGa         <  2:274470/146‑1 (MQ=255)
                                                                                                                                  gCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc     >  3:155348/1‑148 (MQ=255)
                                                                                                                                    ctcttCGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaa  >  4:171635/4‑149 (MQ=255)
                                                                                                                                    cTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCcgcccgcc                                   >  2:127974/1‑116 (MQ=255)
                                                                                                                                    cTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCcgcccgcc                                   <  1:127974/116‑1 (MQ=255)
                                                                                                                                      cttCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaa  >  6:316541/2‑147 (MQ=255)
                                                                                                                                      cttCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaa  <  5:316541/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
GGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAA  >  NC_002947/1637181‑1637461

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: ‑