Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 439,419 A→T 7.9% N185K (AAT→AAA PP_0361 ← competence protein ComF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947439,4190AT 7.9% 261.0 / 10.0 88N185K (AAT→AAAPP_0361competence protein ComF
Reads supporting (aligned to +/- strand):  ref base A (42/39);  new base T (3/4);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCCAGGCAATACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTCGCGCAGTC  >  NC_002947/439270‑439567
                                                                                                                                                     |                                                                                                                                                    
gCCAGGCAATACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTAGCGCaaa                                                                                                                                                      <  2:1424537/150‑1 (MQ=255)
  cAGGCAATACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAAtt                                                                                                                                                    <  1:677158/150‑1 (MQ=255)
  cAGGCAATACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCAACGGTGCTACAGGCAGCAGCGACCATGGCGACTTACTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGACGCAAAtt                                                                                                                                                    >  2:152764/1‑150 (MQ=255)
    ggCAATACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTgc                                                                                                                                                  <  1:460020/150‑1 (MQ=255)
         tACCCATCCCCGCGCCCCGCACCTGCCTTGCGCTGTGCCCGTGCGATAGCCTGCGAATTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGCACGCCTGGCGCAATTTGcgccag                                                                                                                                             <  1:1643779/150‑3 (MQ=255)
         tACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAATTTGcgccag                                                                                                                                             <  2:710677/150‑3 (MQ=255)
         tACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCCGGCGCAATTTGcgccag                                                                                                                                             <  1:404632/150‑3 (MQ=255)
          acacATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCgccgcc                                                                                                                                            >  2:1361928/1‑150 (MQ=255)
            acaTCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGt                                                                                                                                          >  1:826420/1‑150 (MQ=255)
               tCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTgc                                                                                                                                                  <  1:1372016/139‑1 (MQ=255)
               tCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTgc                                                                                                                                                  <  1:1372066/139‑1 (MQ=255)
               tCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTgc                                                                                                                                                  >  2:1372016/1‑139 (MQ=255)
               tCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTgc                                                                                                                                                  >  2:1372066/1‑139 (MQ=255)
                  aCCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTg                                                                                                                                    <  1:879457/150‑1 (MQ=255)
                  aCCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTg                                                                                                                                    >  1:582969/1‑150 (MQ=255)
                  aCCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTg                                                                                                                                    >  2:155407/1‑150 (MQ=255)
                      gccgcGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGgcgc                                                                                                                                >  2:185688/1‑150 (MQ=255)
                         gctcACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCa                                                                                                                             <  1:185688/147‑1 (MQ=255)
                         gcgcACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCa                                                                                                                             <  2:1110104/150‑1 (MQ=255)
                            cACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGAt                                                                                                                          <  1:606378/150‑1 (MQ=255)
                            cACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGAt                                                                                                                          <  1:60502/150‑1 (MQ=255)
                            cACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGAt                                                                                                                          <  2:826420/150‑1 (MQ=255)
                                tGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTg                                                                                                                      <  1:1430392/150‑1 (MQ=255)
                                tGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTCCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTg                                                                                                                      <  1:1870080/150‑1 (MQ=255)
                                 gCCTTGCGCAGTGCCTGGGCGATAGCCTGGGCAGTGGCACCGGTGGGCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGc                                                                                                                     <  2:924008/150‑1 (MQ=255)
                                  ccTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCt                                                                                                                    >  1:1214987/1‑150 (MQ=255)
                                         cAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACACCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCgg                                                                                                             <  1:8619/150‑1 (MQ=255)
                                              ccTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGtct                                                                                                        >  1:367848/1‑150 (MQ=255)
                                               cTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGtctc                                                                                                       <  1:231267/150‑1 (MQ=255)
                                                tGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTcg                                                                                                      >  1:1111696/1‑150 (MQ=255)
                                                         gCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCgg                                                                                             >  2:1077914/1‑150 (MQ=255)
                                                           cTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGc                                                                                           >  1:790917/1‑150 (MQ=255)
                                                            tGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGCACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGcc                                                                                          >  1:1067915/1‑150 (MQ=255)
                                                             gcgcAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCa                                                                                         <  1:188298/150‑1 (MQ=255)
                                                                     gCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGc                                                                                 >  2:673140/1‑150 (MQ=255)
                                                                       aCCGGTGGTCAACACATCTGCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGCTGTTGCTGTGCGGGTGTCTCGCGGGTGCGGGTGTGCAAGCGCTc                                                                               <  2:594373/150‑1 (MQ=255)
                                                                       aCCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTTGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCa                                                                                                                             <  2:1071337/104‑1 (MQ=255)
                                                                       aCCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCATTTTCCCTCTCCTTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTc                                                                               >  1:1278637/1‑150 (MQ=255)
                                                                       aCCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCa                                                                                                                             >  1:1071337/1‑104 (MQ=255)
                                                                       aCCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTc                                                                               >  1:1390463/1‑150 (MQ=255)
                                                                        ccGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCa                                                                                         <  2:1072122/139‑1 (MQ=255)
                                                                        ccGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCa                                                                                         >  1:1072122/1‑139 (MQ=255)
                                                                             ggtCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCa                                                                         <  1:1361928/150‑1 (MQ=255)
                                                                               tCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAgg                                                                       >  1:668004/1‑150 (MQ=255)
                                                                                 aaCACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGtt                                                                     >  1:917236/1‑150 (MQ=255)
                                                                                     catcatCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTAccc                                                                 <  2:1111696/150‑1 (MQ=255)
                                                                                      atcatcGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTAcccc                                                                <  1:155407/150‑1 (MQ=255)
                                                                                         atcGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGt                                                                                                                        >  2:626726/1‑91 (MQ=255)
                                                                                         atcGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGt                                                                                                                        <  1:626726/91‑1 (MQ=255)
                                                                                         atcGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGAc                                                             >  2:1472609/1‑150 (MQ=255)
                                                                                          tcGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACg                                                            >  2:1754176/1‑150 (MQ=255)
                                                                                             acgaTGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGc                                                         <  1:1237660/150‑1 (MQ=255)
                                                                                                aTGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGcgac                                                      <  1:1429252/150‑1 (MQ=255)
                                                                                                   gCGACATGCTTGCCTTCCAGTTCGCCCTGTCCCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCgacgac                                                   <  1:1530237/150‑1 (MQ=255)
                                                                                                      aCATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAac                                                <  1:503721/150‑1 (MQ=255)
                                                                                                      aCATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAac                                                <  2:1070200/150‑1 (MQ=255)
                                                                                                            ttGCCTTCCAGTTCGCCCTGTCCCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGc                                          <  1:304342/150‑1 (MQ=255)
                                                                                                             tGCCTTCCAGTTCGCCCTGTACCTCGAACGCCTGGCGCAAATTGCACCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACGACCGcc                                         >  1:290838/1‑150 (MQ=255)
                                                                                                             tGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGcc                                         >  1:851428/1‑150 (MQ=255)
                                                                                                             tGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGcc                                         >  2:761401/1‑150 (MQ=255)
                                                                                                                cTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCaa                                      >  1:1568381/1‑150 (MQ=255)
                                                                                                                cTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCaa                                      <  2:1722278/150‑1 (MQ=255)
                                                                                                                 ttCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTc                                                                               >  2:297943/1‑108 (MQ=255)
                                                                                                                 ttCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTc                                                                               <  1:297943/108‑1 (MQ=255)
                                                                                                                 ttCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGcga                                                       >  1:1737899/1‑132 (MQ=255)
                                                                                                                 ttCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGcga                                                       <  2:1737899/132‑1 (MQ=255)
                                                                                                                  tCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGc                                    >  2:1549020/1‑150 (MQ=255)
                                                                                                                   ccAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGACCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCa                                   >  1:761524/1‑150 (MQ=255)
                                                                                                                     aGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATc                                 <  1:1754176/150‑1 (MQ=255)
                                                                                                                     aGTTCGCCCGGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATc                                 <  2:898681/150‑1 (MQ=255)
                                                                                                                           cccTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGcc                           >  2:1610175/1‑150 (MQ=255)
                                                                                                                            ccTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCt                          >  2:1750674/1‑150 (MQ=255)
                                                                                                                            ccTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCt                          >  1:772619/1‑150 (MQ=255)
                                                                                                                                tACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGtt                      <  2:668004/150‑1 (MQ=255)
                                                                                                                                tACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGtt                      <  1:301396/150‑1 (MQ=255)
                                                                                                                                  ccGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGa                    >  1:1340483/1‑150 (MQ=255)
                                                                                                                                  ccGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGa                    >  2:284644/1‑150 (MQ=255)
                                                                                                                                    gcgAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGTCGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAAc                  >  1:1563249/1‑150 (MQ=255)
                                                                                                                                    gcgAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAAc                  >  1:411469/1‑150 (MQ=255)
                                                                                                                                    gcgAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAAc                  >  1:1010573/1‑150 (MQ=255)
                                                                                                                                       aaCGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCGAGCATCCCCGCCTGGTTGAAcccc               <  2:1288349/150‑1 (MQ=255)
                                                                                                                                         acggcggCGCAATTTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAAccccc              >  1:636152/6‑149 (MQ=255)
                                                                                                                                         acggcggCGCAATTTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAAccccc              <  2:636152/144‑1 (MQ=255)
                                                                                                                                          gCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGc            <  1:1610175/150‑1 (MQ=255)
                                                                                                                                            cTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTc          <  1:1136315/150‑1 (MQ=255)
                                                                                                                                             tGGCGCAAATTGCGCCGCCGTGCCTTAGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGACCAGGCATCCCCTCCTGGGTGAACCCCCGCTcc         >  2:957966/1‑149 (MQ=255)
                                                                                                                                              ggCGCAATTTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACCACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTcgc        >  2:1568376/1‑150 (MQ=255)
                                                                                                                                               gcgcAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGAGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTcgcg       >  2:146015/1‑150 (MQ=255)
                                                                                                                                                 gcAAATTGCGCCGCCGTGCCTTGACGACAAAGGGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTCGCGCa     >  1:1790292/1‑150 (MQ=255)
                                                                                                                                                   aaaTTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCTCCGCCCATGCATCCCCTTCTGGTCCAGCCCCCGCTCGCGCAGt   <  1:152764/150‑1 (MQ=255)
                                                                                                                                                   aaaTTGCGCCGCCGTGCCTTGGCGAACAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTCGCGCAGt   <  1:293250/150‑1 (MQ=255)
                                                                                                                                                    aaTTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTGGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTCGCGCAGTc  >  1:554024/1‑150 (MQ=255)
                                                                                                                                                    aaTTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTCGCGCAGTc  <  1:872025/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCAGGCAATACACATCCACCCGCCGCGCACCTGCCTTGCGCAGTGCCTGGGCGATAGCCTGCGCAGTGGCACCGGTGGTCAACACATCATCGACGATGGCGACATGCTTGCCTTCCAGTTCGCCCTGTACCGCGAACGCCTGGCGCAAATTGCGCCGCCGTGCCTTGGCGCCCAGATGTTGCTGTGCCGGTGTCTCGCGGGTGCGGGCCAGCAAACGCTCGTTGCAGGTTACCCCGACGTGCGACGACAACCACCGCCCAAGCATCCCCGCCTGGTTGAACCCCCGCTCGCGCAGTC  >  NC_002947/439270‑439567

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: