Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,297 A→G 5.3% W13R (TGG→CGG)  PP_3717 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,2970AG 5.3% 228.6 / 2.1 75W13R (TGG→CGG) PP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (40/31);  new base G (2/2);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACGGTGGCG  >  NC_002947/4240148‑4240445
                                                                                                                                                        |                                                                                                                                                    
ggCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCTGCAGTCTAGCCTGCCa                                                                                                                                                      <  1:1246021/150‑1 (MQ=255)
 gCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGAAGCGCGCCCATTGCGAGGCTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAg                                                                                                                                                     <  1:237581/150‑1 (MQ=255)
   aGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGc                                                                                                                                                   <  1:1780815/150‑1 (MQ=255)
     cGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGAGGTAAGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAgatcg                                                                                                                                                 >  1:990964/1‑146 (MQ=255)
         cacaTACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGccc                                                                                                                                             <  1:1178958/150‑1 (MQ=255)
          acaTACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCcg                                                                                                                                            >  2:888152/1‑150 (MQ=255)
             tcagccAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCgcgc                                                                                                                                         >  1:1542936/5‑150 (MQ=255)
                tCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGcaca                                                                                                                                      <  1:1580471/150‑1 (MQ=255)
                   aGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGcc                                                                                                                                   <  2:1542936/150‑1 (MQ=255)
                        cTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTg                                                                                                                              >  1:1329815/1‑150 (MQ=255)
                        cTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTg                                                                                                                              >  1:922548/1‑150 (MQ=255)
                         tCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTgc                                                                                                                             >  1:1516383/1‑150 (MQ=255)
                             gggTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCa                                                                                                                         >  2:1396099/1‑150 (MQ=255)
                             cccgcttaccctcacaatGGGCGCGCTGCAGCAGTGGCACGTGGC‑‑‑GTTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCa                                                                                                                         <  1:258868‑M1/132‑1 (MQ=255)
                                    aCCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCAt                                                                                                                  >  2:1054698/1‑150 (MQ=255)
                                           tcacaatgggcgcgctgcagcagtggcacgtggcgttgcagcgcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGc                                                                                                              <  2:258868‑M2/102‑1 (MQ=255)
                                                   GCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGCCGggcatga                                                                                                                                               >  2:1941586‑M1/1‑103 (MQ=255)
                                                   GCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGCCGggcatga                                                                                                                                               <  1:1941586‑M1/106‑4 (MQ=255)
                                                         CAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAgc                                                                                             >  2:1516383‑M1/1‑150 (MQ=255)
                                                            cagcagtggcacgtggcgttgcagcgcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAgc                                                                                             >  1:1625383‑M2/32‑150 (MQ=255)
                                                            c‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGcc                                                                                          <  2:1329815/150‑1 (MQ=255)
                                                                 ccACGTGCCACTGTTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGt                                                                                        >  1:15095/1‑150 (MQ=255)
                                                                 ccACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCCCACAAACTTTTGCGCCAGGTTCATCGGCCTTGGCACAGGCCAAGCGCCGt                                                                                        <  1:1901249/150‑1 (MQ=255)
                                                                   aCGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTg                                                                                                          <  1:1561968/130‑1 (MQ=255)
                                                                   aCGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTg                                                                                                          >  2:1561968/1‑130 (MQ=255)
                                                                       gCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTc                                                                                  >  2:761727/1‑150 (MQ=255)
                                                                          aCTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcg                                                                               >  1:672530/1‑150 (MQ=255)
                                                                          aCTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcg                                                                               >  1:1115718/1‑150 (MQ=255)
                                                                              ctgcAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcgcgcg                                                                           <  1:761727/150‑1 (MQ=255)
                                                                                gcagcgcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGt                                                                         >  2:476643‑M2/12‑150 (MQ=255)
                                                                                     gcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                                    <  1:1054698‑M2/144‑1 (MQ=255)
                                                                                     gcgcccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                                    <  1:1396099‑M2/144‑1 (MQ=255)
                                                                                         ccATTGTGAGGGTAAGCGGGGTATCGGGAGGGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                                                >  1:1724341‑M2/3‑150 (MQ=255)
                                                                                         ccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAat                                                                   >  1:1704604‑M2/3‑147 (MQ=255)
                                                                                         ccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAat                                                                   <  2:1704604‑M2/145‑1 (MQ=255)
                                                                                          cATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGTGCCAGGTTCATCGGCCTTg                                                                                                          <  1:1034515‑M2/106‑1 (MQ=255)
                                                                                          cATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGTGCCAGGTTCATCGGCCTTg                                                                                                          >  2:1034515‑M2/2‑107 (MQ=255)
                                                                                           aTTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGc                                                              >  1:1211834/1‑150 (MQ=255)
                                                                                            ttGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCa                                                             >  1:1145190/1‑150 (MQ=255)
                                                                                             tgtgAGGGTAAGCGGGGTATCGGGAGGGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCAc                                                            <  2:1724341/150‑1 (MQ=255)
                                                                                             tgtgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCAc                                                            >  1:135203/1‑150 (MQ=255)
                                                                                              gtgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                           >  2:757227/1‑150 (MQ=255)
                                                                                               tgAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                                          <  2:1625383/150‑1 (MQ=255)
                                                                                                    gTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCa                                                     >  1:349212/1‑150 (MQ=255)
                                                                                                    gTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCa                                                     >  2:656102/1‑150 (MQ=255)
                                                                                                     tAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAg                                                    <  2:187644/150‑1 (MQ=255)
                                                                                                           gggTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGtctc                                              <  1:888152/150‑1 (MQ=255)
                                                                                                                  gggAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCTACGAtc                                       <  1:476643/150‑1 (MQ=255)
                                                                                                                    gagaGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtctc                                     >  1:327096/1‑150 (MQ=255)
                                                                                                                       agGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                                          <  2:868836/126‑1 (MQ=255)
                                                                                                                       agGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTg                                                          >  1:868836/1‑126 (MQ=255)
                                                                                                                         gCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtc                                       >  1:868945/1‑143 (MQ=255)
                                                                                                                         gCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtc                                       <  2:868945/143‑1 (MQ=255)
                                                                                                                             gggTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTg                            >  1:742301/1‑150 (MQ=255)
                                                                                                                                tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                                    <  2:796906/107‑1 (MQ=255)
                                                                                                                                tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                                    >  1:796906/1‑107 (MQ=255)
                                                                                                                                    gCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATgg                     >  2:184384/1‑150 (MQ=255)
                                                                                                                                    gCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATgg                     <  1:561940/150‑1 (MQ=255)
                                                                                                                                     cccGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATggg                    >  2:1010639/1‑150 (MQ=255)
                                                                                                                                      ctggCAGGCTAGACTGCCGGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcgcgcg                                                                           <  2:1918634/92‑1 (MQ=255)
                                                                                                                                      ctggCAGGCTAGACTGCCGGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCAcgcgcg                                                                           >  1:1918634/3‑94 (MQ=255)
                                                                                                                                      ccGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGt                   >  2:938945/1‑150 (MQ=255)
                                                                                                                                        ggCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatga                                                                 >  2:1897874/1‑102 (MQ=255)
                                                                                                                                        ggCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatga                                                                 <  1:1897874/102‑1 (MQ=255)
                                                                                                                                         gCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGtgat                <  2:135203/150‑1 (MQ=255)
                                                                                                                                             tCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATc            >  2:325545/1‑150 (MQ=255)
                                                                                                                                             tCTAGCCTGCCAGGCACGCCCGCGCACAGCACTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                           >  1:126134/1‑103 (MQ=255)
                                                                                                                                             tCTAGCCTGCCAGGCACGCCCGCGCACAGCACTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                           <  2:126134/103‑1 (MQ=255)
                                                                                                                                              cTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                                                >  2:691949/1‑97 (MQ=255)
                                                                                                                                              cTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAatgat                                                                <  1:691949/97‑1 (MQ=255)
                                                                                                                                                aGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAAc         <  2:349212/150‑1 (MQ=255)
                                                                                                                                                 gCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGtc                                                >  1:966586/1‑110 (MQ=255)
                                                                                                                                                 gCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGtc                                                <  2:966586/110‑1 (MQ=255)
                                                                                                                                                 gCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACg        <  1:184384/150‑1 (MQ=255)
                                                                                                                                                    tGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTACATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACGGTg     >  1:1570666/1‑150 (MQ=255)
                                                                                                                                                       cAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACGGTGGCg  >  1:1162040/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATCAACGGTGGCG  >  NC_002947/4240148‑4240445

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