Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,058,264:1 +T 5.9% intergenic (‑190/+30) PP_0914 ← / ← sodB GGDEF domain‑containing protein/superoxide dismutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,058,2641.T 5.9% 279.0 / 10.2 68intergenic (‑190/+30)PP_0914/sodBGGDEF domain‑containing protein/superoxide dismutase
Reads supporting (aligned to +/- strand):  ref base . (32/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCGTGTTC  >  NC_002947/1058119‑1058413
                                                                                                                                                  |                                                                                                                                                     
gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCCGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGCTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                                   <  1:1151227/150‑1 (MQ=255)
gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                                   >  2:287799/1‑150 (MQ=255)
gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                                   <  2:749695/150‑1 (MQ=255)
 cAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCgg                                                                                                                                                  <  1:599356/150‑1 (MQ=255)
  aGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCggg                                                                                                                                                 >  1:816960/1‑150 (MQ=255)
   gtgtATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTAACCAACTTCCAGAGGTTTGATGTTAGCGTGGTCTACTTCTTTAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCCATGTGCCGGGt                                                                                                                                                <  1:270783/150‑1 (MQ=255)
       aTCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtct                                                                                                                                            >  1:44432/1‑150 (MQ=255)
        tCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTACAGCTACACATTGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtctc                                                                                                                                           <  2:31334/150‑1 (MQ=255)
                        gCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGa                                                                                                                           <  1:948987/150‑1 (MQ=255)
                                     gcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGaa                                                                                                              <  1:556772/150‑1 (MQ=255)
                                       tGTATCAATCTTCCAGAGGTTTGATGTTCGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGAGAGATGGCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAgg                                                                                                            >  1:62579/1‑150 (MQ=255)
                                           tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                        >  2:491882/1‑150 (MQ=255)
                                              aTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACt                                                                                                                         <  1:777611/130‑1 (MQ=255)
                                              aTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACt                                                                                                                         >  2:777611/1‑130 (MQ=255)
                                                     agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                                     <  1:583527/143‑1 (MQ=255)
                                                     agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                                     >  2:583527/1‑143 (MQ=255)
                                                     agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAAc                                                                                              >  1:246455/1‑150 (MQ=255)
                                                        ggttataTGTTAGCGTGGTATGCGGCCTAAGCAAGCAAACGCGTGCCAGCCTGTTTGAGGCAACACACGCACAAAAAAACCCGGCACATGGACGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGc                                                                                           <  2:476790/144‑1 (MQ=255)
                                                         gTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGGCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTc                                                                                                  <  2:721670/142‑1 (MQ=255)
                                                         gTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGGCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTc                                                                                                  >  1:721670/1‑142 (MQ=255)
                                                         gTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACt                                                                                                                         >  1:437198/1‑119 (MQ=255)
                                                         gTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACt                                                                                                                         <  2:437198/119‑1 (MQ=255)
                                                            tGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCgg                                                                                       <  2:44432/150‑1 (MQ=255)
                                                               tGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGc                                                                                                                                  >  2:632438/1‑104 (MQ=255)
                                                               tGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGc                                                                                                                                  <  1:632438/104‑1 (MQ=255)
                                                               tGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCaa                                                                                    >  1:897403/1‑150 (MQ=255)
                                                                    gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                               >  2:982772/1‑150 (MQ=255)
                                                                    gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                               >  1:781693/1‑150 (MQ=255)
                                                                    gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                               >  2:479225/1‑150 (MQ=255)
                                                                    gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                               >  1:686817/1‑150 (MQ=255)
                                                                    gCGTGGTACGCGGCCTGACCACGCAACACTGTGCCAGTCTGTTTGCGGAAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACCGCTCGGCAACGAAg                                                                               <  2:62579/150‑1 (MQ=255)
                                                                        ggTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGccc                                                                           <  1:738276/150‑1 (MQ=255)
                                                                         gTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCa                                                                          >  1:1001099/1‑150 (MQ=255)
                                                                          tATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAg                                                                         <  2:246455/150‑1 (MQ=255)
                                                                            tGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGtt                                                                       <  2:816960/150‑1 (MQ=255)
                                                                                gCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCATGCATGGCCCATGTGAcc                                                                   <  2:897403/150‑1 (MQ=255)
                                                                                       caagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGAACAGGTTcc                                                            >  2:1157545/1‑150 (MQ=255)
                                                                                        aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGGCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa                                                           >  1:764749/1‑150 (MQ=255)
                                                                                          gcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGa                                                         >  1:296186/1‑150 (MQ=255)
                                                                                           caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGaa                                                        <  1:287799/150‑1 (MQ=255)
                                                                                                 cGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTc                                                  <  2:25025/150‑1 (MQ=255)
                                                                                                   tGCCAGCCTGTTTGAGGCAAAACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGa                                                >  2:595953/1‑150 (MQ=255)
                                                                                                   gtgttgCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGa                                                >  2:71862/6‑150 (MQ=255)
                                                                                                    gCCAGCCTGTTTGCGGCAACACAGGCAAAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGAc                                               <  2:1001099/150‑1 (MQ=255)
                                                                                                       aGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTa                                            <  1:491882/150‑1 (MQ=255)
                                                                                                       aGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTa                                            <  2:686817/150‑1 (MQ=255)
                                                                                                        gCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                                     >  2:399488/1‑92 (MQ=255)
                                                                                                        gCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                                     <  1:399488/92‑1 (MQ=255)
                                                                                                             tttGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGCCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCgg                                      >  2:781216/1‑150 (MQ=255)
                                                                                                              ttGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGa                                     <  1:982772/150‑1 (MQ=255)
                                                                                                                    cAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAgg                               >  2:588741/1‑150 (MQ=255)
                                                                                                                            gCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                         <  1:1019807/84‑1 (MQ=255)
                                                                                                                            gCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                         >  2:1019807/1‑84 (MQ=255)
                                                                                                                            gCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTa                       <  1:479225/150‑1 (MQ=255)
                                                                                                                             cacaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAAc                                                       >  1:908413/1‑117 (MQ=255)
                                                                                                                             cacaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAAc                                                       <  2:908413/117‑1 (MQ=255)
                                                                                                                              acaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGt                     <  1:1157545/150‑1 (MQ=255)
                                                                                                                                aaaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg                   >  1:731163/1‑150 (MQ=255)
                                                                                                                                 gagaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                        <  1:811983/61‑1 (MQ=18)
                                                                                                                                 gagaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                        >  2:811983/4‑64 (MQ=21)
                                                                                                                                 aaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGa                  >  2:646099/1‑150 (MQ=255)
                                                                                                                                  aaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt                 <  2:651496/150‑1 (MQ=255)
                                                                                                                                      aCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGGTGCGGTAGTCGAAgtac             >  2:130536/1‑149 (MQ=255)
                                                                                                                                          ggCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc                                                            >  1:800595/1‑99 (MQ=255)
                                                                                                                                          ggCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc                                                            <  2:800595/99‑1 (MQ=255)
                                                                                                                                           gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                         >  2:545798/1‑133 (MQ=255)
                                                                                                                                           gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                         <  1:545798/133‑1 (MQ=255)
                                                                                                                                           gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGc        >  2:401981/1‑150 (MQ=255)
                                                                                                                                                 ggCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCGTGTTc  >  1:355004/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCGTGTTC  >  NC_002947/1058119‑1058413

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