Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,253,435 G→C 6.5% intergenic (‑75/+32) 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,4350GC 6.5% 291.5 / 10.0 92intergenic (‑75/+32)PP_1096/metGelectron transport complex protein/methionine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (47/39);  new base C (3/3);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GGAGAGTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCC  >  NC_002947/1253286‑1253582
                                                                                                                                                       |                                                                                                                                                   
ggAGGGTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCggg                                                                                                                                                     <  1:1133808/150‑1 (MQ=255)
  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)
                     gCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCggg                                                                                                                                                     <  2:719851/129‑1 (MQ=255)
                     gCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCggg                                                                                                                                                     >  1:719851/1‑129 (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                                                                                                                         >  2:554983/1‑143 (MQ=255)
                                   tCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAggg                                                                                                                         <  1:554983/143‑1 (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                                                                                                                                       >  1:2088065/1‑114 (MQ=255)
                                                  tCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCcac                                                                                                                                       <  2:2088065/114‑1 (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:414950/1‑150 (MQ=255)
                                                         cAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTgg                                                                                            >  2:872231/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                                                                                                                                                  <  2:2341201/86‑1 (MQ=255)
                                                                   tcattcatGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAg                                                                                                                                                  >  1:2341201/1‑86 (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                                                          <  1:2033935/150‑1 (MQ=255)
                                                                                           ggAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttcttc                                                          >  2:1302252/1‑150 (MQ=255)
                                                                                           ggAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATttcttc                                                          >  1:2321672/1‑150 (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:2780111/1‑150 (MQ=255)
                                                                                                                   tGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGAc                                  >  1:483007/1‑150 (MQ=255)
                                                                                                                     aCATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCa                                <  2:314829/150‑1 (MQ=255)
                                                                                                                     aCATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCa                                >  1:2788530/1‑150 (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)
                                                                                                                                                       |                                                                                                                                                   
GGAGAGTGGCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGC‑‑GTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCC  >  NC_002947/1253286‑1253582

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: