Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,382,821 G→A 11.7% 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,8210GA11.7% 177.7 / 8.7 68intergenic (‑32/+230)PP_3857/PP_3858hypothetical protein/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (4/4);  total (35/33)
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

GTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCG  >  NC_002947/4382677‑4382960
                                                                                                                                                  |                                                                                                                                           
gTATCGACACCGCCGTTTACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGCCATAAAAAAACccgccgatgc                                                                                                                                         >  8:194893/1‑145 (MQ=255)
gTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGc                                                                                                                                         >  4:285730/1‑149 (MQ=255)
       cacCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGtttt                                                                                                                                  <  5:280096/149‑2 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 <  5:321977/149‑1 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 <  4:47190/149‑1 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 <  7:268939/149‑1 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 >  6:321977/1‑149 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 <  3:285730/149‑1 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttt                                                                                                                                 <  6:190574/149‑3 (MQ=255)
          cgccgTTGACATAGACATGCTGCATCTCCGCCGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGCCAATTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 <  1:282151/149‑1 (MQ=255)
          cgccgTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGt                                                                                                                               >  4:109762/1‑149 (MQ=255)
          cgccgTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATACGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttttg                                                                                                                               <  7:47372/149‑5 (MQ=255)
                   cATAGACATGCTGCAGCTCCGACGCGCTGTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATTGGC‑‑TTCCAAAGCTGATATATCAGTCATTCACTTTTCTCCAGGCATCTTAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAt                                                                                                                       <  1:276910/148‑1 (MQ=255)
                   cATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAtt                                                                                                                      >  6:200816/1‑149 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGt                                                                                                                               <  8:224278/139‑1 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGt                                                                                                                               >  7:224278/1‑139 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                     <  7:13531/149‑1 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                     <  5:200816/149‑1 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                     >  4:63714/1‑149 (MQ=255)
                      aGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTgg                                                                                                                   >  4:8685/1‑148 (MQ=255)
                                   cTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAc                                                                                                      >  8:304319/1‑149 (MQ=255)
                                     ccGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttttg                                                                                                                               <  3:51959/122‑5 (MQ=255)
                                     ccGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttttg                                                                                                                               >  4:51959/1‑118 (MQ=255)
                                      cGACGCGCTCTTGGATTGCATTGCCGCTTCGTACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGTGTCAAACGATTTGAGCAGGCTCGTTAcgc                                                                                                    <  7:194893/148‑1 (MQ=255)
                                            gcTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGc                                                                                             >  5:608/1‑149 (MQ=255)
                                              tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAACGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGGGGCGGGTTTGGGGCAACCGATTTGCGCAGGCTCGTTACGCGTCACgat                                                                                            >  4:36941/1‑146 (MQ=255)
                                               ctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGTATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                   >  4:133035/1‑108 (MQ=255)
                                               ctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGTATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                   <  3:133035/108‑1 (MQ=255)
                                                ttGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                   <  8:15828/107‑1 (MQ=255)
                                                ttGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                   >  7:15828/1‑107 (MQ=255)
                                                  ggATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTa                                                                                         <  6:608/147‑1 (MQ=255)
                                                   gATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 <  6:236349/106‑1 (MQ=255)
                                                   gATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                 >  5:236349/1‑106 (MQ=255)
                                                       gcattgcCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                   <  8:333036/148‑1 (MQ=255)
                                                         attgcCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                           >  5:277994/1‑106 (MQ=255)
                                                         attgcCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                           <  6:277994/106‑1 (MQ=255)
                                                                   cGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAg                                                                       >  8:43644/1‑148 (MQ=255)
                                                                        aacaGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGcatcat                                                                 >  6:95596/1‑149 (MQ=255)
                                                                          caGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATcg                                                               >  8:5792/1‑149 (MQ=255)
                                                                           aGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATcgg                                                              >  2:349269/1‑149 (MQ=255)
                                                                             gCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATcggcg                                                            >  1:68861/1‑149 (MQ=255)
                                                                                    tAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCt                                                     <  1:349269/149‑1 (MQ=255)
                                                                                       gC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAg                                                       >  5:31034/1‑144 (MQ=255)
                                                                                       gC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAg                                                       <  6:31034/144‑1 (MQ=255)
                                                                                           ttCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGtttttt                                                                                                                                >  3:109462/1‑66 (MQ=255)
                                                                                           ttCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGtttttt                                                                                                                                <  4:109462/69‑4 (MQ=255)
                                                                                              cAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAcc                                              <  7:304319/148‑1 (MQ=255)
                                                                                              cAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCa                                             <  3:109762/149‑1 (MQ=255)
                                                                                                aaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGa                                                                              <  1:47885/114‑1 (MQ=255)
                                                                                                aaGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGa                                                                              >  2:47885/1‑114 (MQ=255)
                                                                                                      aaaTATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCa                                      >  4:1132/1‑148 (MQ=255)
                                                                                                      aaaTATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCa                                      >  1:46961/1‑148 (MQ=255)
                                                                                                             aGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTca                              >  5:351092/1‑149 (MQ=255)
                                                                                                             aGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTca                              <  7:5792/149‑1 (MQ=255)
                                                                                                                     aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTg                                                    <  4:218196/119‑1 (MQ=255)
                                                                                                                     aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTg                                                    >  3:218196/1‑119 (MQ=255)
                                                                                                                     aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCa                       <  3:1132/148‑1 (MQ=255)
                                                                                                                       ttttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccatc                        <  6:245286/145‑1 (MQ=255)
                                                                                                                       ttttCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatccatc                        >  5:245286/1‑145 (MQ=255)
                                                                                                                          tctcCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCt                  <  2:46961/148‑1 (MQ=255)
                                                                                                                                ggCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGAt                                                                                  >  6:30670/1‑78 (MQ=255)
                                                                                                                                  cATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc         <  2:68861/149‑1 (MQ=255)
                                                                                                                                     cAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCa      >  8:45851/1‑149 (MQ=255)
                                                                                                                                      cccaaaCCCGCCACGGCGGGTTTTGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAt     >  1:325129/4‑149 (MQ=255)
                                                                                                                                       ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt    >  1:143518/3‑149 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt    >  4:281180/1‑149 (MQ=255)
                                                                                                                                        aacaCCCGCCGTGGCTGGTTTTTGTCAAACGATTTGAGCAGGCTCTTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc   <  3:36941/146‑1 (MQ=255)
                                                                                                                                        aaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc   >  7:204989/1‑149 (MQ=255)
                                                                                                                                         aaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCg  >  6:138489/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGC‑‑TTCCAAAGCTGAAATATCAGTCATTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCG  >  NC_002947/4382677‑4382960

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: