Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,382,821 G→A 9.4% intergenic (‑32/+230) 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,8210GA 9.4% 168.4 / 5.3 63intergenic (‑32/+230)PP_3857/PP_3858hypothetical protein/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base A (4/2);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

GTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAG  >  NC_002947/4382677‑4382964
                                                                                                                                                |                                                                                                                                                 
gTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGc                                                                                                                                               >  2:217637/1‑149 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                       <  6:160710/149‑1 (MQ=255)
          cgccgTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGt                                                                                                                                     >  4:99295/1‑149 (MQ=255)
             cgTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttttgatg                                                                                                                                  <  8:249894/149‑8 (MQ=255)
               ttGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAAc                                                                                                                                <  2:310071/149‑1 (MQ=255)
                 gACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGa                                                                                                                              >  7:199148/1‑149 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                           >  6:173191/1‑149 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTCGGGTCAAACGAttt                                                                                                                           <  7:319549/149‑1 (MQ=255)
                      aGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                         >  8:304004/1‑133 (MQ=255)
                      aGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                         <  7:304004/133‑1 (MQ=255)
                       gACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGa                                                                                                                         <  2:372440/148‑1 (MQ=255)
                       gACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGa                                                                                                                         <  5:187718/148‑1 (MQ=255)
                        aCATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAg                                                                                                                        >  5:53159/1‑148 (MQ=255)
                         cATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAg                                                                                                                        >  3:278085/1‑147 (MQ=255)
                               gcAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCg                                                                                                                >  2:120394/1‑149 (MQ=255)
                                cAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTc                                                                                                                 <  3:31559/147‑1 (MQ=255)
                                  gCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACGCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTa                                                                                                             <  3:923/149‑1 (MQ=255)
                                             ctctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGct                                                                                                 >  1:192/1‑149 (MQ=255)
                                              tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGctc                                                                                                >  2:103682/1‑149 (MQ=255)
                                              tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttt                                                                                                                                       <  6:201015/111‑3 (MQ=255)
                                              tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttt                                                                                                                                       >  5:201015/1‑109 (MQ=255)
                                              tattcgATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCAAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACgat                                                                                                 >  6:289189/6‑146 (MQ=255)
                                                 tGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGACTCGTTACGCGT‑CACGCTCTAc                                                                                             >  8:235295/1‑149 (MQ=255)
                                                   gATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑gg                                                                                          >  7:182127/1‑149 (MQ=255)
                                                    attgcattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTAATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGAGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGc                                                                                         >  5:164846/1‑149 (MQ=255)
                                                        cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGAtt                                                                                     <  8:182127/149‑1 (MQ=255)
                                                        cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGAtt                                                                                     <  4:45066/149‑1 (MQ=255)
                                                        cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGAtt                                                                                     <  8:199148/149‑1 (MQ=255)
                                                        cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGAtt                                                                                     <  8:27841/149‑1 (MQ=255)
                                                        cattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGACCACCCTCCTTACGCGT‑CACGCTCTAC‑GGCGAtt                                                                                     >  2:278071/1‑149 (MQ=255)
                                                          ttgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTgg                                                                                   <  3:99295/149‑1 (MQ=255)
                                                            gcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGAt                                                                                 >  7:86634/1‑149 (MQ=255)
                                                                cTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGcc                                                                             <  4:341678/149‑1 (MQ=255)
                                                                cTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGc                                                                              >  3:154152/1‑148 (MQ=255)
                                                                cTTCGGACAACAGGCCTGAATAGGCTACCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCACGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGc                                                                              >  8:31775/1‑148 (MQ=255)
                                                                   cGGACAACAGGCCTGAATAGGCTTCCAACGCTGAACAATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCa                                                                            <  8:340495/147‑1 (MQ=255)
                                                                     gACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGca                                                                         >  5:90934/1‑148 (MQ=255)
                                                                     gACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGca                                                                         >  8:49835/1‑148 (MQ=255)
                                                                        aacaGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGcatcat                                                                     >  3:179887/1‑149 (MQ=255)
                                                                                      ggCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGt                                                       <  5:173191/149‑1 (MQ=255)
                                                                                       gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGtt                                                      >  7:344657/1‑149 (MQ=255)
                                                                                       gCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGtt                                                      >  4:14148/1‑149 (MQ=255)
                                                                                           ccAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAATACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CTCGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                   <  3:208806/148‑1 (MQ=255)
                                                                                             aaaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCAt                                                >  5:103137/1‑149 (MQ=255)
                                                                                                  tGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCACCAGCATCATCGGCGGCGAGCTGTTTACCATTACg                                            >  4:128243/1‑148 (MQ=255)
                                                                                                           aGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGTGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTca                                  >  6:15892/1‑149 (MQ=255)
                                                                                                           aGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTc                                   <  8:86634/148‑1 (MQ=255)
                                                                                                              cATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTACGGGAGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcat                                 <  1:103682/148‑1 (MQ=255)
                                                                                                                   aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCa                           <  6:90934/148‑1 (MQ=255)
                                                                                                                     ttttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCCTCCATCAATg                        >  8:208135/1‑149 (MQ=255)
                                                                                                                          tcCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCa                   >  6:148344/1‑149 (MQ=255)
                                                                                                                            cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                    <  7:49835/146‑1 (MQ=255)
                                                                                                                               gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACg              >  8:224904/1‑149 (MQ=255)
                                                                                                                               gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACg              >  5:97620/1‑149 (MQ=255)
                                                                                                                               gCATCAAAAAACCAGCAGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGATCGTTACGCGG‑CACGCTCTAC‑GGCGATTGGATCGACAGAATCATCGGCGGCGAGCTGTTTAACATTACGCAAAACGTCATGCATCAATGCTTCATGCCg              >  4:270612/1‑149 (MQ=255)
                                                                                                                                cATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc             <  2:53152/149‑1 (MQ=255)
                                                                                                                                 aTCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCt            >  1:284795/1‑149 (MQ=255)
                                                                                                                                     ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTGCACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTG‑TTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt        >  1:293186/3‑149 (MQ=255)
                                                                                                                                     ccaaaCCCGCCACGGCGGGGTTGGGTCAAACGATTTGCGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCAGCATCGGCGGCCAGCTGTGTAGCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAt         >  1:308559/3‑148 (MQ=255)
                                                                                                                                      caaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc       >  4:58842/2‑149 (MQ=255)
                                                                                                                                      aaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc       <  6:103137/149‑1 (MQ=255)
                                                                                                                                          cccGCCGTGGCTGGTTTGGGTCACACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATGGCCAGCATCATCGGCGGCGAGCTGTTTACCTTTACGCAAAACGTCATCCATCTATGCTTCATGACGCTCATTCGTGa   <  6:164846/149‑1 (MQ=255)
                                                                                                                                          cccGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGa   <  3:14148/149‑1 (MQ=255)
                                                                                                                                           ccgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAg  >  5:122191/1‑149 (MQ=255)
                                                                                                                                           ccgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAg  <  4:179887/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGT‑CACGCTCTAC‑GGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAG  >  NC_002947/4382677‑4382964

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