Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,821:1 +G 6.9% coding (494/1503 nt) nhaB → Na(+)/H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,541,8211.G 6.9% 525.4 / 32.3 130coding (494/1503 nt)nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base . (72/49);  new base G (2/7);  total (74/56)
Fisher's exact test for biased strand distribution p-value = 3.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

Only 100 of 140 total aligned reads displayed.

GTTCCTGTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCC  >  NC_002947/4541672‑4541970
                                                                                                                                                      |                                                                                                                                                     
gttcctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGgt                                                                                                                                                        <  1:270647/150‑1 (MQ=255)
gttcctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCGATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGgt                                                                                                                                                        <  2:216372/150‑1 (MQ=255)
 ttcctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGgtg                                                                                                                                                      >  2:731427/1‑150 (MQ=255)
  tcctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGgtgt                                                                                                                                                     <  1:2232998/150‑1 (MQ=255)
  tcctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGgtgt                                                                                                                                                     >  1:1420443/1‑150 (MQ=255)
   cctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTa                                                                                                                                                    >  2:1593260/1‑150 (MQ=255)
    ctgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTac                                                                                                                                                   >  1:931524/1‑150 (MQ=255)
      gttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGgt                                                                                                                                                        <  1:1684971/144‑1 (MQ=255)
      gttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTacca                                                                                                                                                 >  1:2101021/1‑150 (MQ=255)
            aCGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgt                                                                                                                                           >  1:2448675/1‑150 (MQ=255)
            aCGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgt                                                                                                                                           >  1:1364980/1‑150 (MQ=255)
            aCGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgt                                                                                                                                           <  2:2727495/150‑1 (MQ=255)
             cGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgta                                                                                                                                          <  2:1611908/150‑2 (MQ=255)
             cGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtc                                                                                                                                          <  1:2545483/150‑1 (MQ=255)
             cGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCACTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtc                                                                                                                                          <  2:1953170/150‑1 (MQ=255)
               cgttttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgc                                                                                                                                        >  1:2745136/3‑150 (MQ=255)
                aTTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgtagaa                                                                                                                                       <  1:1525840/150‑5 (MQ=255)
                aTTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgc                                                                                                                                             >  2:792223/1‑144 (MQ=255)
                aTTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgc                                                                                                                                             <  1:792223/144‑1 (MQ=255)
                 ttttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgtagaag                                                                                                                                      <  1:731427/150‑6 (MQ=255)
                   ttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgtagaagcc                                                                                                                                    <  1:1313929/150‑8 (MQ=255)
                      cTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                                 <  2:2101021/150‑1 (MQ=255)
                      cTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                                 <  1:2883984/150‑1 (MQ=255)
                       tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                                >  2:2571194/1‑150 (MQ=255)
                       tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                                >  2:1293275/1‑150 (MQ=255)
                       tCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                                >  2:2722659/1‑150 (MQ=255)
                        cGGGGTGCGTGCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTGCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGc                                                                                                                               <  2:171766/150‑1 (MQ=255)
                          gggTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCa                                                                                                                             <  2:1132139/150‑1 (MQ=255)
                                 ttCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgg                                                                                                                                  <  1:2935094/138‑1 (MQ=255)
                                 ttCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgg                                                                                                                                  >  2:2935094/1‑138 (MQ=255)
                                 ttCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGCCGTGTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcgcg                                                                                                                      >  1:3010883/1‑150 (MQ=255)
                                  tCCAAGGCCATTCTGGCGCTGCTGTTTTTCGTACTGTCAGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATTATCAGCGCTGCGGGCTGCTTCTACGCTGTGTACCACCGCGACGCGTCCGCGGCAACTCggcgc                                                                                                                     >  2:71027/1‑150 (MQ=255)
                                    cAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGa                                                                                                                   >  1:1133883/1‑150 (MQ=255)
                                     aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATTATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                  >  1:2060203/1‑150 (MQ=255)
                                     aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                  <  1:2722659/150‑1 (MQ=255)
                                      aGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgtagaa                                                                                                                                       <  2:1736987/128‑5 (MQ=255)
                                      aGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgcgtagaa                                                                                                                                       >  1:1736987/1‑124 (MQ=255)
                                            ttCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGAcgacg                                                                                                           >  1:1287302/1‑146 (MQ=255)
                                            ttCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACa                                                                                                              >  2:963124/1‑147 (MQ=255)
                                              cTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCa                                                                                                         <  2:2448675/150‑1 (MQ=255)
                                                 gcgcTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTg                                                                                                      <  2:2745136/150‑1 (MQ=255)
                                                 gcgcTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATCCCCTGACCGTCACGGCGGTGATTATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTg                                                                                                      <  2:2060203/150‑1 (MQ=255)
                                                   gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATTCGCGCGAGGACAGCGCACTGga                                                                                                    <  1:577981/150‑1 (MQ=255)
                                                     tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTacca                                                                                                                                                 <  2:2241539/103‑1 (MQ=255)
                                                     tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTacca                                                                                                                                                 >  1:2241539/1‑103 (MQ=255)
                                                     tgctgTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGgaga                                                                                                  >  2:2425632/1‑150 (MQ=255)
                                                                 tACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAg                                                                                      >  1:1515950/1‑150 (MQ=255)
                                                                      tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                                 <  2:2069810/102‑1 (MQ=255)
                                                                      tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGcc                                                                                 >  1:1862363/1‑149 (MQ=255)
                                                                         gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCAc                                                                              <  1:207132/150‑1 (MQ=255)
                                                                           gTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCt                                                                            >  2:221116/1‑150 (MQ=255)
                                                                           gTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCt                                                                            >  1:1933206/1‑150 (MQ=255)
                                                                                 aTCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCg                                                                      <  2:654304/150‑1 (MQ=255)
                                                                                  tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcgcgc                                                                                                                     >  2:1333989/1‑102 (MQ=255)
                                                                                  tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCgg                                                                     >  2:1803890/1‑150 (MQ=255)
                                                                                  tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCgg                                                                     >  1:2754851/1‑150 (MQ=255)
                                                                                  tCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGAAAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCgg                                                                     >  1:1272556/1‑150 (MQ=255)
                                                                                   ccGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCggg                                                                    >  2:2289183/1‑150 (MQ=255)
                                                                                    cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTGCCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAg                                                                                          >  1:2832297/1‑127 (MQ=255)
                                                                                    cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAg                                                                                          <  2:2832297/127‑1 (MQ=255)
                                                                                    cgcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGa                                                                   >  1:919074/1‑150 (MQ=255)
                                                                                     gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                  >  2:1601129/1‑150 (MQ=255)
                                                                                     gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                                  <  1:2425632/150‑1 (MQ=255)
                                                                                       gTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGa                                                                <  1:2571194/150‑1 (MQ=255)
                                                                                         ttCTCGATGCCCGGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGAcc                                                              <  2:1255287/150‑1 (MQ=255)
                                                                                              gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGAc                                                         <  1:2289183/150‑1 (MQ=255)
                                                                                              gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGAc                                                         <  2:2070701/150‑1 (MQ=255)
                                                                                               aTGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTTCGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGAcc                                                        >  2:2221898/1‑150 (MQ=255)
                                                                                               aTGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCg                                                                                   >  1:1774439/1‑123 (MQ=255)
                                                                                               aTGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGAcc                                                        >  2:3043283/1‑150 (MQ=255)
                                                                                                 gccgtGACGGTCAGGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGAc                                                                                             <  2:292666/107‑1 (MQ=255)
                                                                                                   ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGtt                                                    >  2:1876140/1‑150 (MQ=255)
                                                                                                       aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                                >  1:2791089/1‑150 (MQ=255)
                                                                                                        ccGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                  <  2:65132/83‑1 (MQ=255)
                                                                                                        ccGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAg                                                                                                                  >  1:65132/1‑83 (MQ=255)
                                                                                                          gTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGAcc                                                              >  2:1848908/1‑133 (MQ=255)
                                                                                                                 cggTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGc                                      <  2:1364980/150‑1 (MQ=255)
                                                                                                                    tGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                   <  2:2023814/150‑1 (MQ=255)
                                                                                                                         atcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGc                                 >  1:1520275/1‑147 (MQ=255)
                                                                                                                         atcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCACCTTCCTGCGTAGc                                 <  2:1520275/147‑1 (MQ=255)
                                                                                                                          tcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                             >  1:1930068/1‑150 (MQ=255)
                                                                                                                            aGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGCCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACAAGGAAGAGCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgctg                           >  1:2503759/1‑150 (MQ=255)
                                                                                                                             gcgcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTTGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGa                          >  1:478026/1‑150 (MQ=255)
                                                                                                                               gcTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATg                        >  1:1199853/1‑150 (MQ=255)
                                                                                                                                     gTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgc                  >  2:2868026/1‑150 (MQ=255)
                                                                                                                                          cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt             <  1:1293275/150‑1 (MQ=255)
                                                                                                                                          cTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt             >  2:605817/1‑150 (MQ=255)
                                                                                                                                          cTTCTACGCGGTGTACCACCGCGCAGCGCCCGGGGCAAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgt             <  2:1272556/150‑1 (MQ=255)
                                                                                                                                           ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGGGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg            >  2:2210422/1‑150 (MQ=255)
                                                                                                                                           ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg            >  1:1549345/1‑150 (MQ=255)
                                                                                                                                           ttCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCCTAGCCTGCTGATGCATGTCGCtgtg            >  1:72495/1‑150 (MQ=255)
                                                                                                                                               aCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGAt                         <  2:497260/133‑1 (MQ=255)
                                                                                                                                               aCGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTa        >  2:1074391/1‑150 (MQ=255)
                                                                                                                                               aCGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTa        >  1:111960/1‑150 (MQ=255)
                                                                                                                                                 gcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAg                                                      <  2:1123068/102‑1 (MQ=255)
                                                                                                                                                 gcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAg                                                      >  1:1123068/1‑102 (MQ=255)
                                                                                                                                                 gcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTAcc      >  1:1853809/1‑150 (MQ=255)
                                                                                                                                                   ggTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                   >  2:1558062/1‑119 (MQ=255)
                                                                                                                                                   ggTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                   <  1:1558062/119‑1 (MQ=255)
                                                                                                                                                     tgtACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGccc  <  2:931524/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTTCCTGTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCGCCC  >  NC_002947/4541672‑4541970

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: