Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,382,822 T→C 8.5% intergenic (‑33/+229) 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,8220TC 8.5% 189.6 / 7.0 70intergenic (‑33/+229)PP_3857/PP_3858hypothetical protein/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base C (4/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

TCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGC  >  NC_002947/4382680‑4382968
                                                                                                                                               |                                                                                                                                                  
tCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCggg                                                                                                                                              >  4:263719/1‑149 (MQ=255)
        gccgTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTc                                                                                                                                      <  7:246420/149‑1 (MQ=255)
        gccgTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGAGGGTTTGGGt                                                                                                                                       >  2:261222/1‑148 (MQ=255)
            ttGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAAc                                                                                                                                  <  5:376269/149‑1 (MQ=255)
                 aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                             <  1:48603/149‑1 (MQ=255)
                    gACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGa                                                                                                                           <  2:254268/148‑1 (MQ=255)
                             cAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCg                                                                                                                  <  1:349617/148‑1 (MQ=255)
                               gCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTa                                                                                                               >  2:71859/1‑149 (MQ=255)
                                cTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTCAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTATACAGGCATCAAAAAACCCGCAGTGGCGGGTTTGGGACAAACGATTTGAGCAGGCTCGTTAc                                                                                                              >  8:34216/1‑149 (MQ=255)
                                   cGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcgc                                                                                                            <  2:242575/148‑1 (MQ=255)
                                   cGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCAGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcgcg                                                                                                           >  4:55644/1‑149 (MQ=255)
                                     aCGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGt                                                                                                          >  4:103639/1‑148 (MQ=255)
                                        cgcTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACg                                                                                                      <  3:179824/149‑1 (MQ=255)
                                        cgcTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACg                                                                                                      <  3:29606/149‑1 (MQ=255)
                                           tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTc                                                                                                                   >  7:63603/1‑133 (MQ=255)
                                           tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTc                                                                                                                   <  8:63603/133‑1 (MQ=255)
                                           tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGctc                                                                                                   >  8:179469/1‑149 (MQ=255)
                                             ttGGATTGCATTGCCGCTTCGGACAACAGGCCTGAGATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGctc                                                                                                   >  1:55420/1‑148 (MQ=255)
                                                gATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACgg                                                                                              >  6:303938/1‑149 (MQ=255)
                                                gATTGAATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTGACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACgg                                                                                              >  7:235576/1‑149 (MQ=255)
                                                    gcattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGAt                                                                                          >  6:17526/1‑149 (MQ=255)
                                                     cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCACACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGAtt                                                                                         <  1:71859/149‑1 (MQ=255)
                                                      attgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTg                                                                                        >  7:338651/1‑149 (MQ=255)
                                                      attgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTg                                                                                        <  7:57630/149‑1 (MQ=255)
                                                                   acaacaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGcatc                                                                           <  5:17526/149‑1 (MQ=255)
                                                                          gCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATcggcg                                                                    >  2:317753/1‑149 (MQ=255)
                                                                            cTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCggcggc                                                                  >  5:297224/1‑149 (MQ=255)
                                                                                         ccAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACTCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                       <  8:235576/148‑1 (MQ=255)
                                                                                         ccAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                       <  3:263719/148‑1 (MQ=255)
                                                                                          cAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCa                                                     >  4:229164/1‑149 (MQ=255)
                                                                                            aaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCAtt                                                   <  1:195516/149‑1 (MQ=255)
                                                                                              gCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTAc                                                 <  7:179469/149‑1 (MQ=255)
                                                                                                tGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACg                                                >  4:205137/1‑148 (MQ=255)
                                                                                                tGAAATATCAGTCATTCACTTTTCTCCAGGAATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATATGAGCAGGCTCGTTAAGCGTCACGCACTACGGCGATTGGATCGCCAGCATCAACGGCGTCGAGCTGATTACAATTACg                                                >  5:52979/1‑148 (MQ=255)
                                                                                                       tCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGt                                        >  2:79966/1‑149 (MQ=255)
                                                                                                          gTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCTTCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcat                                     >  7:279361/1‑149 (MQ=255)
                                                                                                          gTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcat                                     <  8:338651/149‑1 (MQ=255)
                                                                                                           tCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatc                                    <  2:55420/149‑1 (MQ=255)
                                                                                                             aTTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTa                                                        >  4:29845/1‑127 (MQ=255)
                                                                                                             aTTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTa                                                        <  3:29845/127‑1 (MQ=255)
                                                                                                             aTTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatcca                                  >  4:196661/1‑149 (MQ=255)
                                                                                                             aTTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatcca                                  >  6:19302/1‑149 (MQ=255)
                                                                                                             aTTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatcca                                  >  3:98451/1‑149 (MQ=255)
                                                                                                                 aattttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGCTCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATGaa                              >  4:229355/3‑149 (MQ=255)
                                                                                                                  cTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAAt                             >  8:135330/1‑149 (MQ=255)
                                                                                                                   ttttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAAt                             >  4:288991/1‑148 (MQ=255)
                                                                                                                     ttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTa                                                  <  2:361236/125‑1 (MQ=255)
                                                                                                                     ttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTa                                                  >  1:361236/1‑125 (MQ=255)
                                                                                                                       ctcCAGGCATCAAAAAACCCTCCGTGGCGGTTTTGGGTCAAACTATTTGAGCAGTCTCGTTACTCGTCACGCTCTACTGCGATTGGATCGCCATCATCATGGTCGGCTAGCTGTTTTCCATTACGCAAAACGTCATCCATCAATGCtt                         <  7:34216/148‑1 (MQ=255)
                                                                                                                        tcCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCa                       >  5:155484/1‑149 (MQ=255)
                                                                                                                          cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                        <  3:103639/146‑1 (MQ=255)
                                                                                                                          cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                        <  3:288991/146‑1 (MQ=255)
                                                                                                                            ggCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGAc                   <  6:297224/149‑1 (MQ=255)
                                                                                                                             gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGttt                                                         >  4:131477/1‑110 (MQ=255)
                                                                                                                             gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCAGCGGCGGCGAGCTGttt                                                         <  3:131477/110‑1 (MQ=255)
                                                                                                                               aTCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCt                <  1:317753/149‑1 (MQ=255)
                                                                                                                                  cccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc                 <  7:303791/142‑1 (MQ=255)
                                                                                                                                  cccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGAGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc                 >  8:303791/4‑145 (MQ=255)
                                                                                                                                  aaaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                        >  2:329562/1‑138 (MQ=255)
                                                                                                                                  aaaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                        <  1:329562/138‑1 (MQ=255)
                                                                                                                                    caaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGCGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc           >  5:74181/2‑149 (MQ=255)
                                                                                                                                    caaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc           >  5:267122/2‑149 (MQ=255)
                                                                                                                                    aaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc           <  6:155484/149‑1 (MQ=255)
                                                                                                                                     aaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCtt                         <  6:180424/134‑1 (MQ=255)
                                                                                                                                     aaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCtt                         >  5:180424/1‑134 (MQ=255)
                                                                                                                                      aaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGt         <  1:79966/149‑1 (MQ=255)
                                                                                                                                         ccgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAg      >  2:231028/1‑149 (MQ=255)
                                                                                                                                          cgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAt             >  7:83270/1‑141 (MQ=255)
                                                                                                                                          cgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAt             <  8:83270/141‑1 (MQ=255)
                                                                                                                                             cgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGc  <  3:55644/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGCGGC  >  NC_002947/4382680‑4382968

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: