Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,253,437 A→T 6.5% intergenic (‑77/+30) PP_1096 ← / ← metG electron transport complex protein/methionine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,253,4370AT 6.5% 266.6 / 8.4 93intergenic (‑77/+30)PP_1096/metGelectron transport complex protein/methionine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (49/37);  new base T (3/3);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

AGAGTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCGGGG  >  NC_002947/1253288‑1253586
                                                                                                                                                       |                                                                                                                                                     
agagTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCa                                                                                                                                                       <  1:916875/150‑1 (MQ=255)
  agTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGGAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAgg                                                                                                                                                     >  1:3039603/1‑150 (MQ=255)
    tGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTACTTGTGGGAGCGGGCAGGcc                                                                                                                                                   >  2:2984091/1‑150 (MQ=255)
      gCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGc                                                                                                                                                 <  1:2853413/150‑1 (MQ=255)
       cTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCt                                                                                                                                                >  1:1316947/1‑150 (MQ=255)
        tGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGCCTGCCCGCTc                                                                                                                                               >  2:1145944/1‑150 (MQ=255)
         gTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGCCTGCCCGCTcc                                                                                                                                              <  1:1145944/150‑1 (MQ=255)
         gTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGGTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGGCGCTcc                                                                                                                                              >  1:98751/1‑150 (MQ=255)
                 aGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATgg                                                                                                                                      >  2:2993070/1‑150 (MQ=255)
                    cacaGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCCGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGccc                                                                                                                                   <  1:851299/150‑1 (MQ=255)
                     acaGGTGATTGATCAGTGCGGCGATGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCg                                                                                                                                  <  1:2568298/150‑1 (MQ=255)
                       aGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCa                                                                                                                                >  1:1266137/1‑150 (MQ=255)
                       aGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCAGGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGGGCCATTGTGGGAGCGGCCTGCCCGCTCCCACAAttacgca                                                                                                                                  <  2:1501934/150‑7 (MQ=25)
                            cctgatCAGCGCCGCACTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCg                                                                                                                           >  2:2697097/3‑150 (MQ=255)
                                 tCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAggg                                                                                                                             <  1:554983/143‑1 (MQ=255)
                                 tCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAggg                                                                                                                             >  2:554983/1‑143 (MQ=255)
                                   aGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACtt                                                                                                                    >  1:210103/1‑150 (MQ=255)
                                   aGTGCCGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGAATATAGCCGGGAATCGGTGCGATTTTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACtt                                                                                                                    <  1:2743178/150‑1 (MQ=255)
                                    gTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTg                                                                                                                   <  1:405727/150‑1 (MQ=255)
                                     tGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGa                                                                                                                  >  1:1244003/1‑150 (MQ=255)
                                       ccgcacTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATg                                                                                                                <  1:2697097/145‑1 (MQ=255)
                                          cgcTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGCCTGCCCGCTCCCACAAtta                                                                                                                                      <  1:445248/125‑3 (MQ=39)
                                          cgcTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGCCTGCCCGCTCCCACAAtta                                                                                                                                      >  2:445248/1‑123 (MQ=39)
                                           gcTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGct                                                                                                            <  2:2551329/150‑1 (MQ=255)
                                            cTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGGAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctg                                                                                                           <  2:3039603/150‑1 (MQ=255)
                                             tGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctgg                                                                                                          >  2:1027498/1‑150 (MQ=255)
                                               aTCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctggcc                                                                                                        >  1:2877859/1‑150 (MQ=255)
                                                tCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCcac                                                                                                                                           <  2:2088065/114‑1 (MQ=255)
                                                tCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCcac                                                                                                                                           >  1:2088065/1‑114 (MQ=255)
                                                tCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctggcct                                                                                                       >  2:2645426/1‑150 (MQ=255)
                                                    gACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTAc                                                                                                                      >  1:2699622/1‑131 (MQ=255)
                                                    gACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTAc                                                                                                                      <  2:2699622/131‑1 (MQ=255)
                                                     aCCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCtt                                                                                                  <  1:1035477/150‑1 (MQ=255)
                                                       cAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTgg                                                                                                >  2:872231/1‑150 (MQ=255)
                                                       cAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTgg                                                                                                >  2:414950/1‑150 (MQ=255)
                                                         ggATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgcg                                                                                              <  1:872231/150‑1 (MQ=255)
                                                          gATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgcgc                                                                                             <  2:606790/150‑1 (MQ=255)
                                                              tAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCt                                                                                         <  2:1233071/150‑1 (MQ=255)
                                                              tAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCt                                                                                         >  2:1754360/1‑150 (MQ=255)
                                                                 tcattcatGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAg                                                                                                                                                      >  1:2341201/1‑86 (MQ=255)
                                                                 tcattcatGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAg                                                                                                                                                      <  2:2341201/86‑1 (MQ=255)
                                                                      catGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGc                                                                                 <  1:1027498/150‑1 (MQ=255)
                                                                       atGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATgcg                                                                                                              <  2:2177647/120‑1 (MQ=255)
                                                                       atGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATgcg                                                                                                              >  1:2177647/1‑120 (MQ=255)
                                                                          agaTGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTcag                                                                             <  2:1316947/150‑1 (MQ=255)
                                                                            aTGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTcagca                                                                           >  2:2033935/1‑150 (MQ=255)
                                                                              gACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAgg                                                                         >  1:422155/1‑150 (MQ=255)
                                                                              gACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAgg                                                                         <  2:1244003/150‑1 (MQ=255)
                                                                              gACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAgg                                                                         <  2:54634/150‑1 (MQ=255)
                                                                                    agacGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGAtt                                                                   >  2:920/4‑150 (MQ=255)
                                                                                      cggggAAGCGGGGGGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttc                                                                 <  2:123872/148‑1 (MQ=255)
                                                                                       cGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttct                                                                >  1:2188641/1‑150 (MQ=255)
                                                                                       cGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGAGGCGCTGGCCTGGCGTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttct                                                                >  1:2201117/1‑150 (MQ=255)
                                                                                         ggAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttcttc                                                              >  2:1302252/1‑150 (MQ=255)
                                                                                         ggAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttcttc                                                              >  1:2321672/1‑150 (MQ=255)
                                                                                         ggAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttcttc                                                              <  1:2033935/150‑1 (MQ=255)
                                                                                           aaTCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCAc                                                            <  1:414950/150‑1 (MQ=255)
                                                                                                gggCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctggcct                                                                                                       <  2:1058523/102‑1 (MQ=255)
                                                                                                gggCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctggcct                                                                                                       >  1:1058523/1‑102 (MQ=255)
                                                                                                    gATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGc                                                   >  2:339868/1‑150 (MQ=255)
                                                                                                         tCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCCGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAg                                                      <  1:2457822/142‑1 (MQ=255)
                                                                                                         tCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAg                                                      >  2:2457822/1‑142 (MQ=255)
                                                                                                          ccGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTACCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAg                                             >  2:2152420/1‑150 (MQ=255)
                                                                                                           cGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGc                                            >  1:2019043/1‑150 (MQ=255)
                                                                                                            gggTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGcc                                           <  2:2019043/150‑1 (MQ=255)
                                                                                                            gggTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGcc                                           <  1:1904930/150‑1 (MQ=255)
                                                                                                            gggTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCAGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGcc                                           >  1:537284/1‑150 (MQ=255)
                                                                                                                 tGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGAc                                      >  1:483007/1‑150 (MQ=255)
                                                                                                                 tGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGAc                                      >  1:2780111/1‑150 (MQ=255)
                                                                                                                   aCATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCa                                    >  1:2788530/1‑150 (MQ=255)
                                                                                                                   aCATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCa                                    <  2:314829/150‑1 (MQ=255)
                                                                                                                    cATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGAc                                      >  2:957143/1‑147 (MQ=255)
                                                                                                                    cATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGAc                                      <  1:957143/147‑1 (MQ=255)
                                                                                                                     aTGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATg                                  >  2:2313742/1‑150 (MQ=255)
                                                                                                                     aTGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATg                                  >  2:2522900/1‑150 (MQ=255)
                                                                                                                      tGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGc                                 >  2:3041349/1‑150 (MQ=255)
                                                                                                                          ggggCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGTTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTc                             <  2:1349905/150‑1 (MQ=255)
                                                                                                                          ggggCCTGC‑‑GTAATTGTGGGAGCGGGCAGCCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCCGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTc                             <  2:223677/150‑1 (MQ=255)
                                                                                                                           gggCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCg                            <  2:210103/150‑1 (MQ=255)
                                                                                                                           gggCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCg                            <  1:2993070/150‑1 (MQ=255)
                                                                                                                             gCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCCCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGa                          <  2:2208069/150‑1 (MQ=255)
                                                                                                                              ccTGC‑‑GTAATTGTGGGAGCGGGCAGGACTTTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGac                         <  1:920/150‑1 (MQ=255)
                                                                                                                                tGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGacac                       >  1:1912297/1‑150 (MQ=255)
                                                                                                                                 gC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACg                      >  1:775529/1‑150 (MQ=255)
                                                                                                                                        aTTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCaaa                 <  2:98751/150‑1 (MQ=255)
                                                                                                                                           gtgGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATg                                                                                                                <  2:523526/52‑1 (MQ=255)
                                                                                                                                           gtgGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATg                                                                                                                >  1:523526/1‑52 (MQ=255)
                                                                                                                                            tgGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCGGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAgcgc             >  1:2756525/1‑150 (MQ=255)
                                                                                                                                            tgGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAgcgc             >  1:2744266/1‑150 (MQ=255)
                                                                                                                                            tgGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGGTACTTGCTGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCGGGTAGATTTCTTCACAGCCCGGGCCGGCAGCCCTGTCCATGCCTTTGGCCACGCCAAAgcgc             >  1:1270811/1‑150 (MQ=255)
                                                                                                                                             gggAGCGGCCTGCCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCGTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGATCATGCCTTCGGACACGCCAAAGCGCa            >  1:488007/1‑150 (MQ=255)
                                                                                                                                               gAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATc          >  2:529986/1‑150 (MQ=255)
                                                                                                                                                   gggCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGTGCACGCCAAAGCGCATCTccc      >  2:27034/1‑150 (MQ=255)
                                                                                                                                                   gggCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTcc      >  2:1726151/1‑150 (MQ=255)
                                                                                                                                                       aGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCgggg  >  2:919808/1‑150 (MQ=255)
                                                                                                                                                       aGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCgggg  <  2:2877859/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AGAGTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCGGGG  >  NC_002947/1253288‑1253586

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