Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,187,974 T→C 7.9% intergenic (+79/‑131) PP_3683 → / → PP_3684 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,187,9740TC 7.9% 180.5 / 3.9 63intergenic (+79/‑131)PP_3683/PP_3684hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base C (2/3);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTCGCGT  >  NC_002947/4187828‑4188122
                                                                                                                                                   |                                                                                                                                                    
tAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCt                                                                                                                                                      >  6:200144/1‑147 (MQ=255)
tAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCt                                                                                                                                                      >  7:83287/1‑147 (MQ=255)
  gCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGc                                                                                                                                                  <  2:330349/149‑1 (MQ=255)
   cGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCg                                                                                                                                                 <  2:181741/149‑1 (MQ=255)
     ttggttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGc                                                                                                                                               >  4:225355/1‑149 (MQ=255)
        gttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTg                                                                                                                                            <  5:172019/149‑1 (MQ=255)
        gttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTg                                                                                                                                            >  3:67788/1‑149 (MQ=255)
        gttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTg                                                                                                                                            <  7:75787/149‑1 (MQ=255)
         ttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTgg                                                                                                                                           >  7:14138/1‑149 (MQ=255)
         ttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCTTgg                                                                                                                                           <  1:95703/149‑1 (MQ=255)
         ttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCAGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTgg                                                                                                                                           <  4:23556/149‑1 (MQ=255)
                cTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAAc                                                                                                                                    >  7:182297/1‑149 (MQ=255)
                cTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAAc                                                                                                                                    >  6:234115/1‑149 (MQ=255)
                 tCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACg                                                                                                                                   >  8:140282/1‑149 (MQ=255)
                   aTCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACggg                                                                                                                                 <  5:146850/149‑1 (MQ=255)
                      ggTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCtt                                                                                                                              <  5:234115/149‑1 (MQ=255)
                        tAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGTAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTc                                                                                                                            <  2:314834/149‑1 (MQ=255)
                              ccGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGTTTTCATc                                                                                                                         <  3:225355/146‑1 (MQ=255)
                              ccGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATcc                                                                                                                        <  8:274565/147‑1 (MQ=255)
                              ccGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCt                                                                                                                       <  4:248572/148‑1 (MQ=255)
                              ccGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTc                                                                                                                      <  7:140282/149‑1 (MQ=255)
                                gATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTcgc                                                                                                                    >  8:188918/1‑149 (MQ=255)
                                gATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCTTg                                                                                                                                            <  4:205930/125‑1 (MQ=255)
                                gATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGCAAGCGGCTTg                                                                                                                                            >  3:205930/1‑125 (MQ=255)
                                    aCGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCt                                                                                                                >  2:102222/1‑149 (MQ=255)
                                      gATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTcc                                                                                                              >  4:160826/1‑149 (MQ=255)
                                          ccAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGc                                                                                                           <  8:14138/148‑1 (MQ=255)
                                              ggATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAgg                                                                                                      <  8:83287/149‑1 (MQ=255)
                                                   aaGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTgc                                                                                                 <  7:31248/149‑1 (MQ=255)
                                                   aaGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTg                                                                                                  <  8:254768/148‑1 (MQ=255)
                                                     gACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTgcgc                                                                                               >  5:291737/1‑149 (MQ=255)
                                                              tCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTTTCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCt                                                                                       <  1:102222/149‑1 (MQ=255)
                                                              tCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCtt                                                                                      <  8:182297/149‑1 (MQ=255)
                                                                cTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGc                                                                                    >  6:340483/1‑149 (MQ=255)
                                                                   gCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGcc                                                                                   <  2:154810/147‑1 (MQ=255)
                                                                   gCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCa                                                                                  >  6:191584/1‑148 (MQ=255)
                                                                         gCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCaa                                                                            <  8:34710/148‑1 (MQ=255)
                                                                           tttCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAg                                                                         >  3:233556/1‑149 (MQ=255)
                                                                               tcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTc                                                                      >  4:274726/1‑148 (MQ=255)
                                                                                 aGCTCTCTGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCAtc                                                                   <  8:180986/149‑1 (MQ=255)
                                                                                 aGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCAtc                                                                   >  5:89822/1‑149 (MQ=255)
                                                                                 aGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCAtc                                                                   >  7:136460/1‑149 (MQ=255)
                                                                                        ggCGTTACCTCATGGCGCGCCCCTTGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGa                                                             <  5:200144/148‑1 (MQ=255)
                                                                                         gCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATa                                                           <  7:270522/149‑1 (MQ=255)
                                                                                           gTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGc                                                         <  6:89822/149‑1 (MQ=255)
                                                                                            ttACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCt                                                        <  8:136460/149‑1 (MQ=255)
                                                                                                cTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTcc                                                      >  2:171638/1‑147 (MQ=255)
                                                                                                 tCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCaaaa                                                                          <  4:35362/126‑1 (MQ=255)
                                                                                                 tCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCaaaa                                                                          >  3:35362/1‑126 (MQ=255)
                                                                                                   aTGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTa                                                 >  6:6224/1‑149 (MQ=255)
                                                                                                   aTGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACAGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTa                                                 >  3:283466/1‑149 (MQ=255)
                                                                                                      gcgcgcCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACaa                                              <  7:188918/149‑1 (MQ=255)
                                                                                                      gcgcgcCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACaa                                              <  4:233556/149‑1 (MQ=255)
                                                                                                      gcgcgcCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACa                                               <  5:6224/148‑1 (MQ=255)
                                                                                                                 gCGGGGCCGGCTGGCGTGAACCAAGGCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTc                                   <  3:274726/149‑1 (MQ=255)
                                                                                                                 gCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACAGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTc                                   <  4:283466/149‑1 (MQ=255)
                                                                                                                          gCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGAt                          >  5:112532/1‑149 (MQ=255)
                                                                                                                               cGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGc                     <  3:212214/149‑1 (MQ=255)
                                                                                                                                gTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGcc                    >  5:95979/1‑149 (MQ=255)
                                                                                                                                   aaCCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCtt                  <  5:191584/148‑1 (MQ=255)
                                                                                                                                    gccAAGCCGCTTAGCAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGt                                                  <  7:186033/114‑1 (MQ=255)
                                                                                                                                    gccAAGCCGCTTAGCAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGt                                                  >  8:186033/2‑115 (MQ=255)
                                                                                                                                    aCCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTgc                >  6:316982/1‑149 (MQ=255)
                                                                                                                                               ttGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTc      >  6:195305/1‑149 (MQ=255)
                                                                                                                                                 gCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTcgc    >  8:3441/1‑149 (MQ=255)
                                                                                                                                                   tAAGCGGCTTGGCCCTACCGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTCGCGt  >  5:262041/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTGCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTCGCGT  >  NC_002947/4187828‑4188122

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