Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,604,010 G→C 6.9% H229Q (CAC→CAG)  PP_3179 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,604,0100GC 6.9% 202.3 / 5.0 67H229Q (CAC→CAG) PP_3179LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (24/38);  new base C (2/3);  total (26/41)
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.84e-01

CATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAA  >  NC_002947/3603861‑3604159
                                                                                                                                                     |                                                                                                                                                     
cATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGgtg                                                                                                                                                       <  1:1817487/150‑1 (MQ=255)
  tttCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGtgtg                                                                                                                                                     <  1:1001716/150‑1 (MQ=255)
  tttCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGAGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCACGCCACCAGACCCAGCAGGgtgag                                                                                                                                                     >  2:638973/1‑148 (MQ=255)
     cATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCACACCCTGCTGGGTCTggtg                                                                                                                                                  <  2:12676/150‑3 (MQ=255)
       tGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCt                                                                                                                                                <  2:662881/150‑1 (MQ=255)
         tACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTgc                                                                                                                                              >  2:1470097/1‑150 (MQ=255)
             gCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGcgccg                                                                                                                                          >  2:11133/1‑150 (MQ=255)
              cccAAGCCCTGGGGGAGGCAGGTCATCGATTGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCAGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGcggcgc                                                                                                                                         >  2:1115298/1‑150 (MQ=255)
                    cccTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTg                                                                                                                                   >  2:127775/1‑150 (MQ=255)
                     ccTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGc                                                                                                                                  >  1:1821019/1‑150 (MQ=255)
                     ccTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGc                                                                                                                                  >  1:1094862/1‑150 (MQ=255)
                       tGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCAc                                                                                                                                >  2:1683417/1‑150 (MQ=255)
                         ggggTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGa                                                                                                                              <  2:911876/150‑1 (MQ=255)
                            ggggCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCt                                                                                                                           >  2:1920465/3‑150 (MQ=255)
                              ggCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCCCCCGAGCCACGCCAAAGCCGCAGGCTTCCGGACCGATCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTg                                                                                                                         <  1:1470097/150‑1 (MQ=255)
                              ggCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTg                                                                                                                         <  1:1623693/150‑1 (MQ=255)
                              ggCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTg                                                                                                                         <  1:1045540/150‑1 (MQ=255)
                              ggCTGGTCATCGATAGGGCGGATGGGGAAACCAAACACGATTTCGCGCGGGGTTTAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTg                                                                                                                         <  1:888092/150‑1 (MQ=255)
                               gCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTgc                                                                                                                                              <  2:1473052/128‑1 (MQ=255)
                               gCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTgc                                                                                                                                              >  1:1473052/1‑128 (MQ=255)
                                ctgggtaTCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGAGCCGCCTCCTGCACGACCTTggg                                                                                                                       <  2:943769/144‑1 (MQ=255)
                                cTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTggg                                                                                                                       <  2:513709/150‑1 (MQ=255)
                                cTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTggg                                                                                                                       <  1:1426297/150‑1 (MQ=255)
                                 tGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTgggg                                                                                                                      >  2:1795615/1‑150 (MQ=255)
                                   gTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTa                                                                                                                    >  1:1113314/1‑150 (MQ=255)
                                    tCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAg                                                                                                                   <  1:732618/150‑1 (MQ=255)
                                       tCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAAc                                                                                                                <  1:1683417/150‑1 (MQ=255)
                                        cGAGAGGGCGGAAGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAAcc                                                                                                               <  2:423696/150‑1 (MQ=255)
                                           tAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCACACCCTGCTGGGTCTggt                                                                                                                                                   <  2:223033/111‑2 (MQ=255)
                                           tAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCACACCCTGCTGGGTCTggt                                                                                                                                                   >  1:223033/1‑110 (MQ=255)
                                                     ggggCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGca                                                                                                  <  1:742165/150‑1 (MQ=255)
                                                     ggggCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGca                                                                                                  <  2:1514793/150‑1 (MQ=255)
                                                      gggCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGcac                                                                                                 <  1:1157077/150‑1 (MQ=255)
                                                        gCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCTGCTGGGTCTggtgg                                                                                                                                                 >  2:251013/1‑97 (MQ=255)
                                                        gCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCTGCTGGGTCTggtgg                                                                                                                                                 <  1:251013/100‑4 (MQ=255)
                                                         cAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAgc                                                                                              <  2:1132089/150‑1 (MQ=255)
                                                             aCTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGcaa                                                                                          <  1:1892804/150‑1 (MQ=255)
                                                                  aCGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGa                                                                                                                              >  1:899848/1‑109 (MQ=255)
                                                                  aCGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGa                                                                                                                              <  2:899848/109‑1 (MQ=255)
                                                                      gTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTgggg                                                                                                                      <  2:641997/113‑1 (MQ=255)
                                                                      gTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTgggg                                                                                                                      >  1:641997/1‑113 (MQ=255)
                                                                         cgcgcgGGGTGGAGCGGACGAAGGACGCCGGCACCATGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGGGTGGATCTGCGCCGCCTCCTGCACGCCCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGt                                                                              <  2:23917/150‑1 (MQ=255)
                                                                             cgGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTa                                                                          <  2:1821019/150‑1 (MQ=255)
                                                                               gggTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGc                                                                        <  1:1655163/150‑1 (MQ=255)
                                                                                     aGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGc                                                                  <  1:1827254/150‑1 (MQ=255)
                                                                                               ggACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAgc                                                                                              <  1:1904691/112‑1 (MQ=255)
                                                                                               ggACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAgc                                                                                              >  2:1904691/1‑112 (MQ=255)
                                                                                                 aCGCGGGCACCAGGGCAACGCCAAGGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTTGGGTAGAGCCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCg                                                      <  1:638973/150‑1 (MQ=255)
                                                                                                 aCGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCg                                                      >  1:848466/1‑150 (MQ=255)
                                                                                                    cgggCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCgggg                                                   <  2:37393/148‑1 (MQ=255)
                                                                                                    ccGGTACCAGGGCGACGCCATAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCgggg                                                   <  2:37304/150‑1 (MQ=255)
                                                                                                     cGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGt                                                  >  1:1877540/1‑150 (MQ=255)
                                                                                                     cGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGt                                                  >  1:1748706/1‑150 (MQ=255)
                                                                                                     cGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGt                                                  >  2:887597/1‑150 (MQ=255)
                                                                                                           cAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAacac                                            >  1:1347648/1‑150 (MQ=255)
                                                                                                             gggCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGa                                          <  2:1007296/150‑1 (MQ=255)
                                                                                                              ggCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGaa                                         <  1:862642/150‑1 (MQ=255)
                                                                                                                ccACGCCAAAGCCGCAGGACTGGAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGt                                       <  1:1533626/150‑1 (MQ=255)
                                                                                                                 cACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCCCCCAGCGCTTGCGGGGTGAACACGAAGTc                                      >  2:723812/1‑150 (MQ=255)
                                                                                                                  aCGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCTACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTcc                                     <  2:1094862/150‑1 (MQ=255)
                                                                                                                  aCGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGcaaccaa                                                                                      >  1:602723/1‑101 (MQ=255)
                                                                                                                  aCGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGcaaccaa                                                                                      <  2:602723/101‑1 (MQ=255)
                                                                                                                        aaGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGGCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTAGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCg                               <  2:170792/150‑1 (MQ=255)
                                                                                                                            cGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCCCGCACCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCaa                           >  1:423290/1‑150 (MQ=255)
                                                                                                                                ggCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTc                       <  2:1579337/150‑1 (MQ=255)
                                                                                                                                   caccaGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTcgcg                    >  1:1283911/1‑150 (MQ=255)
                                                                                                                                      caGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCaa                 <  1:11133/150‑1 (MQ=255)
                                                                                                                                           ccAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGa            >  1:1908097/1‑150 (MQ=255)
                                                                                                                                             agcagGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCCCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACg          <  1:723812/150‑1 (MQ=255)
                                                                                                                                                  ggTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTc     <  2:1283911/150‑1 (MQ=255)
                                                                                                                                                    tgtgGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGa   >  2:1526896/1‑150 (MQ=255)
                                                                                                                                                     gtgGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCGGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGaa  >  2:1163619/1‑150 (MQ=255)
                                                                                                                                                     gtgGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGaa  >  2:1240741/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACGGCCTCGAA  >  NC_002947/3603861‑3604159

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