Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,374,086 G→C 9.1% intergenic (‑49/‑70) csoR ← / → PP_2970 copper utilization transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,374,0860GC 9.1% 634.8 / 40.9 202intergenic (‑49/‑70)csoR/PP_2970copper utilization transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (84/96);  new base C (5/13);  total (90/112)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 213 total aligned reads displayed.

TCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGCACACCT  >  NC_002947/3373938‑3374234
                                                                                                                                                     |                                                                                                                                                    
tCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTgg                                                                                                                                                     <  1:3280180/150‑1 (MQ=255)
  aTGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTgggg                                                                                                                                                   <  1:2622617/150‑1 (MQ=255)
  aTGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTgggg                                                                                                                                                   <  1:690023/150‑1 (MQ=255)
   tGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTggggg                                                                                                                                                  <  2:1482667/150‑1 (MQ=255)
    gCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTgggggg                                                                                                                                                 <  1:2053276/150‑1 (MQ=255)
      tGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTggggg                                                                                                                                                  >  1:3603266/1‑147 (MQ=255)
      tGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTggggg                                                                                                                                                  <  2:3603266/147‑1 (MQ=255)
      tGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGa                                                                                                                                               >  2:3715580/1‑150 (MQ=255)
      tGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGa                                                                                                                                               >  1:3841419/1‑150 (MQ=255)
       gACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGat                                                                                                                                              <  1:4867462/150‑1 (MQ=255)
         cGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatag                                                                                                                                            <  1:1068342/150‑2 (MQ=255)
          gATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatagg                                                                                                                                           <  2:2434536/150‑3 (MQ=255)
           aTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATgg                                                                                                                                          <  2:4363110/150‑1 (MQ=255)
           aTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatagga                                                                                                                                          <  1:4303001/150‑4 (MQ=255)
            tgctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCCGGGGGataggat                                                                                                                                         <  2:1328620/150‑5 (MQ=255)
            tgctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCATATCCCCCCCAGGGGGataggat                                                                                                                                         <  2:2173670/150‑5 (MQ=255)
            tgctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCATATCCCCCCCAGGGGGataggat                                                                                                                                         <  1:1953278/150‑5 (MQ=255)
             gctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTc                                                                                                                                        <  1:2694113/150‑1 (MQ=255)
              ctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCa                                                                                                                                       >  1:4506617/1‑150 (MQ=255)
               tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCTCCTGGGGGGGATATGGTCaa                                                                                                                                      >  2:5028500/1‑150 (MQ=255)
               tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaa                                                                                                                                      >  1:839122/1‑150 (MQ=255)
               tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaa                                                                                                                                      >  2:3293717/1‑150 (MQ=255)
               tccgcAGGTGCCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAAGTACGCGCGCATCCAATCCCCATGGGGGGGATATGGTCaa                                                                                                                                      >  2:538310/3‑150 (MQ=255)
                gcgcAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaaa                                                                                                                                     <  2:127727/150‑1 (MQ=255)
                gcgcAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaaa                                                                                                                                     >  1:2720284/1‑150 (MQ=255)
                  gcAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAAcc                                                                                                                                   >  1:3100381/1‑150 (MQ=255)
                   cAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCa                                                                                                                                  >  1:2883472/1‑150 (MQ=255)
                           ccTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTAcccc                                                                                                                          >  1:59327/1‑150 (MQ=255)
                           ccTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGATATGGTCAAACCAACTACCCCt                                                                                                                         >  1:4165224/1‑150 (MQ=255)
                             tCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCtt                                                                                                                        <  2:3022319/150‑1 (MQ=255)
                               cGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCa                                                                                                                      <  1:4061530/150‑1 (MQ=255)
                                 cACGCTTCAACCTTTTGATGATATCGCTATGACTATGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGa                                                                                                                    <  2:2720284/150‑1 (MQ=255)
                                     cTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGaacaa                                                                                                                <  1:325034/150‑1 (MQ=255)
                                      ttCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGGATATGGTCAAACCAACTACCCCTTCAGaacaa                                                                                                                <  2:2059707/150‑1 (MQ=255)
                                      ttCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATga                                                                                                                                          >  1:2283156/1‑122 (MQ=255)
                                      ttCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATga                                                                                                                                          <  2:2283156/123‑2 (MQ=255)
                                      ttCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGaacaac                                                                                                               >  2:2794597/1‑150 (MQ=255)
                                           ccTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATcc                                                                                                          <  2:4234893/150‑1 (MQ=255)
                                            cTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCg                                                                                                         >  1:2531035/1‑150 (MQ=255)
                                            cTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCg                                                                                                         >  1:4748176/1‑150 (MQ=255)
                                              tttGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCa                                                                                                                      <  1:2802708/135‑1 (MQ=255)
                                               ttGATGATATCGCTATGACTTTGGTGCTGACGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCtt                                                                                                      <  1:1444214/150‑1 (MQ=255)
                                                tgatgaTATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTa                                                                                                     <  2:1606854/150‑1 (MQ=255)
                                                       atCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCa                                                                                              >  2:5450725/1‑150 (MQ=255)
                                                        tCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGat                                                                                                                                              >  2:4529385/1‑101 (MQ=255)
                                                         cGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaa                                                                                            <  1:237368/150‑1 (MQ=255)
                                                          gCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                           <  1:1371587/150‑1 (MQ=255)
                                                          gCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCAGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                           <  2:1598069/150‑1 (MQ=255)
                                                           cTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAg                                                                                          <  2:416570/150‑1 (MQ=255)
                                                                aCTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAaca                                                                                     <  1:4610677/150‑1 (MQ=255)
                                                                 cTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatagg                                                                                                                                           >  1:5199757/1‑93 (MQ=255)
                                                                  tttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACg                                                                                   >  2:3506184/1‑150 (MQ=255)
                                                                  tttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACg                                                                                   >  2:71110/1‑150 (MQ=255)
                                                                  tttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACg                                                                                   >  1:4076039/1‑150 (MQ=255)
                                                                  tttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGataggat                                                                                                                                         <  2:344465/96‑5 (MQ=255)
                                                                  tttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCATATCCCCCCCAGGGGGatagga                                                                                                                                          >  1:1456592/1‑92 (MQ=255)
                                                                  tttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCATATCCCCCCCAGGGGGatagga                                                                                                                                          <  2:1456592/95‑4 (MQ=255)
                                                                   ttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGa                                                                                  <  1:4963698/150‑1 (MQ=255)
                                                                     ggTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaaa                                                                                                                                     >  1:3442220/1‑97 (MQ=255)
                                                                      gTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACAt                                                                               <  2:4756676/150‑1 (MQ=255)
                                                                       tGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTAc                                                                                                                             >  2:2352187/1‑103 (MQ=255)
                                                                          tGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTa                                                                                                     <  1:2470415/124‑1 (MQ=255)
                                                                          tGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTc                                                                           <  2:394121/150‑1 (MQ=255)
                                                                            aTGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAg                                                                         <  1:4989515/150‑1 (MQ=255)
                                                                            aTGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAg                                                                         <  2:4876210/150‑1 (MQ=255)
                                                                                 cGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAAcc                                                                    >  2:4053009/1‑150 (MQ=255)
                                                                                  gtgtTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAAccc                                                                   >  2:1124924/1‑150 (MQ=255)
                                                                                         tGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAAc                                                                                                                                    >  2:2958137/1‑78 (MQ=255)
                                                                                         tGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGc                                                            >  2:4521741/1‑150 (MQ=255)
                                                                                         tGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGc                                                            <  1:4025815/150‑1 (MQ=255)
                                                                                         tGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGc                                                            >  2:1681434/1‑150 (MQ=255)
                                                                                           aTCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCAc                                                          >  2:5073230/1‑150 (MQ=255)
                                                                                             cGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTa                                                        >  2:2107683/1‑150 (MQ=255)
                                                                                               cTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAg                                                      <  1:5276765/150‑1 (MQ=255)
                                                                                               cTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAg                                                      >  2:2318847/1‑150 (MQ=255)
                                                                                                           gTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGcc                                          <  1:1616849/150‑1 (MQ=255)
                                                                                                             ctcATGCATCGATGAATTACGCTCGCATCCTATCCTCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCt                                        <  1:5028500/150‑1 (MQ=255)
                                                                                                                aTGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCa                                                                                                                      <  1:1371972/69‑1 (MQ=255)
                                                                                                                aTGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCa                                                                                                                      >  2:1371972/1‑69 (MQ=255)
                                                                                                                aTGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGc                                     <  1:962432/150‑1 (MQ=255)
                                                                                                                             ttACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACGACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTg                        >  1:2031311/1‑150 (MQ=255)
                                                                                                                              tACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGa                       <  2:5200040/150‑1 (MQ=255)
                                                                                                                               aCGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGAt                      <  1:3902461/150‑1 (MQ=255)
                                                                                                                                cGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATa                     <  2:4326296/150‑1 (MQ=255)
                                                                                                                                      cATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTgg               >  2:3365186/1‑150 (MQ=255)
                                                                                                                                      cATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTgg               >  1:4693731/1‑150 (MQ=255)
                                                                                                                                      cATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTgg               <  2:1937633/150‑1 (MQ=255)
                                                                                                                                       accatATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTggg              >  1:1538503/5‑150 (MQ=255)
                                                                                                                                        tCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGc             >  2:1694049/1‑150 (MQ=255)
                                                                                                                                        ccatATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGc             >  1:1751266/4‑150 (MQ=255)
                                                                                                                                        ccatATCCCCCCCAGGGGGATAGGGTCAATCCAACTAGCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCAGTAAGGTGATTTTTGCCGTTGCTCTGGCGATCTTGATTGCGTgggg             >  2:3012814/4‑149 (MQ=255)
                                                                                                                                           tATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCgg          >  1:1186596/1‑150 (MQ=255)
                                                                                                                                           tATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCgg          >  1:3488216/1‑150 (MQ=255)
                                                                                                                                           tATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCgg          <  2:2869851/150‑1 (MQ=255)
                                                                                                                                            ctCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTACAGAACAACCCTCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCggg         <  1:709646/149‑1 (MQ=255)
                                                                                                                                            aTCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCggg         >  2:4965646/1‑150 (MQ=255)
                                                                                                                                             tCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGc        >  1:1930801/1‑150 (MQ=255)
                                                                                                                                              cccccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGca       <  1:3365186/150‑1 (MQ=255)
                                                                                                                                               ccccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGcac      <  2:2639390/150‑1 (MQ=255)
                                                                                                                                                cccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGcaca     <  1:1248846/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGCACACCT  >  NC_002947/3373938‑3374234

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