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% 321.7 / 6.8 98P676A (CCC→GCC) oplA5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base C (51/42);  new base G (2/3);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCT  >  NC_002947/3989403‑3989700
                                                                                                                                                     |                                                                                                                                                    
ccACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAc                                                                                                                                                      >  1:236714/1‑150 (MQ=255)
 cACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAcc                                                                                                                                                     >  1:2360896/1‑150 (MQ=255)
     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                                                                                                                                     <  1:1630523/150‑1 (MQ=255)
                 cgCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGc                                                                                                                                     <  2:1684432/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:2196925/1‑150 (MQ=255)
                                 gTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCa                                                                                                                     <  1:2207821/150‑1 (MQ=255)
                                  tGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCAc                                                                                                                    >  2:420056/1‑150 (MQ=255)
                                      gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCATACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGATCGCCTGGAGCCTGCGGCCACCTAc                                                                                                                >  2:1261098/1‑150 (MQ=255)
                                      gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCg                                                                                                                         <  1:2260990/141‑1 (MQ=255)
                                      gACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCg                                                                                                                         >  2:2260990/1‑141 (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:2175797/150‑1 (MQ=255)
                                         ggTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCgg                                                                                                             <  2:1551531/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                                                                                                                                      >  1:145121/1‑108 (MQ=255)
                                                      ccACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCAGCCGGGTGGCAggacaa                                                                                                                                      <  2:145121/112‑5 (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)
                                                                          aaCCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCg                                                                                                                                                   <  2:868126/79‑1 (MQ=255)
                                                                          aaCCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCg                                                                                                                                                   >  1:868126/1‑79 (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)
                                                                                cTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCg                                                                                                                                                   >  2:481307/1‑73 (MQ=255)
                                                                                cTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCg                                                                                                                                                   <  1:481307/73‑1 (MQ=255)
                                                                                    tGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCa                                                                  <  1:844430/150‑1 (MQ=255)
                                                                                       cTGGGCAGACCGTGCATGGCCCCTGTGTGCTCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATgg                                                               <  2:117134/150‑1 (MQ=255)
                                                                                        tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                              >  1:1035777/1‑150 (MQ=255)
                                                                                        tGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATggg                                                              <  2:793375/150‑1 (MQ=255)
                                                                                              gACCGTGCATGGCCCCTGTGTGATCGAGGAAGAGACGACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGAt                                                        >  1:2019040/1‑150 (MQ=255)
                                                                                                ccGTGCATGGCCCCTGTGCGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATaa                                                      <  1:1719095/150‑1 (MQ=255)
                                                                                                   tGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                    >  2:2369327/1‑117 (MQ=255)
                                                                                                   tGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCg                                                                                    <  1:2369327/117‑1 (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                                                                                                  >  2:717906/1‑83 (MQ=255)
                                                                                                                       gAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCgg                                                                                                  <  1:717906/83‑1 (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                           >  2:533091/1‑150 (MQ=255)
                                                                                                                           aaGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGc                           <  1:1814246/150‑1 (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:1054952/1‑150 (MQ=255)
                                                                                                                                  caccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGt                    >  2:679557/1‑150 (MQ=255)
                                                                                                                                  caccaATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGt                    >  1:636232/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)
                                                                                                                                                     |                                                                                                                                                    
CCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACACCACCAATATTGTCCTGCCACCCGGCTGGCAGGCTCGCCTGGAGCCTTCGGCCACCTACCGGGTGACACGCGGGGGTTGATCTTGCGGTCGGCGCGCGCCGATCATGGGCCTGATAACCGCCCGGCAGTTTTGCGCTGCCGGGCGGTTCGTACAATAACAAGAGCTGCT  >  NC_002947/3989403‑3989700

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