Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,989,557 C→G 6.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 6.0% 270.0 / 4.8 83G677G (GGC→GGG) oplA5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base C (37/41);  new base G (2/3);  total (39/44)
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.94e-01

AATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGT  >  NC_002947/3989412‑3989705
                                                                                                                                                 |                                                                                                                                                    
aaTCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGc                                                                                                                                                  <  1:1482929/150‑1 (MQ=255)
   ccgtccgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAgg                                                                                                                                               <  2:246377/150‑1 (MQ=255)
    cgtccgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAgga                                                                                                                                              <  1:2040807/150‑2 (MQ=255)
           ccGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTg                                                                                                                                       <  2:1246882/150‑1 (MQ=255)
               tGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGc                                                                                                                                   >  2:36879/1‑150 (MQ=255)
                  tGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCtt                                                                                                                                >  1:870059/1‑150 (MQ=255)
                    ttCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCg                                                                                                                              <  1:591844/150‑1 (MQ=255)
                       aGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGcc                                                                                                                           >  1:2176373/1‑150 (MQ=255)
                             gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTAc                                                                                                                     <  2:1133231/150‑1 (MQ=255)
                             gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTAc                                                                                                                     <  1:970221/150‑1 (MQ=255)
                                ggTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCgg                                                                                                                  >  1:2072419/1‑150 (MQ=255)
                                  tGAATGGAAGCAAACACAGGAGTTCAACGGCAACCACCTGTTGGCTGAGCAGACGGTGCATGGCCCCGGGGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGt                                                                                                                >  2:25512/1‑150 (MQ=255)
                                   gAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTg                                                                                                               >  1:1512772/1‑150 (MQ=255)
                                      tGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGaca                                                                                                            >  2:1443787/1‑150 (MQ=255)
                                             ccACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAAAATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCgggg                                                                                                     <  2:1384614/150‑1 (MQ=255)
                                              cacaCCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGa                                                                                                              >  1:1455767/1‑140 (MQ=255)
                                              cacaCCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGa                                                                                                              <  2:1455767/140‑1 (MQ=255)
                                              cacaCCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCggggg                                                                                                    <  1:829875/150‑1 (MQ=255)
                                               acacCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGt                                                                                                   >  2:1867339/1‑150 (MQ=255)
                                                cacCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGtt                                                                                                  <  2:2176373/150‑1 (MQ=255)
                                                 acCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTg                                                                                                 >  1:2123801/1‑150 (MQ=255)
                                                     gtgtTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGa                                                                                                              >  1:546894/1‑133 (MQ=255)
                                                     gtgtTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGa                                                                                                              <  2:546894/133‑1 (MQ=255)
                                                       gtTCAACGGCAACCACCTGCTGCCTGGGCAGACCGTGCATGGCTCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTg                                                                                           <  2:1587529/150‑1 (MQ=255)
                                                        ttCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGc                                                                                          <  2:870059/150‑1 (MQ=255)
                                                         tCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                         <  2:2194649/150‑1 (MQ=255)
                                                            aCGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCACGCTCCCCTTGAGCCTTAGGCCACCTACTGGGTGACACGCGGGGGTTGATCTTGCGGTc                                                                                      >  2:572179/1‑150 (MQ=255)
                                                              ggCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCgg                                                                                    >  2:2241184/1‑150 (MQ=255)
                                                                    caccTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGgcgcgcg                                                                              >  2:133574/1‑150 (MQ=255)
                                                                          ttGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATc                                                                        >  2:2120169/1‑150 (MQ=255)
                                                                          ttGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATc                                                                        >  1:2313445/1‑150 (MQ=255)
                                                                             ccTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     <  1:1867339/150‑1 (MQ=255)
                                                                             ccTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     <  1:239723/150‑1 (MQ=255)
                                                                             ccTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     <  1:1443787/150‑1 (MQ=255)
                                                                             ccTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGAGCATg                                                                     >  2:239723/1‑150 (MQ=255)
                                                                              cTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAgga                                                                                                                                              <  2:1930744/76‑2 (MQ=39)
                                                                              cTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAgga                                                                                                                                              >  1:1930744/1‑75 (MQ=255)
                                                                              cTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     >  2:26592/1‑149 (MQ=255)
                                                                              cTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     <  1:26592/149‑1 (MQ=255)
                                                                               tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                                   >  2:1441826/1‑150 (MQ=255)
                                                                               tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                                   >  1:1846365/1‑150 (MQ=255)
                                                                                  gCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCt                                                                <  1:1441826/150‑1 (MQ=255)
                                                                                      aCCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATa                                                            >  1:267419/1‑150 (MQ=255)
                                                                                       ccGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGt                                                                                                   >  1:916264/1‑110 (MQ=255)
                                                                                       ccGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGt                                                                                                   <  2:916264/110‑1 (MQ=255)
                                                                                       ccGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATaa                                                           <  2:2265993/150‑1 (MQ=255)
                                                                                       ccGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATaa                                                           <  1:2241184/150‑1 (MQ=255)
                                                                                        cGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAAc                                                          <  2:1846365/150‑1 (MQ=255)
                                                                                                  cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCg                                                    <  2:1164384/146‑1 (MQ=255)
                                                                                                  cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCg                                                    >  1:1164384/1‑146 (MQ=255)
                                                                                                  cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAg                                                >  2:603920/1‑150 (MQ=255)
                                                                                                  cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAg                                                >  2:1589669/1‑150 (MQ=255)
                                                                                                      gtgtgATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGtttt                                            <  2:2123801/150‑1 (MQ=255)
                                                                                                        gtgATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTgc                                          <  1:36879/150‑1 (MQ=255)
                                                                                                        gtgATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTgc                                          <  1:1589669/150‑1 (MQ=255)
                                                                                                          gATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     >  2:2006889/1‑121 (MQ=255)
                                                                                                          gATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     <  1:2006889/121‑1 (MQ=255)
                                                                                                           aTCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCt                                       <  1:2120169/150‑1 (MQ=255)
                                                                                                            tCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTg                                      >  1:1428023/1‑150 (MQ=255)
                                                                                                                   aGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCg                               <  1:603920/150‑1 (MQ=255)
                                                                                                                    gACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCgg                              >  1:2110797/1‑150 (MQ=255)
                                                                                                                          accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGt                                                                                                                <  1:254935/62‑1 (MQ=255)
                                                                                                                          accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGt                                                                                                                >  2:254935/1‑62 (MQ=255)
                                                                                                                          accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGcc                                                                 >  2:996494/1‑109 (MQ=255)
                                                                                                                          accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGcc                                                                 <  1:996494/109‑1 (MQ=255)
                                                                                                                          accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTa                        <  2:2313445/150‑1 (MQ=255)
                                                                                                                          accaATATTGTCCTGCCACCCGGCTGGCAGGATCCCCTGGAGAACCAGGCCACCTACCGGGTGACAGGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTa                        >  2:883168/1‑150 (MQ=255)
                                                                                                                           ccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     <  2:2023770/104‑1 (MQ=255)
                                                                                                                           ccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                                     >  1:2023770/1‑104 (MQ=255)
                                                                                                                               tatTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATa                   >  1:400154/1‑150 (MQ=255)
                                                                                                                                    tCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGAt                                                                         <  2:1752690/91‑1 (MQ=255)
                                                                                                                                    tCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGAt                                                                         >  1:1752690/1‑91 (MQ=255)
                                                                                                                                    tCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAag              <  2:2110797/150‑1 (MQ=255)
                                                                                                                                      cTGCCAGCCGGGTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAagag            >  2:1187448/1‑150 (MQ=255)
                                                                                                                                        gCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCACGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAgct          <  2:543788/150‑1 (MQ=255)
                                                                                                                                          cagccGGGTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAgctgc        <  1:1187448/147‑1 (MQ=255)
                                                                                                                                          cACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAagag            >  1:991633/1‑146 (MQ=255)
                                                                                                                                          cACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAagag            <  2:991633/146‑1 (MQ=255)
                                                                                                                                          cACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAgctgc        <  1:133574/150‑1 (MQ=255)
                                                                                                                                             ccGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAAt                    <  2:478939/135‑1 (MQ=255)
                                                                                                                                             ccGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAAt                    >  1:478939/1‑135 (MQ=255)
                                                                                                                                             ccGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTg     >  1:409099/1‑150 (MQ=255)
                                                                                                                                                gCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGt  <  2:1512772/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCTTGAGT  >  NC_002947/3989412‑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: ‑