Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,604,002 A→T 7.1% L232Q (CTG→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,0020AT 7.1% 202.3 / 6.7 70L232Q (CTG→CAG) PP_3179LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (27/38);  new base T (2/3);  total (29/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 = 8.37e-01

CGTTCCATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACG  >  NC_002947/3603856‑3604151
                                                                                                                                                  |                                                                                                                                                     
tgtgccATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCcagca                                                                                                                                                    <  1:191510/146‑1 (MQ=255)
  ttCCATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGCTGAGCGGACGACGGACCCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCCgg                                                                                                                                                  >  2:53679/1‑150 (MQ=255)
   tCCATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCcag                                                                                                                                                      <  2:1604019/145‑1 (MQ=255)
   tCCATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCcag                                                                                                                                                      >  1:1604019/1‑145 (MQ=255)
   tCCATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAggg                                                                                                                                                 >  1:1132089/1‑150 (MQ=255)
     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:1094862/1‑150 (MQ=255)
                          ccTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGc                                                                                                                          >  1:1821019/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                                                                                                                                      >  1:1473052/1‑128 (MQ=255)
                                    gCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTgc                                                                                                                                      <  2:1473052/128‑1 (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                                                                                                                                           >  1:223033/1‑110 (MQ=255)
                                                tAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCACACCCTGCTGGGTCTggt                                                                                                                                           <  2:223033/111‑2 (MQ=255)
                                                          ggggCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGca                                                                                          <  2:1514793/150‑1 (MQ=255)
                                                          ggggCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGca                                                                                          <  1:742165/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                                                                                                                      <  2:899848/109‑1 (MQ=255)
                                                                       aCGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGa                                                                                                                      >  1:899848/1‑109 (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                                          >  2:887597/1‑150 (MQ=255)
                                                                                                          cGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGt                                          >  1:1877540/1‑150 (MQ=255)
                                                                                                          cGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGt                                          >  1:1748706/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)
                                                                                                                                                  |                                                                                                                                                     
CGTTCCATTTCATGTACAGCCCAAGCCCTGGGGGTGGCTGGTCATCGATAGGGCGGATGGGGCAAAACTGCACGCGTTCGCGCGGGGTTGAGCGGACGAAGGACGCCGGCACCAGGGCCACGCCAAAGCCGCAGGCCACCAGACCCAGCAGGGTGTGGATCTGCGCCGCCTCCTGCACGACCTTGGGGTAGAACCCGGCCGACTCGCACAGCGCAACCAACTGTGCGTGGTAGCCCCCGCCCAGCGCTTGCGGGGTGAACACGAAGTCCTCATCGGCAAACTCGCGCAAGTCGACG  >  NC_002947/3603856‑3604151

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