Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,549 T→G 6.7% intergenic (+38/+44) PP_4704 → / ← PP_4705 hypothetical protein/glyoxalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,5490TG 6.7% 174.9 / 4.5 60intergenic (+38/+44)PP_4704/PP_4705hypothetical protein/glyoxalase
Reads supporting (aligned to +/- strand):  ref base T (33/23);  new base G (2/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGTTGAAC  >  NC_002947/5350401‑5350698
                                                                                                                                                      |                                                                                                                                                     
gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTa                                                                                                                                                      >  2:1323889/1‑150 (MQ=255)
 ggtgacggcgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCTCCCACAGGCCTGCGCACTAc                                                                                                                                                     >  2:731697‑M1/124‑150 (MQ=255)
      aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                                                >  2:861421/1‑150 (MQ=255)
       aGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                                               <  1:125165/150‑1 (MQ=255)
         cgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGa                                                                                                                                             >  2:1249711‑M1/116‑150 (MQ=255)
           cAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGc                                                                                                                                           <  1:659021/150‑1 (MQ=255)
            aaGCTTGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCg                                                                                                                                          <  1:340934/150‑1 (MQ=255)
                 ggAAGTCGACTAGGTGTACTTCATTGATACTAGCCTGGGGCAGCCTCAAAGAGATGACCGGGTAGTGTCACCGATGTGCCTGGATACAGAGTCATCCCCTCTGGCGGGACAGCCCCGCGAGGCCTGCGCGCTATCTGTGGGAGGGGGCGt                                                                                                                                     >  2:681911/1‑150 (MQ=255)
                 gcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGt                                                                                                                                     >  2:782840‑M1/108‑150 (MQ=255)
                               tgtACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                                                       <  1:1323889/150‑1 (MQ=255)
                                 tACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                                                     >  1:762406/1‑150 (MQ=255)
                                  ccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                                                >  1:363118‑M1/91‑122 (MQ=255)
                                  ccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                                                <  2:363118‑M1/32‑1 (MQ=255)
                                     gccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                                               <  1:451920‑M1/33‑1 (MQ=255)
                                     gccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                                               >  2:451920‑M1/88‑120 (MQ=255)
                                       aTCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACgg                                                                                                               >  2:1396755/1‑150 (MQ=255)
                                        tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACggg                                                                                                              >  1:1344005/1‑150 (MQ=255)
                                        tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACggg                                                                                                              <  2:606593/150‑1 (MQ=255)
                                                           gCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCCCTGcc                                                                                           >  1:141852/1‑150 (MQ=255)
                                                             cTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcccc                                                                                         >  2:1254468/1‑150 (MQ=255)
                                                              tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTTGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCg                                                                                        >  1:336708/1‑150 (MQ=255)
                                                               cAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGa                                                                                       <  1:1254468/150‑1 (MQ=255)
                                                                 acccggcggtgtcggggcagcggttgttttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                                                     >  2:903558‑M1/60‑118 (MQ=255)
                                                                      aTGGCCGGGTTGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACAccgcc                                                                                >  1:1265637/1‑150 (MQ=255)
                                                                          ccGGGTGGTGACACCGATGCACCTGGATGCCGAGTCAGCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                                            <  2:336708/150‑1 (MQ=255)
                                                                                     cacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCAAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                                                 <  2:346893/150‑1 (MQ=255)
                                                                                      acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                                                                >  1:33114/1‑150 (MQ=255)
                                                                                          aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                                            <  1:8298/150‑1 (MQ=255)
                                                                                          aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                                            <  2:1344005/150‑1 (MQ=255)
                                                                                                cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                                      >  2:618929/1‑150 (MQ=255)
                                                                                                cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                                      >  1:992257/1‑150 (MQ=255)
                                                                                                cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                                      <  2:762406/150‑1 (MQ=255)
                                                                                                        cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccg                                              >  2:124273/1‑150 (MQ=255)
                                                                                                        cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccg                                              <  2:1210988/150‑1 (MQ=255)
                                                                                                          tggcctcttcgcgggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                                            >  1:903558‑M1/19‑118 (MQ=255)
                                                                                                               tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCt                                                                                              >  1:1149929/1‑95 (MQ=255)
                                                                                                               tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCt                                                                                              <  2:1149929/95‑1 (MQ=255)
                                                                                                                 cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca                                     >  2:1215650/1‑150 (MQ=255)
                                                                                                                 cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca                                     >  1:217633/1‑150 (MQ=255)
                                                                                                                  ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagcag                                    >  1:125291/1‑150 (MQ=255)
                                                                                                                   ctctGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAgg                                   >  1:58632/1‑150 (MQ=255)
                                                                                                                     cgggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCCTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTg                                 <  2:731720‑M1/143‑1 (MQ=255)
                                                                                                                      gggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                                           >  2:705306‑M1/7‑139 (MQ=255)
                                                                                                                      gggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                                           <  1:705306‑M1/133‑1 (MQ=255)
                                                                                                                         cAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                                                   >  2:1031191/1‑112 (MQ=255)
                                                                                                                         cAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                                                   <  1:1031191/112‑1 (MQ=255)
                                                                                                                          aGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGc                            >  2:1211884/1‑150 (MQ=255)
                                                                                                                          aGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGc                            >  1:940848/1‑150 (MQ=255)
                                                                                                                            GCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCtg                          <  1:851214‑M1/150‑1 (MQ=255)
                                                                                                                                 CCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCtgcgtgc                     >  1:1396755‑M1/1‑150 (MQ=255)
                                                                                                                                    ccACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGcc                  >  2:1081903/1‑150 (MQ=255)
                                                                                                                                      acaGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGt                >  1:711041/1‑150 (MQ=255)
                                                                                                                                        aGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCa              <  2:58632/150‑1 (MQ=255)
                                                                                                                                           ccTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCa           <  1:1215650/150‑1 (MQ=255)
                                                                                                                                                 cgcACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGt       <  1:124273/150‑1 (MQ=255)
                                                                                                                                                  gcACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGtt      <  2:940848/150‑1 (MQ=255)
                                                                                                                                                  gcACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGtt      <  1:861421/150‑1 (MQ=255)
                                                                                                                                                    aCTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                                                 <  1:1161723/89‑1 (MQ=255)
                                                                                                                                                    aCTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                                                 >  2:1161723/1‑89 (MQ=255)
                                                                                                                                                    aCTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCAAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGTTGa    >  1:771688/1‑150 (MQ=255)
                                                                                                                                                      tACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGTTGAAc  >  1:66760/1‑150 (MQ=255)
                                                                                                                                                      tACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGTTGAAc  <  2:217633/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCATGGTTGAAC  >  NC_002947/5350401‑5350698

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