Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,174 A→G 7.7% E159G (GAG→GGG)  PP_3452 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,912,1740AG 7.7% 256.9 / 8.9 91E159G (GAG→GGG) PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (38/45);  new base G (3/4);  total (42/49)
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.56e-01

CCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912025‑3912322
                                                                                                                                                     |                                                                                                                                                    
cccGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGa                                                                                                                                                      <  2:1403538/150‑1 (MQ=255)
  cGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGc                                                                                                                                                    <  1:2380573/150‑1 (MQ=255)
   gATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGcc                                                                                                                                                   >  2:1498083/1‑150 (MQ=255)
        cccTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTTCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCACCGAGCCCGATg                                                                                                                                              >  2:775383/1‑150 (MQ=255)
         ccTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGTTCATTTTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGACCTCGAtat                                                                                                                                             <  2:1201853/150‑3 (MQ=255)
           tGCGGCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGGGCTCGATgcgg                                                                                                                                           <  1:216097/150‑2 (MQ=255)
           tGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGGGCTCGATgcgg                                                                                                                                           <  2:1863562/150‑2 (MQ=255)
           tGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATgcgc                                                                                                                                           <  1:496660/150‑1 (MQ=255)
              gcCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACAACACCCGCATCGAGCCCGATGCGCTCg                                                                                                                                        >  1:219959/1‑150 (MQ=255)
                ccGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGGGCTCGATgcgg                                                                                                                                           <  2:53473/145‑2 (MQ=255)
                ccGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGGGCTCGATgcgg                                                                                                                                           >  1:53473/1‑144 (MQ=255)
                ccGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGaa                                                                                                                                      >  1:2396838/1‑150 (MQ=255)
                ccGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGaa                                                                                                                                      <  2:1366068/150‑1 (MQ=255)
                          tggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGc                                                                                                                            <  1:211960/150‑1 (MQ=255)
                          tggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGc                                                                                                                            <  1:180831/150‑1 (MQ=255)
                          tggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACAACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGc                                                                                                                            <  2:219959/150‑1 (MQ=255)
                               cTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATc                                                                                                                       >  2:1254411/1‑150 (MQ=255)
                                  ggggTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGc                                                                                                                    >  1:1808138/1‑150 (MQ=255)
                                       gtCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGc                                                                                                               <  2:618932/150‑1 (MQ=255)
                                        tCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGcc                                                                                                              >  2:692500/1‑150 (MQ=255)
                                          cACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGGTTCTTAATCCGTTCTACCACAACCCAACCGCGCCCGATGCTCTCGAAGGCGCCATGCTGATCAGCTATGCCgt                                                                                                            >  1:1738976/1‑150 (MQ=255)
                                                   gTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAg                                                                                                   >  1:15009/1‑150 (MQ=255)
                                                    tATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGc                                                                                                  <  2:992933/150‑1 (MQ=255)
                                                     aTGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCg                                                                                                 <  2:2396838/150‑1 (MQ=255)
                                                              cGGTGCTACCGGTGCGGCGCTGGTCATTCATCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGATCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTca                                                                                        <  1:556806/150‑1 (MQ=255)
                                                                   cTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctac                                                                                   <  2:1808138/150‑1 (MQ=255)
                                                                        ggTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctacagcta                                                                              >  1:1051/1‑150 (MQ=255)
                                                                         gTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:773760/150‑1 (MQ=255)
                                                                             ggcgCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCCTGCGCTGCAGCGATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGc                                                                         <  1:1870511/150‑1 (MQ=255)
                                                                             ggcgCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATg                                                                                                                <  1:2212094/111‑1 (MQ=255)
                                                                             ggcgCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATg                                                                                                                >  2:2212094/1‑111 (MQ=255)
                                                                              gcgaTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCa                                                                        >  1:1426818/1‑150 (MQ=255)
                                                                                gcTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            >  2:2026546/1‑128 (MQ=255)
                                                                                gcTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            <  1:2026546/128‑1 (MQ=255)
                                                                                    gTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCg                                                                  >  2:1926893/1‑150 (MQ=255)
                                                                                       aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               >  1:138628/1‑150 (MQ=255)
                                                                                       aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               >  2:1064552/1‑150 (MQ=255)
                                                                                        ttcttcTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctaca                                                                                  >  2:733512/1‑130 (MQ=255)
                                                                                        ttcttcTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctaca                                                                                  <  1:733512/130‑1 (MQ=255)
                                                                                        ttcttcTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAg                                                              <  2:1426818/150‑1 (MQ=255)
                                                                                         tcttctGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               >  1:1935011/1‑148 (MQ=255)
                                                                                         tcttctGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               <  2:1935011/148‑1 (MQ=255)
                                                                                           ttctGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGcc                                                           >  2:1284572/1‑150 (MQ=255)
                                                                                                 gcCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTg                                                     >  2:2102813/1‑150 (MQ=255)
                                                                                                   attaagcatcaCCCAGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGGTCCGCAAGGCCAAGCTGca                                                   >  2:1543908/7‑150 (MQ=255)
                                                                                                       agcatcaCCCAGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGGTCCGCAAGGCCAAGCTGCACAAt                                               <  1:1543908/148‑1 (MQ=255)
                                                                                                        tcatcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCt                                                                                                      <  1:1430255/94‑1 (MQ=255)
                                                                                                        tcatcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCt                                                                                                      >  2:1430255/1‑94 (MQ=255)
                                                                                                        gcatcaCCCAGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAAtt                                              >  1:1797855/2‑150 (MQ=255)
                                                                                                        gcatcaCCCAGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAAtt                                              >  1:1813624/2‑150 (MQ=255)
                                                                                                           tcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCt                                           <  2:2064214/150‑1 (MQ=255)
                                                                                                             aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         >  2:1963622/1‑150 (MQ=255)
                                                                                                             aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  1:1064552/150‑1 (MQ=255)
                                                                                                             aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  2:15009/150‑1 (MQ=255)
                                                                                                             aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  1:1284572/150‑1 (MQ=255)
                                                                                                             aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCAATGCTGATCAGATGCACCGTGTTTCGAATCGGCCGGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  1:1498083/150‑1 (MQ=255)
                                                                                                               cTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCa                                                                                                                      >  2:1616788/1‑71 (MQ=255)
                                                                                                               cTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCa                                                                                                                      <  1:1616788/71‑1 (MQ=255)
                                                                                                                      aTGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAg                                <  1:184911/150‑1 (MQ=255)
                                                                                                                       tGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAgg                               >  1:2380625/1‑150 (MQ=255)
                                                                                                                       tGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAgg                               <  1:692500/150‑1 (MQ=255)
                                                                                                                          ttAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGgtgt                            >  2:1488365/1‑150 (MQ=255)
                                                                                                                          ttAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGgtgt                            >  2:284964/1‑150 (MQ=255)
                                                                                                                           tAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGtt                           <  1:1723724/150‑1 (MQ=255)
                                                                                                                              tCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCt                        >  2:326048/1‑150 (MQ=255)
                                                                                                                               ccGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTg                       <  1:1170423/150‑1 (MQ=255)
                                                                                                                                 gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:1275625/94‑1 (MQ=255)
                                                                                                                                 gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:1275625/1‑94 (MQ=255)
                                                                                                                                 gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGa                     <  1:775383/150‑1 (MQ=255)
                                                                                                                                  ttCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAg                    <  1:1926893/150‑1 (MQ=255)
                                                                                                                                     tACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCaa                 <  1:1963622/150‑1 (MQ=255)
                                                                                                                                       ccACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGc               <  1:326048/150‑1 (MQ=255)
                                                                                                                                        cacaCCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGcc              >  2:52923/1‑150 (MQ=255)
                                                                                                                                        cacaCCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGAAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGcc              >  2:442707/1‑150 (MQ=255)
                                                                                                                                           gagcGCATCGGGCTCGATGCGCTCGAAGGCGCCATGGTGATCAGGTATGCCGTGGTTCTAAGCGGCGTGGTCAGCTACAGCTACCTTCAGATCCGGAAGGCCAAGCTTCACAATTTCTACCTTTCCAAGGTGTTGCTGGAGCAAGCCTAc           >  1:1948566/4‑150 (MQ=255)
                                                                                                                                           ccccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTAc           <  2:1738976/149‑1 (MQ=255)
                                                                                                                                           acCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTAc           <  2:1474813/150‑1 (MQ=255)
                                                                                                                                            agcGCATCGGGCTCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:258117/3‑83 (MQ=255)
                                                                                                                                            agcGCATCGGGCTCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:258117/81‑1 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctaca                                                                                  >  2:1855087/1‑77 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               >  1:1598122/1‑96 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               <  2:1598122/96‑1 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAg                                >  2:1666643/1‑127 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAg                                <  1:1666643/127‑1 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  2:1149685/1‑150 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  1:1805955/1‑150 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  1:350407/1‑150 (MQ=255)
                                                                                                                                             ccGCATCGAGCCCGATGCGCTCGAAGGCGCCAGGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctaca                                                                                  <  1:1855087/77‑1 (MQ=255)
                                                                                                                                              cGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTACTGGAGGTATTCTACGTc        <  2:1909144/150‑1 (MQ=255)
                                                                                                                                                 aTCGAGCCAGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGAt     <  1:1042003/150‑1 (MQ=255)
                                                                                                                                                    gAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGcc  <  1:1149685/150‑1 (MQ=255)
                                                                                                                                                    gAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGcc  <  2:1813624/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912025‑3912322

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: ‑