Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,382,882 G→C 6.5% intergenic (‑93/+169) PP_3857 ← / ← PP_3858 hypothetical protein/tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,382,8820GC 6.5% 226.8 / 4.1 73intergenic (‑93/+169)PP_3857/PP_3858hypothetical protein/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (2/3);  total (37/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 = 5.65e-01

TTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGAC  >  NC_002947/4382735‑4383028
                                                                                                                                                   |                                                                                                                                                  
ttgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTg                                                                                                                                                    >  4:63169/1‑148 (MQ=255)
         cGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCa                                                                                                                                            <  5:234813/147‑1 (MQ=255)
         cGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCa                                                                                                                                            <  3:63169/147‑1 (MQ=255)
          ggACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGc                                                                                                                                          <  4:117751/148‑1 (MQ=255)
           gACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATCCAATCGccgtaga                                                                                                                                        <  3:79299/149‑6 (MQ=255)
              aacaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTTGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGcatc                                                                                                                                       <  4:55931/147‑1 (MQ=255)
              aacaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGcatcat                                                                                                                                     >  7:227184/1‑149 (MQ=255)
               acaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCatcatc                                                                                                                                    <  8:227184/149‑1 (MQ=255)
                caGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATcg                                                                                                                                   >  6:153362/1‑149 (MQ=255)
                   gCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATcggcg                                                                                                                                >  5:307568/1‑149 (MQ=255)
                     cTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCggcggc                                                                                                                              <  5:153362/149‑1 (MQ=255)
                             gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGtt                                                                                                                      >  6:48374/1‑149 (MQ=255)
                             gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATCCAATCGccgta                                                                                                                                          <  7:313167/129‑4 (MQ=255)
                             gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATCCAATCGccgta                                                                                                                                          >  8:313167/1‑126 (MQ=255)
                               ttCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCCAACGATTTGAGCAGGCTCGTTACGCGTCCCGCTCTACGGCGCTTGGATCGCCAGCAGCATCGGCGGCGCGCTGGTTa                                                                                                                    >  4:120536/1‑149 (MQ=255)
                               ttCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTTGATCGCCAGCATCATCGGCGGCGAGCTGTTTa                                                                                                                    >  4:115574/1‑149 (MQ=255)
                               ttCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGttt                                                                                                                     >  3:331234/1‑148 (MQ=255)
                                 ccAAAGCTGAAATCTCAGTCATTCACTTTTCTCCAGGCATAAAAAAACCGTCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                                                                                   <  7:287686/148‑1 (MQ=255)
                                 ccAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                                                                                   <  4:331234/148‑1 (MQ=255)
                                  cAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCa                                                                                                                 <  2:54570/149‑1 (MQ=255)
                                   aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTg                                                                                                                        >  2:298315/1‑141 (MQ=255)
                                   aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTg                                                                                                                        <  1:298315/141‑1 (MQ=255)
                                   aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCAt                                                                                                                >  5:225967/1‑149 (MQ=255)
                                      gCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGTATCGCCAGCATCATCGGAGGCGAGCTGTTTACCATTAc                                                                                                             >  4:183160/1‑149 (MQ=255)
                                        tGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACg                                                                                                            >  3:133763/1‑148 (MQ=255)
                                         gAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGc                                                                                                           >  3:84809/1‑148 (MQ=255)
                                             tatCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAAc                                                                                                      <  4:18963/149‑1 (MQ=255)
                                             tatCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAAc                                                                                                      <  4:84809/149‑1 (MQ=255)
                                               tcattcatTCACTTTTCTCCAGGCATCAAAAAACCCTCAGTTGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACTCTCTACGTCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACg                                                                                                     <  5:64946/148‑1 (MQ=255)
                                               tCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACg                                                                                                     >  5:174275/1‑148 (MQ=255)
                                                 aGTCATTCAGTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATCCAATCGccgtaga                                                                                                                                        >  5:249603/1‑106 (MQ=255)
                                                 aGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGGGATTGGATCGCCAGCATCATCGGCGGCGCGCTGTTTACCATTACGCAAGACGtga                                                                                                  >  8:315428/1‑147 (MQ=255)
                                                 aGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATCCAATCGccgtaga                                                                                                                                        <  6:249603/111‑6 (MQ=255)
                                                      ttCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccat                                                                                             <  3:267458/149‑1 (MQ=255)
                                                      ttCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccat                                                                                             >  6:162884/1‑149 (MQ=255)
                                                      ttCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAAGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatcaa                                                                                              >  7:200806/1‑146 (MQ=255)
                                                       tCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccat                                                                                             >  5:323402/1‑148 (MQ=255)
                                                        cACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTAc                                                                                                             >  4:158128/1‑131 (MQ=255)
                                                        cACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTAc                                                                                                             <  3:158128/131‑1 (MQ=255)
                                                        cACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccatc                                                                                            <  3:79460/148‑1 (MQ=255)
                                                         aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCaa                                                                                          <  5:48374/149‑1 (MQ=255)
                                                             ttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGc                                                                                       <  3:183160/148‑1 (MQ=255)
                                                               ctcCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCtt                                                                                     <  6:307568/148‑1 (MQ=255)
                                                                tcCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCa                                                                                   >  7:192489/1‑149 (MQ=255)
                                                                      cATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc                                                                             <  6:174275/149‑1 (MQ=255)
                                                                      cATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc                                                                             <  6:225967/149‑1 (MQ=255)
                                                                        tCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCt                                                                            >  1:242295/1‑148 (MQ=255)
                                                                          aaaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAt                                                                         >  6:45867/1‑149 (MQ=255)
                                                                           aaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt                                                                        >  4:111634/1‑149 (MQ=255)
                                                                            aaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc                                                                       >  8:32528/1‑149 (MQ=255)
                                                                               aCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTg                                                                    <  6:323402/149‑1 (MQ=255)
                                                                                cccGCCGTGGCGGGTTTTGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGa                                                                   <  8:200806/149‑1 (MQ=255)
                                                                                  cgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGc                                                                 <  8:192489/149‑1 (MQ=255)
                                                                                         gCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGaa                                                          <  5:162884/149‑1 (MQ=255)
                                                                                          cGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCg                                                            >  2:221660/1‑146 (MQ=255)
                                                                                                      aaaCGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGc                                             <  7:32528/149‑1 (MQ=255)
                                                                                                      aaaCGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGc                                             >  8:312810/1‑149 (MQ=255)
                                                                                                      aaaCGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGc                                             >  4:173516/1‑149 (MQ=255)
                                                                                                                    ggCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCgg                                                               >  6:30008/1‑117 (MQ=255)
                                                                                                                    ggCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCgg                                                               <  5:30008/117‑1 (MQ=255)
                                                                                                                    ggCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGaa                               <  7:312810/149‑1 (MQ=255)
                                                                                                                       tCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACAt                            >  7:46708/1‑149 (MQ=255)
                                                                                                                         gTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCg                          >  3:268526/1‑149 (MQ=255)
                                                                                                                             cgcgTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAAt                      <  3:173516/149‑1 (MQ=255)
                                                                                                                                gTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGga                   <  1:221660/149‑1 (MQ=255)
                                                                                                                                 tCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGgaa                  <  5:45867/149‑1 (MQ=255)
                                                                                                                                       tctACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                                                                                <  4:141966/81‑1 (MQ=255)
                                                                                                                                       tctACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                                                                                >  3:141966/1‑81 (MQ=255)
                                                                                                                                           cGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGt        >  1:209790/1‑149 (MQ=255)
                                                                                                                                           cGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGt        >  5:234882/1‑149 (MQ=255)
                                                                                                                                            ggCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGTc       >  3:311314/1‑149 (MQ=255)
                                                                                                                                            ggCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGTc       >  6:33200/1‑149 (MQ=255)
                                                                                                                                                aTTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGa   <  7:189073/149‑1 (MQ=255)
                                                                                                                                                 ttGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGAc  >  7:301343/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGCCGAAGATCATCAGGTGCGTCATCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGAC  >  NC_002947/4382735‑4383028

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