Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,293,255 T→G 5.8% A286A (GCT→GCG)  PP_1128 → OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,293,2550TG 5.8% 207.5 / 2.1 65A286A (GCT→GCG) PP_1128OmpA family protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base G (2/2);  total (34/31)
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

GCGCCGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAG  >  NC_002947/1293107‑1293400
                                                                                                                                                    |                                                                                                                                                   
gcgcCGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCt                                                                                                                                                     <  3:384394/149‑1 (MQ=255)
   ccGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGc                                                                                                                                                   >  5:410822/1‑148 (MQ=255)
     gAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCaa                                                                                                                                                <  4:228869/149‑1 (MQ=255)
      aGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCGCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGccttg                                                                                                                                               >  4:323060/1‑146 (MQ=255)
       gggCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑gg                                                                                                                                             <  4:110174/149‑1 (MQ=255)
          ccGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCg                                                                                                                                          >  7:14162/1‑149 (MQ=255)
             cgcCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtg                                                                                                                                       <  2:12735/149‑1 (MQ=255)
             cgcCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtg                                                                                                                                       <  2:6070/149‑1 (MQ=255)
               cccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGACCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtggt                                                                                                                                     <  1:175567/149‑1 (MQ=255)
                ccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGCCTT‑GGCAGCtggtg                                                                                                                                    <  3:323060/149‑1 (MQ=255)
                 cTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtggtgg                                                                                                                                   <  4:149877/149‑1 (MQ=255)
                  tGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtggtgg                                                                                                                                   >  7:339391/1‑148 (MQ=255)
                      ccgccgCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAAc                                                                                                                              <  6:410822/149‑1 (MQ=255)
                          cgcgTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGc                                                                                                                          <  8:14162/149‑1 (MQ=255)
                                gAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTc                                                                                                                      >  7:123573/1‑147 (MQ=255)
                                gAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTc                                                                                                                      >  4:38691/1‑147 (MQ=255)
                                  aaTGTTGCTGACCCCGCAAAACACCATGGTCGTGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAgg                                                                                                                  >  2:217763/1‑149 (MQ=255)
                                        gCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGc                                                                                                                                            >  8:29934/1‑117 (MQ=255)
                                         cTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTgcgg                                                                                                            <  2:103389/148‑1 (MQ=255)
                                            aCCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCga                                                                                                         <  8:339391/148‑1 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTTGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGCCTT‑GGCAGCtg                                                                                                                                       <  4:369948/116‑1 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtgg                                                                                                                                      >  4:293045/1‑117 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCtgg                                                                                                                                      <  3:293045/117‑1 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGcc                                                                                                                         <  5:87509/130‑1 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGcc                                                                                                                         >  6:87509/1‑130 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaag                                                                                                       >  6:158736/1‑148 (MQ=255)
                                              cccGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGCCTT‑GGCAGCtg                                                                                                                                       >  3:369948/1‑116 (MQ=255)
                                                   aaaaCACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAAc                                                                                                   >  4:396182/1‑147 (MQ=255)
                                                      acacCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTa                                                                                              <  3:396182/149‑1 (MQ=255)
                                                              gTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAgg                                                                                       <  2:384833/148‑1 (MQ=255)
                                                              gTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGc                                                                                      >  4:147082/1‑149 (MQ=255)
                                                                   gctgctGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCaaaa                                                                                 >  5:314805/1‑149 (MQ=255)
                                                                        tgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAACGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGaaa                                                                                                    <  1:285479/126‑1 (MQ=255)
                                                                        tgGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGaaa                                                                                                    >  2:285479/1‑125 (MQ=255)
                                                                         ggCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTCCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTCCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTg                                                                           >  2:78818/1‑149 (MQ=255)
                                                                          gCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaa                                                                                                        <  3:97174/119‑1 (MQ=255)
                                                                          gCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCgaa                                                                                                        >  4:97174/1‑119 (MQ=255)
                                                                          gCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCt                                                                            >  6:32916/1‑147 (MQ=255)
                                                                          gCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCt                                                                            <  5:32916/147‑1 (MQ=255)
                                                                             aaaTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCGGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCg                                                                         >  6:375498/1‑147 (MQ=255)
                                                                                tACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGc                                                                              <  6:141108/139‑1 (MQ=255)
                                                                                tACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGc                                                                              >  5:141108/1‑139 (MQ=255)
                                                                                   cagcagGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAg                                                                                        >  4:140629/1‑126 (MQ=255)
                                                                                   cagcagGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAg                                                                                        <  3:140629/126‑1 (MQ=255)
                                                                                         gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGcc                                                                                                                         <  1:38195/87‑1 (MQ=255)
                                                                                         gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGcc                                                                                                                         >  2:38195/1‑87 (MQ=255)
                                                                                           cgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGc                                                                                            <  2:73522/114‑1 (MQ=255)
                                                                                                    aCCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGcc                                                             >  8:286585/1‑136 (MQ=255)
                                                                                                    aCCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGcc                                                             <  7:286585/136‑1 (MQ=255)
                                                                                                    aCCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTg                                               >  2:198443/1‑149 (MQ=255)
                                                                                                      ccGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTgcg                                             >  6:345021/1‑149 (MQ=255)
                                                                                                       cGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTgcgc                                            <  3:147082/149‑1 (MQ=255)
                                                                                                        gCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGcc                                           >  2:318179/1‑149 (MQ=255)
                                                                                                          cccGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCgg                                         >  1:95059/1‑149 (MQ=255)
                                                                                                                  cGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCt                                 >  3:174029/1‑149 (MQ=255)
                                                                                                                     tggtggCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTgcc                              <  6:314805/149‑1 (MQ=255)
                                                                                                                        tggCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAACCCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTgccgc                            >  7:286573/1‑148 (MQ=255)
                                                                                                                                  ggATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCTAAGACTGCAGcc                                                             >  4:240230/1‑106 (MQ=255)
                                                                                                                                  ggATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGcc                                                             <  3:240230/106‑1 (MQ=255)
                                                                                                                                    aTGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGcc               <  1:217763/149‑1 (MQ=255)
                                                                                                                                         aaGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGt             >  3:229802/1‑146 (MQ=255)
                                                                                                                                         aaGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGt             <  4:229802/146‑1 (MQ=255)
                                                                                                                                             cACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACaa       <  1:78818/149‑1 (MQ=255)
                                                                                                                                                cAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTg                                                                           <  2:241924/78‑1 (MQ=255)
                                                                                                                                                cAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTg                                                                           >  1:241924/1‑78 (MQ=255)
                                                                                                                                                 aGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAg  >  1:349157/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCCGAGGGCCGCGCCCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAA‑GGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAG‑CTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAG  >  NC_002947/1293107‑1293400

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