Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,909 C→T 15.0% S338N (AGC→AAC)  PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,9090CT15.0% 155.1 / 18.5 68S338N (AGC→AAC) PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base T (4/6);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

TAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCC  >  NC_002947/3865762‑3866058
                                                                                                                                                   |                                                                                                                                                     
tAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGttc                                                                                                                                                     <  1:2069934/150‑1 (MQ=255)
      cGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGCCCGATCATTCCTGTGGGAGTCGCGCGAGACGGCCTTGGCGGGCACTCTCAGGTGCTCTCgcgc                                                                                                                                               >  1:34295/1‑150 (MQ=255)
         gTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCAcc                                                                                                                                            >  2:1085153/1‑150 (MQ=255)
               ccGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGc                                                                                                                                      <  1:993071/150‑1 (MQ=255)
                    aaCACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCaaa                                                                                                                                 <  1:1788836/150‑1 (MQ=255)
                    aaCACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCaaa                                                                                                                                 <  2:1204913/150‑1 (MQ=255)
                     acacCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAg                                                                                                                                <  1:1034411/150‑1 (MQ=255)
                     acacCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAg                                                                                                                                <  1:1824047/150‑1 (MQ=255)
                                  cAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCctgct                                                                                                                   <  1:1745812/150‑1 (MQ=255)
                                         aTGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCAtt                                                                                                            >  2:437179/1‑150 (MQ=255)
                                              ccGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCa                                                                                                       >  1:1808168/1‑150 (MQ=255)
                                                       cTGCTTCGCAGGAAAACCCGCTGCAGCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGCCGGGCTCGCTCAGGTGGTTTCGGGGACCATCAGCTCGAAGCCCAGAGGCTGCTGCGGATTCGTCACCCGCTTcc                                                                                              <  2:34295/150‑1 (MQ=255)
                                                       cTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTcc                                                                                              >  2:778860/1‑150 (MQ=255)
                                                         gCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCa                                                                                            >  2:876890/1‑150 (MQ=255)
                                                           ttCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATc                                                                                          <  1:862042/150‑1 (MQ=255)
                                                             cGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCa                                                                                        <  2:760034/150‑1 (MQ=255)
                                                              gCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCGCTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAg                                                                                       >  2:1324699/1‑150 (MQ=255)
                                                                 ggAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCaa                                                                                    <  2:1436661/150‑1 (MQ=255)
                                                                 ggAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCaa                                                                                    >  1:2235171/1‑150 (MQ=255)
                                                                   aaaaCCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACg                                                                                  >  1:1280030/1‑150 (MQ=255)
                                                                        ccGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCcagc                                                                             >  2:254452/1‑150 (MQ=255)
                                                                         cGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCcagca                                                                            <  2:1853082/150‑1 (MQ=255)
                                                                         cGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCCGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCcagca                                                                            <  1:219645/150‑1 (MQ=255)
                                                                             gCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAgcgc                                                                        <  2:1808168/150‑1 (MQ=255)
                                                                                   tCGTTTCGGGCCAATCATTACTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCGTCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGc                                                                  <  2:1582971/150‑1 (MQ=255)
                                                                                      tttCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCcgc                                                               >  2:208444/1‑150 (MQ=255)
                                                                                          gggCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGcgccc                                                           >  2:2178091/1‑150 (MQ=255)
                                                                                          gggCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGcgccc                                                           >  1:972177/1‑150 (MQ=255)
                                                                                             ccAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTc                                                                                                           <  1:435137/99‑1 (MQ=255)
                                                                                             ccAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTc                                                                                                           >  2:435137/1‑99 (MQ=255)
                                                                                             ccAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGcgccccgc                                                        >  2:2137433/1‑150 (MQ=255)
                                                                                                 tCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCa                                                                                                                                              >  1:359713/1‑60 (MQ=255)
                                                                                                 tCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCa                                                                                                                                              <  2:359713/60‑1 (MQ=255)
                                                                                                   aTTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCtgctgc                                                                                                                 >  1:1100578/1‑87 (MQ=255)
                                                                                                   aTTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCtgctgc                                                                                                                 <  2:1100578/87‑1 (MQ=255)
                                                                                                      ccTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCa                                               <  1:437179/150‑1 (MQ=255)
                                                                                                       cTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAg                                              >  2:204682/1‑150 (MQ=255)
                                                                                                        tgcgggAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGc                                             >  1:1242489/4‑150 (MQ=255)
                                                                                                          tgGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTc                                                                                               >  1:2121353/1‑98 (MQ=255)
                                                                                                          tgGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTc                                                                                               <  2:2121353/98‑1 (MQ=255)
                                                                                                           gggAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCACCGCCCCACCGCAGCCcg                                          >  2:1306778/1‑150 (MQ=255)
                                                                                                            ggAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCctg                                                                                                                     <  2:2230168/74‑1 (MQ=255)
                                                                                                            ggAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCctg                                                                                                                     >  1:2230168/1‑74 (MQ=255)
                                                                                                            ggAGTAGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGGCCTCTCGGGCAGGATCAGCACAAAGCCGAGATGCAGCTGCGCATTGGTAACGCGCTGGGCATGCAGCAAGGCCAGGAGCGCGTGTGCCGCGCGGCGAGCCAACGCAgagaga                                         >  2:936105/1‑145 (MQ=255)
                                                                                                               gTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgt                                      >  1:492440/1‑150 (MQ=255)
                                                                                                               gTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCGTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgt                                      >  2:1632820/1‑150 (MQ=255)
                                                                                                               gTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgt                                      >  1:1791780/1‑150 (MQ=255)
                                                                                                                   cgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGc                                  <  1:2137433/150‑1 (MQ=255)
                                                                                                                       agaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCCCCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAAt                              <  2:1772398/150‑1 (MQ=255)
                                                                                                                        gaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTc                                                                                               <  1:1017826/84‑1 (MQ=255)
                                                                                                                        gaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTc                                                                                               >  2:1017826/1‑84 (MQ=255)
                                                                                                                          aGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGa                           <  1:1085153/150‑1 (MQ=255)
                                                                                                                           ggCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                          <  2:492440/150‑1 (MQ=255)
                                                                                                                           ggCCTTGGCGGGCACTCTCAGGTGTTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                          <  1:204682/150‑1 (MQ=255)
                                                                                                                            gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCgg                                       <  2:1187462/136‑1 (MQ=255)
                                                                                                                            gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCgg                                       >  1:1187462/1‑136 (MQ=255)
                                                                                                                            gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAgg                         <  2:1280030/150‑1 (MQ=255)
                                                                                                                              cTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTc                       >  1:883681/1‑150 (MQ=255)
                                                                                                                                   cGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCtg                                                                     >  1:69585/1‑99 (MQ=255)
                                                                                                                                   cGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCtg                                                                     <  2:69585/99‑1 (MQ=255)
                                                                                                                                      gCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTgg               <  1:778860/150‑1 (MQ=255)
                                                                                                                                        aCTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCagcagc                                                                          >  1:843988/1‑89 (MQ=255)
                                                                                                                                        aCTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCagcagc                                                                          <  2:843988/89‑1 (MQ=255)
                                                                                                                                          tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGACTGGCGGTCAGCCGTAGCAcc           >  2:660839/1‑150 (MQ=255)
                                                                                                                                          tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTgg               >  1:330521/1‑146 (MQ=255)
                                                                                                                                          tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTgg               <  2:330521/146‑1 (MQ=255)
                                                                                                                                          tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAcc           >  2:1193587/1‑150 (MQ=255)
                                                                                                                                          tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGGCAGCCGTGGCAcc           >  2:1302007/1‑150 (MQ=255)
                                                                                                                                                gtgtTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATgc     >  1:379903/1‑150 (MQ=255)
                                                                                                                                                gtgtTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATgc     >  1:1473609/1‑150 (MQ=255)
                                                                                                                                                gTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATgc     >  1:1591757/1‑150 (MQ=255)
                                                                                                                                                   ctctCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGcc  <  2:2235171/150‑1 (MQ=255)
                                                                                                                                                   ctctCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGcc  <  1:254452/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCC  >  NC_002947/3865762‑3866058

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