Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,407,008 G→C 6.6% intergenic (+220/+34) PP_2110 → / ← PP_2111 ribosomal large subunit pseudouridine synthase/aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,407,0080GC 6.6% 236.2 / 6.2 76intergenic (+220/+34)PP_2110/PP_2111ribosomal large subunit pseudouridine synthase/aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (31/40);  new base C (2/3);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCGGCCAGTTCGG  >  NC_002947/2406861‑2407154
                                                                                                                                                   |                                                                                                                                                  
ccGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGCGCGGGTTCACCCGCGAACACCGGCagag                                                                                                                                                  >  2:625746/1‑150 (MQ=255)
         tATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACTGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCa                                                                                                                                         >  1:1032030/1‑150 (MQ=255)
             cgcgGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTNGTAGGAGCGGGTTCACCCGCGAACACCGGCTCTGCCGGTgttcgcgg                                                                                                                                     <  1:564975/150‑8 (MQ=38)
                     aaCCCGCTCCTCCATGGACCGGGGTGGAGGCGCATCAGTGGTTCGACACAGGCATCGACGGCCGGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGaa                                                                                                                             <  1:695526/150‑1 (MQ=255)
                     aaCCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGaa                                                                                                                             <  2:799958/150‑1 (MQ=255)
                        ccGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCa                                                                                                                                   <  1:1363076/141‑1 (MQ=255)
                        ccGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCa                                                                                                                                   >  2:1363076/1‑141 (MQ=255)
                        ccGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGtc                                                                                                                          >  2:1140163/1‑150 (MQ=255)
                          gCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGtctc                                                                                                                        <  2:751468/150‑1 (MQ=255)
                            tCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACTCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGa                                                                                                                      <  2:1174029/150‑1 (MQ=255)
                              cTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGc                                                                                                                    >  1:819304/1‑150 (MQ=255)
                                aCATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCa                                                                                                                                         <  1:552637/127‑1 (MQ=255)
                                aCATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCa                                                                                                                                         >  2:552637/1‑127 (MQ=255)
                                        tgggCTGGTCGCGCATCAGTTGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGtt                                                                                                          <  2:304506/149‑1 (MQ=255)
                                        cGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGtt                                                                                                          >  1:217596/1‑150 (MQ=255)
                                           gCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGaa                                                                                                       <  2:466316/150‑1 (MQ=255)
                                            cTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAg                                                                                                      <  1:1208742/150‑1 (MQ=255)
                                            cTGGACGCCAAGCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCCTCAGCGGTTGAAg                                                                                                      <  1:1130502/150‑1 (MQ=255)
                                                 cgcgCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTc                                                                                                 <  1:1233005/150‑1 (MQ=255)
                                                 cgcgCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGCTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTc                                                                                                 <  2:118678/150‑1 (MQ=255)
                                                  gcgcATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTcc                                                                                                <  2:385095/150‑1 (MQ=255)
                                                    gcATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCg                                                                                                                                <  2:488456/116‑1 (MQ=255)
                                                    gcATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCg                                                                                                                                >  1:488456/1‑116 (MQ=255)
                                                      aTCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCCGCGGTTGAAGCGCTccacca                                                                                            >  1:918834/1‑150 (MQ=255)
                                                       tCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGATCTTCAGCGGTTGAAGCGCTCCACCAg                                                                                           >  1:1265432/1‑150 (MQ=255)
                                                       tCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAg                                                                                           >  2:819992/1‑150 (MQ=255)
                                                        cAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGc                                                                                          <  2:34863/150‑1 (MQ=255)
                                                        cAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTATCTACCTTCAGCGGTTGAAGCGCTCCACCAGc                                                                                          <  1:260882/150‑1 (MQ=255)
                                                         aGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTg                                                                                                                                            >  2:1051688/1‑99 (MQ=255)
                                                         aGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTg                                                                                                                                            <  1:1051688/99‑1 (MQ=255)
                                                             gTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATa                                                                                     <  1:1240202/150‑1 (MQ=255)
                                                              ttCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATAc                                                                                    <  2:217596/150‑1 (MQ=255)
                                                               tCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACt                                                                                   <  1:1285108/150‑1 (MQ=255)
                                                                cGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGACGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTg                                                                                  <  2:1032030/150‑1 (MQ=255)
                                                                cGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCAGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCTGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTg                                                                                  <  2:916084/150‑1 (MQ=255)
                                                                 gACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGaa                                                                                                                             <  2:1219250/106‑1 (MQ=255)
                                                                 gACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGaa                                                                                                                             >  1:1219250/1‑106 (MQ=255)
                                                                  acacAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGc                                                                                <  1:1234144/150‑1 (MQ=255)
                                                                    acaGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTg                                                                                  >  1:674174/1‑146 (MQ=255)
                                                                    acaGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTg                                                                                  <  2:674174/146‑1 (MQ=255)
                                                                      aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGa                                                                                                                      >  1:73620/1‑108 (MQ=255)
                                                                      aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGa                                                                                                                      <  2:73620/108‑1 (MQ=255)
                                                                      aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAg                                                                                                      <  2:358357/124‑1 (MQ=255)
                                                                      aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAg                                                                                                      >  1:358357/1‑124 (MQ=255)
                                                                      aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTAcc                                                                            >  2:65155/1‑150 (MQ=255)
                                                                          aTCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACccgcc                                                                        >  2:1190422/1‑150 (MQ=255)
                                                                          aTCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACccgcc                                                                        <  1:515185/150‑1 (MQ=255)
                                                                             gACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTg                                                                     >  1:977448/1‑150 (MQ=255)
                                                                               cGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCt                                                                   <  1:1115297/150‑1 (MQ=255)
                                                                                  cccGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGt                                                                >  1:1218665/1‑150 (MQ=255)
                                                                                    cGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCa                                                              <  2:1393254/150‑1 (MQ=255)
                                                                                    cGGCACAGCATCGGTGGGTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCGGTGGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCa                                                              <  1:450275/150‑1 (MQ=255)
                                                                                     ggCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAg                                                             >  1:808125/1‑150 (MQ=255)
                                                                                                ggtggtTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTcc                                                  <  2:977448/150‑2 (MQ=255)
                                                                                                     ttCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGtt                                                                                                          <  1:126775/89‑1 (MQ=255)
                                                                                                     ttCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGtt                                                                                                          >  2:126775/1‑89 (MQ=255)
                                                                                                     ttCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGt                                             <  2:819304/150‑1 (MQ=255)
                                                                                                     ttCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGt                                             >  2:817469/1‑150 (MQ=255)
                                                                                                      tCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTg                                            >  2:67748/1‑150 (MQ=255)
                                                                                                       cGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGCATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGc                                           >  1:1373096/1‑150 (MQ=255)
                                                                                                         gCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAg                                                                                                                                  <  2:363256/61‑1 (MQ=255)
                                                                                                         gCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAg                                                                                                                                  >  1:363256/1‑61 (MQ=255)
                                                                                                            caGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGt                                      >  1:1333152/1‑150 (MQ=255)
                                                                                                                       ggAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGt                           >  1:16579/1‑150 (MQ=255)
                                                                                                                        gggcgGGTTCACCCGCGAACACCGGCGCAGCCGGTGCCATgttccgcgttgcctgcttcgcgggtgaacccgctcctacaaggtcagcggatagctgacgatcaggcgcgtctggtcgaga                                                       >  1:529558‑M1/3‑40 (MQ=255)
                                                                                                                        gggcgGGTTCACCCGCGAACACCGGCGCAGCCGGTGCCATgttccgcgttgcctgcttcgcgggtgaacccgctcctacaaggtcagcggatagctgacgatcaggcgcgtctggtcgaga                                                       <  2:529558‑M1/119‑82 (MQ=255)
                                                                                                                            gggTTCACCCGCGAACACCGGCGCAGCCGGTGCCATgttccgcgttgcctgcttcgcgggtgaacccgcccctacaaggtcagcggatagctgacgatcag                                                                       >  2:1064280‑M1/1‑36 (MQ=255)
                                                                                                                            gggTTCACCCGCGAACACCGGCGCAGCCGGTGCCATgttccgcgttgcctgcttcgcgggtgaacccgcccctacaaggtcagcggatagctgacgatcag                                                                       <  1:1064280‑M1/101‑66 (MQ=255)
                                                                                                                                  aCCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTg                <  1:1140163/150‑1 (MQ=255)
                                                                                                                                   cccGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTgg               <  2:1333152/150‑1 (MQ=255)
                                                                                                                                    ccGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGtgctgc                                 <  1:833493/131‑1 (MQ=255)
                                                                                                                                    ccGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGtgctgc                                 >  2:833493/1‑131 (MQ=255)
                                                                                                                                    ccGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGt              >  1:359793/1‑150 (MQ=255)
                                                                                                                                    ccGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGt              >  1:1301950/1‑150 (MQ=255)
                                                                                                                                      gcgAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCg            <  2:16579/150‑1 (MQ=255)
                                                                                                                                            acCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCGGCCAGt      >  1:451013/1‑150 (MQ=255)
                                                                                                                                                gCAGAGCCGGTGCCATCCCCAGAGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCGCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCGGCCAGTTCgg  <  2:197721/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCGGCCAGTTCGG  >  NC_002947/2406861‑2407154

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