Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,058,264:1 +T 7.7% 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 7.7% 849.3 / 65.0 208intergenic (‑190/+30)PP_0914/sodBGGDEF domain‑containing protein/superoxide dismutase
Reads supporting (aligned to +/- strand):  ref base . (91/101);  new base T (9/7);  total (100/108)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

ATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCGTGTT  >  NC_002947/1058115‑1058412
                                                                                                                                                      |                                                                                                                                                    
    gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                                  >  1:2952461/1‑150 (MQ=255)
    gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                                  >  1:221502/1‑150 (MQ=255)
    gCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTGTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCg                                                                                                                                                  <  1:4352738/150‑1 (MQ=255)
     cAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCgg                                                                                                                                                 >  2:4407050/1‑150 (MQ=255)
       gtgtATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGt                                                                                                                                               >  2:2225282/1‑150 (MQ=255)
       gtgtATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGt                                                                                                                                               <  2:1309099/150‑1 (MQ=255)
          tATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtc                                                                                                                                            >  1:784113/1‑150 (MQ=255)
            tCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtctc                                                                                                                                          <  1:1757651/150‑1 (MQ=255)
             cGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCCATGTGCCGGGttttttt                                                                                                                                         <  2:365950/150‑5 (MQ=255)
               gTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGtt                                                                                                                                       <  2:3816561/150‑1 (MQ=255)
                ttGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTc                                                                                                                                      <  1:397995/150‑1 (MQ=255)
                 tGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCa                                                                                                                                     <  2:2205122/150‑1 (MQ=255)
                  gCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCaa                                                                                                                                    <  1:3161507/150‑1 (MQ=255)
                   ccTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAg                                                                                                                                   <  1:3100623/150‑1 (MQ=255)
                   ccTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAg                                                                                                                                   <  2:3101398/150‑1 (MQ=255)
                    cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAgg                                                                                                                                  >  2:2506569/1‑150 (MQ=255)
                    cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAgg                                                                                                                                  >  1:1419569/1‑150 (MQ=255)
                    cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAgg                                                                                                                                  >  2:580398/1‑150 (MQ=255)
                    cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAgg                                                                                                                                  >  2:957305/1‑150 (MQ=255)
                     tGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGc                                                                                                                                 >  2:2936093/1‑150 (MQ=255)
                     tGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGc                                                                                                                                 >  2:3334448/1‑150 (MQ=255)
                         cccGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGaa                                                                                                                             <  1:1073783/150‑1 (MQ=255)
                                  ttaattagcTTTATCAATCTTACAGATGTTTGATGTTATCGTGGTATGCTGCCTGATCAAGCAAAACCGTGCCATCCTTTTTTAGGCAACAAAGGCACAAAAAAACACTTCAAATTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAgg                                                                                                                    <  1:283484/143‑1 (MQ=255)
                                   gAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGc                                                                                                                   <  1:1493111/150‑1 (MQ=255)
                                       gcgcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTg                                                                                                               <  1:2623557/150‑1 (MQ=255)
                                        cgcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGa                                                                                                              <  2:1062837/150‑1 (MQ=255)
                                         tcTGTAACAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGaa                                                                                                             <  1:3262070/149‑1 (MQ=255)
                                         gcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGaa                                                                                                             <  1:2293595/150‑1 (MQ=255)
                                         gcTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGaa                                                                                                             <  1:2150245/150‑1 (MQ=255)
                                          cTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACTCAGGCACAAAAAAACCCGGCTCATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAg                                                                                                            <  1:5400784/150‑1 (MQ=255)
                                          cTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAg                                                                                                            <  1:4789038/150‑1 (MQ=255)
                                              aTCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTTGTTCAAGGCCGAATGACTCAGGCCTTGaacgaga                                                                                                        <  2:1389540/150‑6 (MQ=255)
                                              aTCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCt                                                                                                        >  1:341931/1‑150 (MQ=255)
                                               tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                       <  2:3189720/150‑1 (MQ=255)
                                               tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                       >  1:2177017/1‑150 (MQ=255)
                                               tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                       >  1:4415855/1‑150 (MQ=255)
                                               tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                                       <  2:203604/150‑1 (MQ=255)
                                                cAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTg                                                                                                      >  1:4198523/1‑150 (MQ=255)
                                                 aaTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGtt                                                                                                                                              >  1:3824180/1‑109 (MQ=255)
                                                 aaTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGc                                                                                                     >  1:1721086/1‑150 (MQ=255)
                                                  aTCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGcc                                                                                                    >  1:2465029/1‑150 (MQ=255)
                                                   tCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCt                                                                                                   <  2:1972851/150‑1 (MQ=255)
                                                   tCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCt                                                                                                   >  1:3813508/1‑150 (MQ=255)
                                                      tCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCg                                                                                                >  2:236501/1‑150 (MQ=255)
                                                       ccAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGa                                                                                               <  1:2742323/150‑1 (MQ=255)
                                                          gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACt                                                                                            <  1:4181559/150‑1 (MQ=255)
                                                          gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACt                                                                                            >  2:4229080/1‑150 (MQ=255)
                                                              tttGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                        >  2:339203/1‑150 (MQ=255)
                                                              tttGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                        <  1:950227/150‑1 (MQ=255)
                                                                  aTGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCa                                                                                    >  2:1555564/1‑150 (MQ=255)
                                                                  aTGTTAGCGTGGTATGCGGCCAGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCa                                                                                    >  2:1915193/1‑150 (MQ=255)
                                                                         cGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAgg                                                                             <  2:2897992/150‑1 (MQ=255)
                                                                          gTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGc                                                                            >  2:5213277/1‑150 (MQ=255)
                                                                           tGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                                        <  2:714465/137‑1 (MQ=255)
                                                                                tGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGtt                                                                      <  2:686934/150‑1 (MQ=255)
                                                                                   ggCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGAc                                                                   >  1:1241548/1‑150 (MQ=255)
                                                                                      cTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAg                                                                <  1:5234918/150‑1 (MQ=255)
                                                                                          gcaagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTc                                                            <  2:3258638/150‑1 (MQ=255)
                                                                                             agcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAg                                                         <  2:4702461/150‑1 (MQ=255)
                                                                                              gcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGa                                                        >  2:5018451/1‑150 (MQ=255)
                                                                                                 aaaCCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACg                                                     <  2:1472972/150‑1 (MQ=255)
                                                                                                 aaaCCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACg                                                     <  1:2031280/150‑1 (MQ=255)
                                                                                                     cGTGCCAGCCTGTGTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTc                                                 <  2:4142856/150‑1 (MQ=255)
                                                                                                       tGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGa                                               >  1:4945140/1‑150 (MQ=255)
                                                                                                       tGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTATCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGa                                               <  2:864540/150‑1 (MQ=255)
                                                                                                        gCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGAc                                              <  2:3993423/150‑1 (MQ=255)
                                                                                                          cAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAg                                                                                                            <  1:4906779/86‑1 (MQ=255)
                                                                                                          cAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAg                                                                                                            >  2:4906779/1‑86 (MQ=255)
                                                                                                              cTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGACTGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACtt                                        <  2:2355163/150‑1 (MQ=255)
                                                                                                               tGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTc                                       <  1:4828259/150‑1 (MQ=255)
                                                                                                               tGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTc                                       >  2:2585333/1‑150 (MQ=255)
                                                                                                                gTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGATCTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCg                                      >  2:4043755/1‑150 (MQ=255)
                                                                                                                gTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCg                                      <  1:3054917/150‑1 (MQ=255)
                                                                                                                 tttGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCgg                                     >  1:1003440/1‑150 (MQ=255)
                                                                                                                 tttGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCgg                                     >  1:4088454/1‑150 (MQ=255)
                                                                                                                           cacaGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTg                           <  2:1715867/150‑1 (MQ=255)
                                                                                                                            acaGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTc                                       >  2:4518413/1‑137 (MQ=255)
                                                                                                                            acaGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTc                                       <  1:4518413/137‑1 (MQ=255)
                                                                                                                            acaGGCACAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCg                         >  1:15911/1‑150 (MQ=255)
                                                                                                                             caGGCACAAAAAAACCCGGCCCATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCg                         <  2:3542903/150‑1 (MQ=255)
                                                                                                                              aGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCgg                        >  2:835743/1‑150 (MQ=255)
                                                                                                                                gCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCTGGTTGCGGtc                      >  1:4694348/1‑149 (MQ=255)
                                                                                                                                gCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTa                      <  1:580398/150‑1 (MQ=255)
                                                                                                                                gCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGATCAAGGCCGAATGACTCAGGCCATGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCCGTAGACCAGGATCCAGAACGCCTCGAAGTCCTTCGGACGCAGGTTGCggga                      >  2:3765663/1‑148 (MQ=255)
                                                                                                                                   caAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTc                   >  2:430961/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg                  >  1:3594962/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg                  >  2:1841365/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg                  >  1:272371/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg                  >  2:4467586/1‑150 (MQ=255)
                                                                                                                                     nagaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGa                 >  1:1672712/4‑150 (MQ=255)
                                                                                                                                     gagaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACGTCGGACGCAGGTTGCGGTAGTCGa                 >  2:3168061/4‑150 (MQ=255)
                                                                                                                                     gagaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGAACGAATGACTCAGGCCTTGAATGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTGCGGACGCAGGTTGCGGTAGACGa                 >  1:5289013/4‑150 (MQ=255)
                                                                                                                                      agaaaCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt                >  2:2297696/3‑150 (MQ=255)
                                                                                                                                        aaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATgt              >  2:1182835/1‑150 (MQ=255)
                                                                                                                                          aCCCGGCCATGTGCCGGGTTTTTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCCGGTTGCGGTAGTCGATgtcg            >  2:5250658/1‑148 (MQ=255)
                                                                                                                                          aCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACAAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATgtag            >  2:5313425/1‑150 (MQ=255)
                                                                                                                                            ccGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATgtagta          <  1:835743/150‑1 (MQ=255)
                                                                                                                                               gCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGc       >  1:4602959/1‑150 (MQ=255)
                                                                                                                                                cacaTGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCg      >  1:1532849/1‑150 (MQ=255)
                                                                                                                                                    tGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCGTGtt  <  1:4229080/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATTAGCAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATGTAGTAGGCGTGTT  >  NC_002947/1058115‑1058412

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