Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,989,552 C→G 5.1% P676A (CCC→GCC)  oplA → 5‑oxoprolinase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,989,5520CG 5.1% 260.4 / 3.6 79P676A (CCC→GCC) oplA5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base G (2/2);  total (37/42)
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.79e-01

CCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGC  >  NC_002947/3989403‑3989699
                                                                                                                                                     |                                                                                                                                                   
ccACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAc                                                                                                                                                     >  2:1482929/1‑150 (MQ=255)
 cACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAcc                                                                                                                                                    >  1:2265993/1‑150 (MQ=255)
   ccGCGCACTCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACGCCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCg                                                                                                                                                  <  2:1810884/150‑1 (MQ=255)
         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                                                                                                               <  1:970221/150‑1 (MQ=255)
                                      gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTAc                                                                                                               <  2:1133231/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                                                                  >  1:2313445/1‑150 (MQ=255)
                                                                                   ttGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATc                                                                  >  2:2120169/1‑150 (MQ=255)
                                                                                      ccTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                               <  1:239723/150‑1 (MQ=255)
                                                                                      ccTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                               <  1:1867339/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                                                               <  1:26592/149‑1 (MQ=255)
                                                                                       cTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                               >  2:26592/1‑149 (MQ=255)
                                                                                        tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                             >  1:1846365/1‑150 (MQ=255)
                                                                                        tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                             >  2:1441826/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                                              >  1:1164384/1‑146 (MQ=255)
                                                                                                           cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCg                                              <  2:1164384/146‑1 (MQ=255)
                                                                                                           cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAg                                          >  2:1589669/1‑150 (MQ=255)
                                                                                                           cccTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAg                                          >  2:603920/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                                                               <  1:2006889/121‑1 (MQ=255)
                                                                                                                   gATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATg                                                               >  2:2006889/1‑121 (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                                                                                                          >  2:254935/1‑62 (MQ=255)
                                                                                                                                   accaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGt                                                                                                          <  1:254935/62‑1 (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)
                                                                                                                                                   cACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAagag      >  1:991633/1‑146 (MQ=255)
                                                                                                                                                   cACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAagag      <  2:991633/146‑1 (MQ=255)
                                                                                                                                                   cACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAgctgc  <  1:133574/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGC  >  NC_002947/3989403‑3989699

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: