Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,269 G→C 10.2% intergenic (+45/+97) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,2690GC10.2% 233.2 / 15.8 79intergenic (+45/+97)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (31/39);  new base C (4/4);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GCTGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGG  >  NC_002947/1707122‑1707418
                                                                                                                                                   |                                                                                                                                                      
gctgTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAAACGCAAGCGCGCGGTGCGCTGCATGGGAGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCGGCGAGATCGAGGGCcgcgcg                                                                                                                                                      >  2:1024164/1‑150 (MQ=255)
    ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCccgcgcg                                                                                                                                                  >  2:817939/1‑150 (MQ=255)
      cGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGCAAACACCCCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGggcggc                                                                                                                                                <  2:1615265/150‑1 (MQ=255)
        ccacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                              <  1:817939/149‑1 (MQ=255)
          acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            <  1:1049405/150‑1 (MQ=255)
          acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            >  2:1049405/1‑150 (MQ=255)
            ccggcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          >  1:2332966/4‑150 (MQ=255)
            ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          <  2:2332966/150‑1 (MQ=255)
            ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          >  2:1259415/1‑150 (MQ=255)
             cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         >  1:691363/1‑150 (MQ=255)
              agcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                        <  1:2080121/150‑1 (MQ=255)
               gcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                       <  2:1224585/150‑1 (MQ=255)
                 agAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     <  1:935749/150‑1 (MQ=255)
                   aaaGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGAGCGAGCGCCGCGCGGGCGGCGCTCGCTCTcaca                                                                                                                                   >  2:1874354/1‑150 (MQ=255)
                        ccTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCcgcgc                                                                                                                                                       <  2:1167897/125‑1 (MQ=255)
                        ccTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCcgcgc                                                                                                                                                       >  1:1167897/1‑125 (MQ=255)
                           gACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGcac                                                                                                                           <  1:133066/150‑1 (MQ=255)
                            aCCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGcaca                                                                                                                          >  2:881906/1‑150 (MQ=255)
                              ccgCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaac                                                                                                                        >  1:1673802/3‑150 (MQ=255)
                                    gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCcgcgcg                                                                                                                                                      >  1:1095429/1‑114 (MQ=255)
                                    gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCcgcgcg                                                                                                                                                      <  2:1095429/114‑1 (MQ=255)
                                    gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCCCACGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATc                                                                                                                  <  1:645236/150‑1 (MQ=255)
                                        gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCg                                                                                                              >  1:1761207/1‑150 (MQ=255)
                                        gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCg                                                                                                              >  2:2148755/1‑150 (MQ=255)
                                           gACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                           >  2:1517187/1‑150 (MQ=255)
                                             cAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaa                                                                                                         >  1:1851855/1‑150 (MQ=255)
                                              ccccGAAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCTGCACAACATTATCGGCGAACaaa                                                                                                        <  1:2238499/148‑1 (MQ=255)
                                                ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                   >  1:1128828/1‑121 (MQ=255)
                                                ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                   <  2:1128828/121‑1 (MQ=255)
                                                     aGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGc                                                                                                 <  2:691363/150‑1 (MQ=255)
                                                      gcgcgcGGTGCGCTGAATGGGAGAGCAAACAGCGAGATACATACATTGAGGAGTGGGCGCATCATGGTGGGCGCCCGCGAGACCGAGCGCCGCGCAGGCCGCACTCGATCCAACAAGCTTCACAACATTATCGGCGAACAAATAAAATcc                                                                                                <  2:1673802/150‑1 (MQ=255)
                                                                  cTGCATGGGCCAGGCAAAACACCGCAAACTGCAGTCAGCTATCGGACTAGAATTTTTGCCTCTCGTGATATCGTGTGTCTCTCGGGCTTCGCTCGTTCTCACAGGTGTCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATcg                                                                                   <  1:2089267/150‑1 (MQ=255)
                                                                   tGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATcgc                                                                                  <  1:112643/150‑1 (MQ=255)
                                                                        gggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGc                                                                             <  1:1259415/150‑1 (MQ=255)
                                                                               gAAACACCGCGGCAGGTGCAGTTAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                      >  2:1668981/1‑150 (MQ=255)
                                                                               gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGCttttctt                                                                      <  2:2059438/150‑6 (MQ=38)
                                                                               gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                      <  1:665698/150‑1 (MQ=255)
                                                                               gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGACTGTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                      <  1:1874354/150‑1 (MQ=255)
                                                                                aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCtt                                                                     <  2:1851855/150‑1 (MQ=255)
                                                                                 aaCACCGCGGCAGGTGCAGTTAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCttt                                                                    <  1:1668981/150‑1 (MQ=255)
                                                                                  acacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                   >  2:1103933/1‑150 (MQ=255)
                                                                                    acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                      <  1:389044/145‑1 (MQ=255)
                                                                                    acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                      >  2:389044/1‑145 (MQ=255)
                                                                                    acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                 >  2:645742/1‑150 (MQ=255)
                                                                                     ccGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggc                                                                                                                                                   <  2:1724816/68‑1 (MQ=255)
                                                                                     ccGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggc                                                                                                                                                   >  1:1724816/1‑68 (MQ=255)
                                                                                       gcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTa                                                              <  1:2148755/150‑1 (MQ=255)
                                                                                              gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCggg                                                                               <  2:2197713/126‑1 (MQ=255)
                                                                                              gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCggg                                                                               >  1:2197713/1‑126 (MQ=255)
                                                                                                  aGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAg                                                   >  2:823955/1‑150 (MQ=255)
                                                                                                   gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                  >  1:1899315/1‑150 (MQ=255)
                                                                                                    tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTCTCAGCTAGCCAGCCTTAgcg                                                 >  1:1304931/1‑150 (MQ=255)
                                                                                                      aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                               <  2:436907/150‑1 (MQ=255)
                                                                                                      aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                               >  1:1597687/1‑150 (MQ=255)
                                                                                                       ctcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                               <  2:1698703‑M2/60‑1 (MQ=255)
                                                                                                       ctcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                               >  1:1698703‑M2/27‑86 (MQ=255)
                                                                                                        tcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                   <  2:2319903‑M2/40‑1 (MQ=255)
                                                                                                        tcgcgggcgccaaaaatgctgcgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                   >  1:2319903‑M2/26‑65 (MQ=255)
                                                                                                        gTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                >  1:284580/1‑115 (MQ=255)
                                                                                                        gTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                <  2:284580/115‑1 (MQ=255)
                                                                                                          ttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAAc                                           <  1:823955/150‑1 (MQ=255)
                                                                                                           tCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCGTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAAcc                                          >  1:1410728/1‑150 (MQ=255)
                                                                                                              GCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                       >  1:881906‑M1/1‑150 (MQ=255)
                                                                                                                  agcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGCCGCCAAAGCGATCCAGTCTTGGCGCTTAGCCTTGCTCa                                   >  2:610843/1‑150 (MQ=25)
                                                                                                                  agcaTTGTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCa                                   <  2:2281086/150‑1 (MQ=255)
                                                                                                                     aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                        >  1:2102110/1‑110 (MQ=255)
                                                                                                                     aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                        <  2:2102110/110‑1 (MQ=255)
                                                                                                                     aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                >  2:1507341/1‑150 (MQ=255)
                                                                                                                             cgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        <  1:1024164‑M2/146‑1 (MQ=255)
                                                                                                                             cgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        <  1:802076‑M2/146‑1 (MQ=255)
                                                                                                                             cgccCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        <  2:1761207‑M2/146‑1 (MQ=255)
                                                                                                                                    gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCg                 <  2:1595949/150‑1 (MQ=255)
                                                                                                                                     agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCgg                <  2:1604251/150‑1 (MQ=255)
                                                                                                                                      gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGc               <  1:1846131/150‑1 (MQ=255)
                                                                                                                                      gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGc               <  1:1507341/150‑1 (MQ=255)
                                                                                                                                        tCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGc                                                                             <  2:1865755/86‑1 (MQ=37)
                                                                                                                                        tCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGc                                                                             >  1:1865755/1‑86 (MQ=37)
                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAg                                                                                                       >  2:2286385/1‑59 (MQ=255)
                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAg                                                                                                       <  1:2286385/59‑1 (MQ=255)
                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                <  1:1375518/81‑1 (MQ=25)
                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                >  2:1375518/1‑81 (MQ=25)
                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt           >  1:1509367/1‑150 (MQ=255)
                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTCTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt          <  2:1304931/150‑1 (MQ=255)
                                                                                                                                             cgccgcCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        >  2:1621589/1‑134 (MQ=255)
                                                                                                                                             cgccgcCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        <  1:1621589/134‑1 (MQ=255)
                                                                                                                                              gccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        >  1:1013685/1‑133 (MQ=255)
                                                                                                                                              gccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                        <  2:1013685/133‑1 (MQ=255)
                                                                                                                                                   gcgggcggCGCTCGATCTCACAGGCGGCACACCATAGTTAGAGGCCGATGAAAAGCCCCGCGTCGATCGCGCGACTGGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTgg  <  1:1103933/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCTGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGG  >  NC_002947/1707122‑1707418

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