Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,058,260 Δ1 bp 5.8% intergenic (‑186/+34) 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,2600A. 5.8% 218.7 / 6.9 69intergenic (‑186/+34)PP_0914/sodBGGDEF domain‑containing protein/superoxide dismutase
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base . (2/2);  total (35/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 = 4.22e-01

CCATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGC  >  NC_002947/1058113‑1058407
                                                                                                                                                   |                                                                                                                                                    
ccATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGca                                                                                                                                                      >  2:1183079/1‑148 (MQ=255)
ccATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGca                                                                                                                                                      <  1:1183079/148‑1 (MQ=255)
 cATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACAt                                                                                                                                                   >  1:1008863/1‑150 (MQ=255)
      gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCCGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGCTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                             <  1:1151227/150‑1 (MQ=255)
      gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                             <  2:749695/150‑1 (MQ=255)
      gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                             >  2:287799/1‑150 (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                                                                                            >  1:721670/1‑142 (MQ=255)
                                                               gTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGGCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTc                                                                                            <  2:721670/142‑1 (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:479225/1‑150 (MQ=255)
                                                                          gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                         >  2:982772/1‑150 (MQ=255)
                                                                          gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                         >  1:686817/1‑150 (MQ=255)
                                                                          gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                         >  1:781693/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                                      <  2:686817/150‑1 (MQ=255)
                                                                                                             aGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTa                                      <  1:491882/150‑1 (MQ=255)
                                                                                                              gCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                               <  1:399488/92‑1 (MQ=255)
                                                                                                              gCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                               >  2:399488/1‑92 (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                                                 <  2:908413/117‑1 (MQ=255)
                                                                                                                                   cacaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAAc                                                 >  1:908413/1‑117 (MQ=255)
                                                                                                                                    acaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGt               <  1:1157545/150‑1 (MQ=255)
                                                                                                                                      aaaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg             >  1:731163/1‑150 (MQ=255)
                                                                                                                                       gagaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                  >  2:811983/4‑64 (MQ=21)
                                                                                                                                       gagaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                  <  1:811983/61‑1 (MQ=18)
                                                                                                                                       aaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGa            >  2:646099/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt           <  2:651496/150‑1 (MQ=255)
                                                                                                                                            aCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGGTGCGGTAGTCGAAgtac       >  2:130536/1‑149 (MQ=255)
                                                                                                                                                ggCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc                                                      <  2:800595/99‑1 (MQ=255)
                                                                                                                                                ggCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc                                                      >  1:800595/1‑99 (MQ=255)
                                                                                                                                                 gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                   >  2:545798/1‑133 (MQ=255)
                                                                                                                                                 gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                   <  1:545798/133‑1 (MQ=255)
                                                                                                                                                 gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGc  >  2:401981/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGC  >  NC_002947/1058113‑1058407

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: