Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,382,821 G→A 15.5% 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,8210GA15.5% 164.4 / 20.4 71intergenic (‑32/+230)PP_3857/PP_3858hypothetical protein/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base A (6/5);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

GTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGC  >  NC_002947/4382677‑4382965
                                                                                                                                                  |                                                                                                                                                 
gTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGc                                                                                                                                               <  4:190801/149‑1 (MQ=255)
  aTCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCAcggcgg                                                                                                                                             >  6:357623/1‑149 (MQ=255)
      acacCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGttt                                                                                                                                         >  4:54134/1‑149 (MQ=255)
       cacCGCCCGTGTCATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCTCTGCCGCGTCGGACAACAGGCCTGAATAGGCTGCCAATGCTGAAATATCAGTCC‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGtttt                                                                                                                                        <  8:132171/149‑2 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                       >  8:137854/1‑149 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                       <  7:102504/149‑1 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttt                                                                                                                                       <  3:221829/149‑3 (MQ=255)
        acCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCACGGCGGGttttt                                                                                                                                       <  5:357623/149‑3 (MQ=255)
              gTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                <  3:54134/149‑1 (MQ=255)
              gTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                <  5:167075/149‑1 (MQ=255)
              gTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                >  7:234646/1‑149 (MQ=255)
              gTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                <  5:287842/149‑1 (MQ=255)
              gTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCAGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCaaa                                                                                                                                >  5:140586/1‑149 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                          <  2:254728/149‑1 (MQ=255)
                    aTAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGAttt                                                                                                                          >  6:180195/1‑149 (MQ=255)
                       gACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGa                                                                                                                        <  5:260777/148‑1 (MQ=255)
                              tgcAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTc                                                                                                                >  5:172231/1‑149 (MQ=255)
                                cAGCTCCGACGCTCGCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCT‑‑TTGTCTTTTCTCCACGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGt                                                                                                              <  6:140586/149‑1 (MQ=255)
                                cAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGt                                                                                                              >  3:212345/1‑149 (MQ=255)
                                 aGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGtt                                                                                                             <  3:326809/149‑1 (MQ=255)
                                    tCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcg                                                                                                          >  6:118094/1‑149 (MQ=255)
                                      cGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcgcg                                                                                                        >  5:227179/1‑149 (MQ=255)
                                      cGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTAcgc                                                                                                         <  5:180195/148‑1 (MQ=255)
                                        aCGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTc                                                                                                      <  7:137854/149‑1 (MQ=255)
                                              tctTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGctc                                                                                                >  3:243655/1‑149 (MQ=255)
                                                  ggATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCAAGCTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGAGTCACGCTCTa                                                                                              >  1:113153/1‑149 (MQ=255)
                                                      tgcattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCg                                                                                         <  4:212345/148‑1 (MQ=255)
                                                       gcattgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGa                                                                                        <  5:221605/148‑1 (MQ=255)
                                                         attgcCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGAtt                                                                                      <  3:295633/148‑1 (MQ=255)
                                                                    ggACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                       >  4:335149/1‑89 (MQ=255)
                                                                    ggACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTgg                                                                                                                                       <  3:335149/89‑1 (MQ=255)
                                                                    ggACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGca                                                                          <  6:227179/149‑1 (MQ=255)
                                                                          caGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCg                                                                                         >  7:259106/1‑128 (MQ=255)
                                                                          caGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCg                                                                                         <  8:259106/128‑1 (MQ=255)
                                                                               cTGATTTGGCTTCCAATTCTGAAATATCAGTCA‑‑TTCACTTTTGTCCAGGCATCAATAATACCGCCGTGGCGGGTTTGTGTCAAACGATTTGAGCAGGCTGGTTACGCGTCACGCTATAAGGCGATTGGATCGCCAGCATAATcgggtgc                                                               <  8:234646/149‑5 (MQ=255)
                                                                                     aGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTg                                                         >  6:99855/1‑149 (MQ=255)
                                                                                       gCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGtt                                                       >  3:187865/1‑149 (MQ=255)
                                                                                         ttCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTa                                                     >  3:48037/1‑149 (MQ=255)
                                                                                           ccAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                    <  4:48037/148‑1 (MQ=255)
                                                                                           ccAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTAc                                                    >  5:35792/1‑148 (MQ=255)
                                                                                           ccAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGCCGAGCTGTTTAc                                                    <  6:172231/148‑1 (MQ=255)
                                                                                                gCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCggcggc                                                               >  4:120144/1‑132 (MQ=255)
                                                                                                gCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCggcggc                                                               <  3:120144/132‑1 (MQ=255)
                                                                                                 cTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACg                                             <  6:35792/149‑1 (MQ=255)
                                                                                                    aaaTATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGc                                            <  2:113153/147‑1 (MQ=255)
                                                                                                     aaTATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAGGACGATTTGAGCAGGCTAGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATAGGCGGCGAGCTGTTTAACATTACGCa                                           >  1:108636/1‑148 (MQ=255)
                                                                                                           aGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTca                                   >  1:331536/1‑149 (MQ=255)
                                                                                                                  ttCACTTTTCTCCAGGCATCAAAAAACCCGCCATGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatcca                               >  1:6713/1‑148 (MQ=255)
                                                                                                                     aCTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCaa                           >  8:32452/1‑149 (MQ=255)
                                                                                                                      cTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAAt                          <  2:331536/149‑1 (MQ=255)
                                                                                                                              cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatc                                 >  2:141401/1‑134 (MQ=255)
                                                                                                                              cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTcatc                                 <  1:141401/134‑1 (MQ=255)
                                                                                                                              cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTc                     <  4:243655/146‑1 (MQ=255)
                                                                                                                              cAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCAt                   >  2:118211/1‑148 (MQ=255)
                                                                                                                               aGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGa                 >  2:366969/1‑149 (MQ=255)
                                                                                                                                 gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCatcatc                                                                     <  6:10502/95‑1 (MQ=255)
                                                                                                                                 gCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCatcatc                                                                     >  5:10502/1‑95 (MQ=255)
                                                                                                                                  cATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGc              <  8:73157/149‑1 (MQ=255)
                                                                                                                                    tCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCa                                                                             <  6:192937/84‑1 (MQ=255)
                                                                                                                                    tCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCa                                                                             >  5:192937/1‑84 (MQ=255)
                                                                                                                                     cAAAAACCCCGCCGTGGCGGGTTTGTGGCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCCTCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCa           <  4:254386/149‑1 (MQ=255)
                                                                                                                                     acccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCtt                      >  6:365295/5‑138 (MQ=255)
                                                                                                                                     acccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCtt                      <  5:365295/134‑1 (MQ=255)
                                                                                                                                     acccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCa           >  4:68705/5‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCa           <  5:99855/148‑1 (MQ=255)
                                                                                                                                       ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt         >  7:111721/3‑149 (MQ=255)
                                                                                                                                       ccaaaCCCGCCACGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt         >  4:318115/3‑149 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCAtt         >  6:380334/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc        >  6:22453/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCGCCATGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTc        <  2:6713/149‑1 (MQ=255)
                                                                                                                                            cccGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGa    <  1:118211/149‑1 (MQ=255)
                                                                                                                                              cgccgTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGc  <  4:187865/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTATCGACACCGCCGTTGACATAGACATGCTGCAGCTCCGACGCGCTCTTGGATTGCATTGCCGCTTCGGACAACAGGCCTGAATAGGCTTCCAAAGCTGAAATATCAGTCA‑‑TTCACTTTTCTCCAGGCATCAAAAAACCCGCCGTGGCGGGTTTGGGTCAAACGATTTGAGCAGGCTCGTTACGCGTCACGCTCTACGGCGATTGGATCGCCAGCATCATCGGCGGCGAGCTGTTTACCATTACGCAAAACGTCATCCATCAATGCTTCATGACGCTCATTCGTGAGC  >  NC_002947/4382677‑4382965

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: