Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,058,260 Δ1 bp 7.5% 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. 7.5% 276.3 / 17.6 93intergenic (‑186/+34)PP_0914/sodBGGDEF domain‑containing protein/superoxide dismutase
Reads supporting (aligned to +/- strand):  ref base A (39/45);  new base . (4/3);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

GACCATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGC  >  NC_002947/1058111‑1058407
                                                                                                                                                     |                                                                                                                                                      
gACCATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGca                                                                                                                                                        >  1:192059/1‑150 (MQ=255)
    aTTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGAAACACAGGCACAAAAAAACCCGGCACATg                                                                                                                                                    <  2:580337/150‑1 (MQ=255)
           gtgtATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGt                                                                                                                                            >  2:1073247/1‑150 (MQ=255)
               attggTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGTTCTGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGGTtct                                                                                                                                        >  1:1327070/4‑150 (MQ=255)
                tCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtctc                                                                                                                                       >  2:5570/1‑150 (MQ=255)
                  ggTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGt                                                                                                                                     >  1:289257/1‑150 (MQ=255)
                   gTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCg                                                                                                                                      <  2:345590/148‑1 (MQ=255)
                   gTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCg                                                                                                                                      >  1:345590/1‑148 (MQ=255)
                   gTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGtt                                                                                                                                    <  2:192059/150‑1 (MQ=255)
                      gCCTGACCCCGTAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTTTTTGCTGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCaa                                                                                                                                 <  2:915437/150‑1 (MQ=255)
                      gCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCaa                                                                                                                                 >  1:1815566/1‑150 (MQ=255)
                        cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAgg                                                                                                                               >  2:138182/1‑150 (MQ=255)
                               cGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATg                                                                                                                        >  1:1718928/1‑150 (MQ=255)
                               cGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCATTATTGGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATg                                                                                                                        >  1:1019916/1‑150 (MQ=255)
                                gCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGa                                                                                                                       <  2:320787/150‑1 (MQ=255)
                                    tttGAATGCGCTGTATCAATCTTCCACAGGTTTGATGTTAGCGCGGTGTGCCGGCTGAGCAAGCAAAACCGTGCCAGCCTTTTTGCGGCAACACAGGCACAAAAACACCCGGCCCATGGCCGGGTTTCTCGTTCAAGGCCGAATTACTCa                                                                                                                 <  2:1389996/149‑1 (MQ=255)
                                     ttGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAg                                                                                                                <  1:5570/150‑1 (MQ=255)
                                     ttGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAg                                                                                                                >  2:275370/1‑150 (MQ=255)
                                     ttGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAg                                                                                                                <  2:289257/150‑1 (MQ=255)
                                     ttGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAg                                                                                                                >  2:844228/1‑150 (MQ=255)
                                            cgcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGa                                                                                                         <  1:275370/150‑1 (MQ=255)
                                            cgcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCCATGTGCCGGGttttt                                                                                                                                        <  1:1235652/121‑3 (MQ=255)
                                            cgcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCCATGTGCCGGGttttt                                                                                                                                        >  2:1235652/1‑119 (MQ=255)
                                              cTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAg                                                                                                       <  1:1441381/150‑1 (MQ=255)
                                                 tATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTAGGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTc                                                                                                    <  1:1383722/150‑1 (MQ=255)
                                                   tCAAGCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGACCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATTGTCGCCTTTCTCGTTCAAGCCCGAATGCCTCACACCTTGAAGGTCtt                                                                                                  <  2:691281/150‑1 (MQ=255)
                                                    cAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTg                                                                                                 >  1:1264827/1‑150 (MQ=255)
                                                     aaTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCGTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGc                                                                                                <  2:818375/150‑1 (MQ=255)
                                                     aaTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGc                                                                                                >  1:1472418/1‑150 (MQ=255)
                                                      aTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGGTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                               >  1:1869029/1‑150 (MQ=255)
                                                       tCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCt                                                                                              >  2:1731793/1‑150 (MQ=255)
                                                            cAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGaa                                                                                         >  2:60763/1‑150 (MQ=255)
                                                             agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAAc                                                                                        <  2:1472418/150‑1 (MQ=255)
                                                              gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACt                                                                                       >  1:1083026/1‑150 (MQ=255)
                                                                  tttGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                   >  1:1636156/1‑150 (MQ=255)
                                                                   ttGATGTTAGCGTAACAATCGGCCTGATCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCg                                                                                  <  2:1869029/150‑1 (MQ=255)
                                                                       tGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCaa                                                                              <  2:1083026/150‑1 (MQ=255)
                                                                         ttAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGaa                                                                                                        <  2:1086488/122‑1 (MQ=255)
                                                                         ttAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGACGAATGACTCAGGCCTTGaa                                                                                                        >  1:1086488/1‑122 (MQ=255)
                                                                            gCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAg                                                                         <  2:529117/150‑1 (MQ=255)
                                                                              gTGGTATGCGTCCTGAGCAAGCAAACCCGTGCCAGCCTGTTTGCTGGAACACAGGCACAAAAAAACCCGGCGCATGACCGTGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGc                                                                       <  1:620211/150‑1 (MQ=255)
                                                                               tGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGcc                                                                      <  1:60763/150‑1 (MQ=255)
                                                                                      cGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGTCTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGa                                                               <  2:1420757/150‑1 (MQ=255)
                                                                                        gCCTGAGCAAGCAAAACCGTGCCAGCCTGTTCGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGAcc                                                             >  2:1194782/1‑150 (MQ=255)
                                                                                          cTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCAGCAACACCGTTAGCAGCAGCTTCGGGACATGGCCGGGTTACTCGTTCAAGGCAGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAg                                                           >  1:1493967/1‑150 (MQ=255)
                                                                                            gAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACAGCCTCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGt                                                         <  1:991123/150‑1 (MQ=255)
                                                                                             agcaagcaAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGtt                                                        >  1:630917/1‑150 (MQ=255)
                                                                                              gcaagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTc                                                       <  1:844228/150‑1 (MQ=255)
                                                                                               caagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCATGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc                                                      <  2:1718928/150‑1 (MQ=255)
                                                                                                aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGaa                                                                          <  1:1431501/129‑1 (MQ=255)
                                                                                                aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGaa                                                                          >  2:1431501/1‑129 (MQ=255)
                                                                                                aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa                                                     >  1:368175/1‑150 (MQ=255)
                                                                                                aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGCCGCTGATATCGAGGGCAAGAAGGGCTTGCCTTCGAACTCCGCAGCTACGAACCGCTAGGTGACCAGGTTCCa                                                     >  1:35102/1‑150 (MQ=255)
                                                                                                aagcaaAACCGTGCAAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa                                                     >  2:477676/1‑150 (MQ=255)
                                                                                                  gcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                  <  2:1612981/103‑1 (MQ=255)
                                                                                                  gcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                  >  1:1612981/1‑103 (MQ=255)
                                                                                                   caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGaa                                                                                         >  1:868887/1‑111 (MQ=255)
                                                                                                   caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGaa                                                                                         <  2:868887/111‑1 (MQ=255)
                                                                                                   caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGaa                                                                          >  2:61886/1‑126 (MQ=255)
                                                                                                   caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGaa                                                                          <  1:61886/126‑1 (MQ=255)
                                                                                                   caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGaa                                                  <  1:1811596/150‑1 (MQ=255)
                                                                                                   caaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGaa                                                  >  2:1811596/1‑150 (MQ=255)
                                                                                                        cgatctcAGCCTGTTTGCGGCAACACAGGCACCAAAGGACCCTGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCt                                             <  2:1610270/144‑1 (MQ=255)
                                                                                                           tGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAAt                                                                                                                         >  1:163500/1‑73 (MQ=255)
                                                                                                           tGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAAt                                                                                                                         <  2:163500/73‑1 (MQ=255)
                                                                                                             ccAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACg                                        >  2:848643/1‑150 (MQ=255)
                                                                                                              cAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                   >  1:1720909/1‑106 (MQ=255)
                                                                                                              cAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                   <  2:1720909/106‑1 (MQ=255)
                                                                                                              cAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGt                                       >  2:70842/1‑150 (MQ=255)
                                                                                                              cAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGt                                       <  1:848643/150‑1 (MQ=255)
                                                                                                               aGCCTGTTTGCGGCAACACAGGCACAAAAAAACCTGACACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAAGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTcagatcg                                      >  1:1610270/1‑144 (MQ=255)
                                                                                                               aGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTa                                      >  2:335017/1‑150 (MQ=255)
                                                                                                               aGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTGCTCGTTCAAGGCCGAATGACTAAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTa                                      <  1:1099092/150‑1 (MQ=255)
                                                                                                                   tGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa                                                     >  2:158572/1‑131 (MQ=255)
                                                                                                                   tGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa                                                     <  1:158572/131‑1 (MQ=255)
                                                                                                                          ggCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTGTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGcg                           <  2:1327070/150‑2 (MQ=255)
                                                                                                                                acaGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGc                     <  2:1636156/150‑1 (MQ=255)
                                                                                                                                 caGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCg                    <  1:1194782/150‑1 (MQ=255)
                                                                                                                                  aGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                   <  2:1264827/150‑1 (MQ=255)
                                                                                                                                  aGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                   <  1:70842/150‑1 (MQ=255)
                                                                                                                                       acgagaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGAc                                                              >  1:177699/6‑102 (MQ=38)
                                                                                                                                       acgagaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGAc                                                              <  2:177699/97‑1 (MQ=38)
                                                                                                                                         aaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGa            <  2:630917/150‑1 (MQ=255)
                                                                                                                                          agaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt           >  1:1265057/3‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt           >  1:777104/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt           <  1:477676/150‑1 (MQ=255)
                                                                                                                                            aaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCAAGTTGACCAGGTTCCAGAACGCCTc                                            >  2:971428/1‑115 (MQ=255)
                                                                                                                                            aaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCAAGTTGACCAGGTTCCAGAACGCCTc                                            <  1:971428/115‑1 (MQ=255)
                                                                                                                                             aaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGCTCAAGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATgta        <  1:138182/150‑1 (MQ=255)
                                                                                                                                               cccGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                   <  2:456317/137‑1 (MQ=255)
                                                                                                                                               cccGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                   >  1:456317/1‑137 (MQ=255)
                                                                                                                                                 cGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGAAAAGAGCGCATTACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGtagtag    <  1:1107348/150‑1 (MQ=255)
                                                                                                                                                   gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGc  <  2:777104/150‑1 (MQ=255)
                                                                                                                                                   gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGc  <  2:1815566/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GACCATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGC  >  NC_002947/1058111‑1058407

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