Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 296,371 C→G 5.7% E19D (GAG→GAC) ‡ PP_0242 ← TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947296,3710CG 5.7% 285.2 / 7.0 90E19D (GAG→GAC) ‡PP_0242TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (36/48);  new base G (2/3);  total (39/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

GCAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCG  >  NC_002947/296222‑296520
                                                                                                                                                       |                                                                                                                                                     
gCAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGACGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGcc                                                                                                                                                       >  2:2323320/1‑150 (MQ=255)
 cAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCt                                                                                                                                                      <  2:584917/150‑1 (MQ=255)
 cAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCt                                                                                                                                                      <  2:168999/150‑1 (MQ=255)
          ccTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCAACCCTTGTTCGGCCAACAGCTCACGCGCAGAC‑‑GCCAGGATCGCCTCGATcctgg                                                                                                                                             >  1:2299054/1‑150 (MQ=255)
           cTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggc                                                                                                                                            >  1:2282087/1‑150 (MQ=255)
           cTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                            <  1:1342979/150‑1 (MQ=255)
            tCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAAGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCGAAGAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcc                                                                                                                                           >  1:633384/1‑150 (MQ=255)
               aGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcctgg                                                                                                                                        >  2:1472030/1‑150 (MQ=255)
               aGCAAGGCCGGCACGCTGACGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcctgg                                                                                                                                        <  1:2323320/150‑1 (MQ=255)
                    ggCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCCAGGCCAGGATCGAGGCGATcctggcggctgcg                                                                                                                                     <  2:99505/150‑8 (MQ=38)
                       cGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGgc                                                                                                                                <  2:184876/150‑1 (MQ=255)
                        ggCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGgcg                                                                                                                               <  2:51679/150‑1 (MQ=255)
                                    aaaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTgcgc                                                                                                                   <  1:100852/150‑1 (MQ=255)
                                    aaaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTgcgc                                                                                                                   >  1:2395343/1‑150 (MQ=255)
                                     aaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGcgcg                                                                                                                  <  1:952632/150‑1 (MQ=255)
                                     aaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGcgcg                                                                                                                  <  2:1222833/150‑1 (MQ=255)
                                      aaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                 >  2:766293/1‑150 (MQ=255)
                                       aaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCggg                                                                                                                <  1:766293/150‑1 (MQ=255)
                                       aaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCggg                                                                                                                <  1:2096317/150‑1 (MQ=255)
                                           gATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtgg                                                                                                            <  2:2282087/150‑1 (MQ=255)
                                            aTACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggt                                                                                                           <  2:2395343/150‑1 (MQ=255)
                                            aTACACCGACGAAGGCGGAATCTGCGCCCGCTCAGCCACGCTGTAGCTATAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggt                                                                                                           >  1:792276/1‑150 (MQ=255)
                                             tACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                          <  1:591501/150‑1 (MQ=255)
                                                acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATGCTGGGCTGACTGTTGTCTCGTGGTTTGCGCGGGgtggtggtg                                                                                                       >  2:998468/1‑150 (MQ=255)
                                                acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                 >  2:1563278/1‑140 (MQ=255)
                                                acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                 <  1:1563278/140‑1 (MQ=255)
                                                acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGgtggtggtg                                                                                                       >  2:1411737/1‑150 (MQ=255)
                                                  cgacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                                     <  2:1083112/150‑1 (MQ=255)
                                                   gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                                     >  2:516477/1‑149 (MQ=255)
                                                   gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                                     <  1:516477/149‑1 (MQ=255)
                                                   gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCGGTAGATAGAACGTGCTGCCACCCCTTATTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCgg                                                                                                    >  1:425939/1‑150 (MQ=255)
                                                         ggCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                              >  1:741419/1‑150 (MQ=255)
                                                         ggCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                              <  1:868170/150‑1 (MQ=255)
                                                          gCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatca                                                                                             >  2:1941509/1‑150 (MQ=255)
                                                            ggAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                           <  2:393655/150‑1 (MQ=255)
                                                              aagcTGCGCCCGCTCGGCCACGCTGTAGATAGAGCCTCCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTGTGCGCGGGGCGGTGGTGCGGTTCATCAActg                                                                                         <  2:633384/147‑1 (MQ=255)
                                                               aTCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAActgc                                                                                        <  1:1533419/150‑1 (MQ=255)
                                                                    cgcCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGCCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGCGGTGGTGGTGCGGTTCATCAACTGCTGAtt                                                                                   <  2:425939/150‑1 (MQ=255)
                                                                        cGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCa                                                                               <  1:956996/150‑1 (MQ=255)
                                                                         gCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGGTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAg                                                                              <  1:924456/150‑1 (MQ=255)
                                                                           tCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                              >  1:821716/1‑132 (MQ=255)
                                                                           tCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                              <  2:821716/132‑1 (MQ=255)
                                                                               ccACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCGATTTCTGATTGCCAGGATGCt                                                                        >  2:2361293/1‑150 (MQ=255)
                                                                               ccACGCTGGAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCATCTGCTGATTGCCAGGATGCt                                                                        <  2:754605/150‑1 (MQ=255)
                                                                                cACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTg                                                                       <  2:608153/150‑1 (MQ=255)
                                                                                  cGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGc                                                                     <  1:1472030/150‑1 (MQ=255)
                                                                                     tGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGtt                                                                                                  <  1:1211865/118‑1 (MQ=255)
                                                                                     tGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGtt                                                                                                  >  2:1211865/1‑118 (MQ=255)
                                                                                      gTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                            <  2:1125176/75‑1 (MQ=38)
                                                                                      gTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                            >  1:1125176/1‑75 (MQ=38)
                                                                                        agatagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTa                                                               <  2:1504465/150‑1 (MQ=255)
                                                                                          atagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcat                                                                                               <  1:736724/116‑1 (MQ=255)
                                                                                          atagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcat                                                                                               >  2:736724/1‑116 (MQ=255)
                                                                                           tagaAAGTGCTGCCACGCTTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGgc                                                            <  2:792276/150‑1 (MQ=255)
                                                                                             gaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAct                                                                                          >  1:904468/1‑118 (MQ=255)
                                                                                             gaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAct                                                                                          <  2:904468/118‑1 (MQ=255)
                                                                                              aaaGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGcc                                                         <  2:1778904/150‑1 (MQ=255)
                                                                                                 gTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCcgac                                                      <  1:2059781/150‑1 (MQ=255)
                                                                                                     tgcCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAAc                                                  >  2:1498485/1‑150 (MQ=255)
                                                                                                         aCCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGt                                              >  2:1562816/1‑150 (MQ=255)
                                                                                                           cccTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCg                                            >  2:1887518/1‑150 (MQ=255)
                                                                                                             cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                          >  2:1197414/1‑86 (MQ=255)
                                                                                                             cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                          <  1:1197414/86‑1 (MQ=255)
                                                                                                             cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGaa                                          <  2:957091/150‑1 (MQ=255)
                                                                                                                gTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                           <  1:943738/98‑1 (MQ=255)
                                                                                                                gTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                           >  2:943738/1‑98 (MQ=255)
                                                                                                                 ttCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCa                                      >  1:2106263/1‑150 (MQ=255)
                                                                                                                    ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGa                                           <  2:1580699/142‑1 (MQ=255)
                                                                                                                    ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGa                                           >  1:1580699/1‑142 (MQ=255)
                                                                                                                    ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCg                                   >  2:278930/1‑150 (MQ=255)
                                                                                                                          cAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCg                             <  2:741419/150‑1 (MQ=255)
                                                                                                                          cAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCg                             >  2:2085567/1‑150 (MQ=255)
                                                                                                                            gCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCGGTGACTTCgac                           <  2:471893/150‑3 (MQ=255)
                                                                                                                            gCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtt                           <  2:2106263/150‑1 (MQ=255)
                                                                                                                             cTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                          <  1:1411737/150‑1 (MQ=255)
                                                                                                                              tCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtttt                         >  1:1182031/1‑150 (MQ=255)
                                                                                                                              tCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGGTCATCAACTGCTGATTGCCATGAGGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtttt                         >  1:1741818/1‑150 (MQ=255)
                                                                                                                                aCGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                          <  2:2105572/147‑1 (MQ=255)
                                                                                                                                aCGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCCGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                          >  1:2105572/1‑147 (MQ=255)
                                                                                                                                  gcgacgcc‑‑gctAGGATCGCTTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTc                     >  2:600258/6‑150 (MQ=255)
                                                                                                                                      aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAg                                                <  2:1235782/119‑1 (MQ=255)
                                                                                                                                      aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAg                                                >  1:1235782/1‑119 (MQ=255)
                                                                                                                                      aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGc                 >  2:2015883/1‑150 (MQ=255)
                                                                                                                                        cc‑‑gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTc               <  1:998468/150‑1 (MQ=255)
                                                                                                                                            gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGaa                                                   >  2:1570321/1‑112 (MQ=255)
                                                                                                                                            gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGaa                                                   <  1:1570321/112‑1 (MQ=255)
                                                                                                                                            gccAGGATCGAGGCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCgg             >  1:183914/1‑150 (MQ=255)
                                                                                                                                             ccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGa            <  1:1941509/150‑1 (MQ=255)
                                                                                                                                                 gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                 <  2:1053330/45‑1 (MQ=255)
                                                                                                                                                 gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                 >  1:1053330/1‑45 (MQ=255)
                                                                                                                                                 gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc        >  1:1042964/1‑150 (MQ=255)
                                                                                                                                                       cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGTGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcg  >  2:767646/1‑150 (MQ=255)
                                                                                                                                                       cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGc                       >  1:109345/1‑129 (MQ=255)
                                                                                                                                                       cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGc                       <  2:109345/129‑1 (MQ=255)
                                                                                                                                                       cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcg  >  2:1764718/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCG  >  NC_002947/296222‑296520

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