Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,286 T→G 7.2% R16S (AGA→AGC 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,2860TG 7.2% 197.6 / 6.5 69R16S (AGA→AGCPP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base G (2/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

AACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATC  >  NC_002947/4240137‑4240435
                                                                                                                                                        |                                                                                                                                                     
aaCTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGt                                                                                                                                                       <  1:203488/150‑1 (MQ=255)
    ggCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAg                                                                                                                                                   <  2:1080292/150‑1 (MQ=255)
     gCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGc                                                                                                                                                  >  2:1780815/1‑150 (MQ=255)
       cGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCt                                                                                                                                                >  2:561940/1‑150 (MQ=255)
           ggCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCTGCAGTCTAGCCTGCCa                                                                                                                                            <  1:1246021/150‑1 (MQ=255)
            gCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCTGGCAGGCTAGACTGccgg                                                                                                                                           <  2:276536/150‑3 (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                                                         <  2:1704604‑M2/145‑1 (MQ=255)
                                                                                                    ccATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAat                                                         >  1:1704604‑M2/3‑147 (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                                           >  2:656102/1‑150 (MQ=255)
                                                                                                               gTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCa                                           >  1:349212/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                             <  2:868945/143‑1 (MQ=255)
                                                                                                                                    gCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGAtc                             >  1:868945/1‑143 (MQ=255)
                                                                                                                                        gggTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTg                  >  1:742301/1‑150 (MQ=255)
                                                                                                                                           tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                          >  1:796906/1‑107 (MQ=255)
                                                                                                                                           tCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCaa                                                          <  2:796906/107‑1 (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                                                 <  2:126134/103‑1 (MQ=255)
                                                                                                                                                        tCTAGCCTGCCAGGCACGCCCGCGCACAGCACTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACt                                                 >  1:126134/1‑103 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
AACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAC‑‑‑GCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGGCCTTGGTACAGGCCAAGCGCCGTGAGTTCACGCGCGGTCTCAATGATGCACTGGCGCAGGGTCTCGACGATCTCGTCCACCTGCGCATGGGTGATGATC  >  NC_002947/4240137‑4240435

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: