Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,989,557 C→G 5.0% G677G (GGC→GGG)  oplA → 5‑oxoprolinase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,989,5570CG 5.0% 317.7 / 6.4 99G677G (GGC→GGG) oplA5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base C (51/43);  new base G (2/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGT  >  NC_002947/3989408‑3989705
                                                                                                                                                     |                                                                                                                                                    
gcgcAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGc                                                                                                                                                      <  1:2185006/150‑1 (MQ=255)
 cgcAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCt                                                                                                                                                     <  1:1608236/150‑1 (MQ=255)
  gcAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTg                                                                                                                                                    <  2:420423/150‑1 (MQ=255)
     aTCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCa                                                                                                                                                 >  2:1630523/1‑150 (MQ=255)
        cgtccgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGc                                                                                                                                              >  1:680126/1‑150 (MQ=255)
          tccgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTc                                                                                                                                            >  2:844430/1‑150 (MQ=255)
            cgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGc                                                                                                                                          <  2:1684432/150‑1 (MQ=255)
            cgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGc                                                                                                                                          <  1:1630523/150‑1 (MQ=255)
              cccGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCt                                                                                                                                        <  1:946023/150‑1 (MQ=255)
              cccGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCt                                                                                                                                        >  2:2207821/1‑150 (MQ=255)
               ccGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGAAGGCTCGCCTg                                                                                                                                       >  2:769433/1‑150 (MQ=255)
                 gATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGa                                                                                                                                     >  2:494840/1‑150 (MQ=255)
                 gATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAAACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGa                                                                                                                                     <  2:1599361/150‑1 (MQ=255)
                      tGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCtt                                                                                                                                <  1:1665590/150‑1 (MQ=255)
                            gTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCa                                                                                                                          <  1:2207821/150‑1 (MQ=255)
                            gTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCa                                                                                                                          >  1:2196925/1‑150 (MQ=255)
                             tGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCAc                                                                                                                         >  2:420056/1‑150 (MQ=255)
                                 gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCATACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGATCGCCTGGAGCCTGCGGCCACCTAc                                                                                                                     >  2:1261098/1‑150 (MQ=255)
                                 gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCg                                                                                                                              >  2:2260990/1‑141 (MQ=255)
                                 gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCg                                                                                                                              <  1:2260990/141‑1 (MQ=255)
                                 gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACAACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAgga                                                                                                                                              <  1:1504750/125‑2 (MQ=255)
                                 gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACAACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAgga                                                                                                                                              >  2:1504750/1‑124 (MQ=255)
                                  aCGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTAcc                                                                                                                    <  1:608658/150‑1 (MQ=255)
                                    ggTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCgg                                                                                                                  <  2:1551531/150‑1 (MQ=255)
                                    ggTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCgg                                                                                                                  <  2:2175797/150‑1 (MQ=255)
                                     gTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCggg                                                                                                                 <  2:236714/150‑1 (MQ=255)
                                     gTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCggg                                                                                                                 >  1:2223092/1‑150 (MQ=255)
                                      tGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCgggg                                                                                                                >  2:2014004/1‑149 (MQ=255)
                                        aaTGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGa                                                                                                              >  2:611774/1‑150 (MQ=255)
                                            gCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACAcg                                                                                                          >  1:1369612/1‑150 (MQ=255)
                                              aGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACAcgcg                                                                                                        <  1:611774/150‑1 (MQ=255)
                                                 ccACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAggacaa                                                                                                                                           <  2:145121/112‑5 (MQ=255)
                                                 ccACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAggacaa                                                                                                                                           >  1:145121/1‑108 (MQ=255)
                                                 ccACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCgggg                                                                                                     <  2:2360896/150‑1 (MQ=255)
                                                  cacaCCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCggggg                                                                                                    >  1:1762703/1‑150 (MQ=255)
                                                   acacCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGt                                                                                                   >  1:2200258/1‑150 (MQ=255)
                                                    cacCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGtt                                                                                                  <  2:1589735/150‑1 (MQ=255)
                                                       tggTGTGCACAGGCAACCACCTGCTGTCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGAt                                                                                               <  2:1960200/149‑1 (MQ=255)
                                                       cGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGAt                                                                                               >  1:598528/1‑150 (MQ=255)
                                                       cGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGAt                                                                                               >  1:598790/1‑150 (MQ=255)
                                                         gtgtTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCt                                                                                             >  2:503028/1‑150 (MQ=255)
                                                            ttCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGc                                                                                          <  2:1387099/150‑1 (MQ=255)
                                                             tCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                         <  2:569881/150‑1 (MQ=255)
                                                              cAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTg                                                                                           >  1:676762/1‑147 (MQ=255)
                                                              cAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTg                                                                                           <  2:676762/147‑1 (MQ=255)
                                                                aCGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTc                                                                                      >  2:2305165/1‑150 (MQ=255)
                                                                 cGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCg                                                                                     <  2:2210097/150‑1 (MQ=255)
                                                                  ggCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCCACGTTGTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCgg                                                                                    >  1:23554/1‑150 (MQ=255)
                                                                        caccTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGaca                                                                                                            >  2:1733813/1‑120 (MQ=255)
                                                                        caccTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGaca                                                                                                            <  1:1733813/120‑1 (MQ=255)
                                                                          ccTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                         <  1:720756/137‑1 (MQ=255)
                                                                          ccTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                         >  2:720756/1‑137 (MQ=255)
                                                                               tGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCa                                                                       <  1:844430/150‑1 (MQ=255)
                                                                                  cTGGGCAGACCGTGCATGGCCCCTGTGTGCTCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATgg                                                                    <  2:117134/150‑1 (MQ=255)
                                                                                   tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                                   <  2:793375/150‑1 (MQ=255)
                                                                                   tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                                   >  1:1035777/1‑150 (MQ=255)
                                                                                         gACCGTGCATGGCCCCTGTGTGATCGAGGAAGAGACGACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGAt                                                             >  1:2019040/1‑150 (MQ=255)
                                                                                           ccGTGCATGGCCCCTGTGCGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATaa                                                           <  1:1719095/150‑1 (MQ=255)
                                                                                              tGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                         <  1:2369327/117‑1 (MQ=255)
                                                                                              tGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                         >  2:2369327/1‑117 (MQ=255)
                                                                                                 aTGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCTTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGccc                                                     >  1:1561741/1‑150 (MQ=255)
                                                                                                    gCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAAc                                                          <  2:777685/142‑1 (MQ=255)
                                                                                                    gCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAAc                                                          >  1:777685/1‑142 (MQ=255)
                                                                                                     ccccTGTGTGATCGAGGAAGACACCACCCTTTTTGTCGTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCa                                                 <  2:2019040/150‑1 (MQ=255)
                                                                                                     ccccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCa                                                 >  1:1602950/1‑150 (MQ=255)
                                                                                                       ccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAgg                                               >  2:1484815/1‑149 (MQ=255)
                                                                                                        cTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGCTGACACGCGGGGGCTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGtt                                              >  2:1303365/1‑150 (MQ=255)
                                                                                                         tgtgtgATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGttt                                             <  2:2200258/150‑1 (MQ=255)
                                                                                                                  gAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCgg                                                                                                       <  1:717906/83‑1 (MQ=255)
                                                                                                                  gAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCgg                                                                                                       >  2:717906/1‑83 (MQ=255)
                                                                                                                   aGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGAAGGGGGTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGccc                                   >  1:747775/1‑149 (MQ=255)
                                                                                                                     gAAGACACCAGCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCggg                                 >  2:1927596/1‑150 (MQ=255)
                                                                                                                     gAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCggg                                 >  2:1235581/1‑150 (MQ=255)
                                                                                                                     gAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGATGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCggg                                 >  2:1689376/1‑150 (MQ=255)
                                                                                                                      aaGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGc                                <  1:1814246/150‑1 (MQ=255)
                                                                                                                      aaGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGc                                >  2:533091/1‑150 (MQ=255)
                                                                                                                      aaGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGc                                <  1:1014977/150‑1 (MQ=255)
                                                                                                                      aaGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGc                                <  2:2223092/150‑1 (MQ=255)
                                                                                                                        gACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCgg                              >  2:1540867/1‑150 (MQ=255)
                                                                                                                           accaccAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTc                           >  1:387552/1‑150 (MQ=255)
                                                                                                                            ccacctATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCg                          <  2:1567934/150‑1 (MQ=255)
                                                                                                                             caccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGt                         >  1:636232/1‑150 (MQ=255)
                                                                                                                             caccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGt                         >  1:1054952/1‑150 (MQ=255)
                                                                                                                             caccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGt                         >  2:679557/1‑150 (MQ=255)
                                                                                                                              accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGTCACCTACCGGGTGACACGCGGGGGTTGATCATGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTa                        >  1:1000982/1‑150 (MQ=255)
                                                                                                                              accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTa                        <  1:420056/150‑1 (MQ=255)
                                                                                                                              accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTa                        >  1:159725/1‑150 (MQ=255)
                                                                                                                              accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTa                        >  1:1232031/1‑150 (MQ=255)
                                                                                                                                        tCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAag              <  1:494840/150‑1 (MQ=255)
                                                                                                                                         ccTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAaga             <  2:1369612/150‑1 (MQ=255)
                                                                                                                                            gCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAgct          >  1:1885049/1‑150 (MQ=255)
                                                                                                                                               aCCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGa                                                              >  1:34961/1‑95 (MQ=255)
                                                                                                                                               aCCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGa                                                              <  2:34961/95‑1 (MQ=255)
                                                                                                                                               aCCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAgctgct       >  1:2155045/1‑150 (MQ=255)
                                                                                                                                                 ccGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTg     >  2:1440583/1‑150 (MQ=255)
                                                                                                                                                  cGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGa    >  2:654494/1‑150 (MQ=255)
                                                                                                                                                  cGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGa    <  1:1440583/150‑1 (MQ=255)
                                                                                                                                                   ggCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAg   <  1:654494/150‑1 (MQ=255)
                                                                                                                                                    gCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGGCAGGTAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGt  <  1:769433/150‑1 (MQ=255)
                                                                                                                                                    gCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGt  <  2:1150030/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGT  >  NC_002947/3989408‑3989705

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