Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,177 C→T 8.7% P160L (CCC→CTC)  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,1770CT 8.7% 243.1 / 10.4 89P160L (CCC→CTC) PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (38/43);  new base T (3/5);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912028‑3912322
                                                                                                                                                     |                                                                                                                                                 
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                                         <  1:1064552/150‑1 (MQ=255)
                                                                                                          aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         >  2:1963622/1‑150 (MQ=255)
                                                                                                          aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  2:15009/150‑1 (MQ=255)
                                                                                                          aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  1:1284572/150‑1 (MQ=255)
                                                                                                          aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCAATGCTGATCAGATGCACCGTGTTTCGAATCGGCCGGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  1:1498083/150‑1 (MQ=255)
                                                                                                            cTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCa                                                                                                                      <  1:1616788/71‑1 (MQ=255)
                                                                                                            cTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCa                                                                                                                      >  2:1616788/1‑71 (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                                                                             >  1:1275625/1‑94 (MQ=255)
                                                                                                                              gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:1275625/94‑1 (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                                                                             <  2:258117/81‑1 (MQ=255)
                                                                                                                                         agcGCATCGGGCTCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:258117/3‑83 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctaca                                                                                  >  2:1855087/1‑77 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               <  2:1598122/96‑1 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               >  1:1598122/1‑96 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAg                                <  1:1666643/127‑1 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAg                                >  2:1666643/1‑127 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  1:350407/1‑150 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  1:1805955/1‑150 (MQ=255)
                                                                                                                                          ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  2:1149685/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)
                                                                                                                                                    cccGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCgtgt                                                                                                          >  1:1840302/1‑43 (MQ=255)
                                                                                                                                                    cccGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCgtgt                                                                                                          <  2:1840302/43‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912028‑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: