Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,407,007 A→G 5.5% intergenic (+219/+35) 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,0070AG 5.5% 230.8 / 2.2 76intergenic (+219/+35)PP_2110/PP_2111ribosomal large subunit pseudouridine synthase/aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (31/40);  new base G (2/2);  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 = 6.81e-01

CTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCGGCCAGTTCGG  >  NC_002947/2406858‑2407154
                                                                                                                                                     |                                                                                                                                                   
ctgcCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCa                                                                                                                                                     <  1:563155/150‑1 (MQ=255)
ctgcCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCa                                                                                                                                                     <  1:941508/150‑1 (MQ=255)
   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                                                                                                                                >  1:488456/1‑116 (MQ=255)
                                                       gcATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCg                                                                                                                                <  2:488456/116‑1 (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                                                                                                                      <  2:73620/108‑1 (MQ=255)
                                                                         aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGa                                                                                                                      >  1:73620/1‑108 (MQ=255)
                                                                         aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAg                                                                                                      >  1:358357/1‑124 (MQ=255)
                                                                         aGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAg                                                                                                      <  2:358357/124‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
CTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACATGGACCGGGCTGGACGCGCATCAGTGGTTCGACACAGGCATCGACGGCCCGGCACAGCATCGGTGGTTCGGCACAGCCGTTGTAGGAGCGGGTTCACCCGCGAACACCGGCAGAGCCGGTGCCATCCCCAGCGGAAGTCTCGAGCTTCAGCGGTTGAAGCGCTCCACCAGCGAATACTGGCTACCCGCCGTGCTGGTCAGTTCTTCGCTCAGCAGTGCCGAGTGCTGCGCCTGTTCGGCGGTCTGGTCGGCCAGTTCGG  >  NC_002947/2406858‑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: