Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,377,112 T→C 7.9% intergenic (‑48/‑71) 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,377,1120TC 7.9% 876.9 / 43.3 292intergenic (‑48/‑71)csoR/PP_2970copper utilization transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (153/113);  new base C (9/14);  total (163/129)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 293 total aligned reads displayed.

GATCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGCACACC  >  NC_002947/3376963‑3377260
                                                                                                                                                     |                                                                                                                                                     
 aTCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTg                                                                                                                                                     >  2:6493296/1‑150 (MQ=255)
 aTCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCACGCATCCTATCCCCCTg                                                                                                                                                     >  1:5926326/1‑150 (MQ=255)
   cATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTggg                                                                                                                                                   <  2:3871325/150‑1 (MQ=255)
    aTGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTgggg                                                                                                                                                  <  1:3139547/150‑1 (MQ=255)
     tGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTggggg                                                                                                                                                 <  1:2885871/150‑1 (MQ=255)
     tGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTggggg                                                                                                                                                 <  1:4947318/150‑1 (MQ=255)
      gCTGACGCTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTGTTGATGATCTCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTgggggg                                                                                                                                                <  1:4676992/150‑1 (MQ=255)
          aCGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGata                                                                                                                                            <  1:6157744/150‑1 (MQ=255)
          aCGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGata                                                                                                                                            <  1:3048301/150‑1 (MQ=255)
            gATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatagg                                                                                                                                          <  1:5849897/150‑3 (MQ=255)
             aTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATgg                                                                                                                                         <  1:5893962/150‑1 (MQ=255)
             aTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatagga                                                                                                                                         <  2:6005350/150‑4 (MQ=255)
             aTGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCCCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCCCAGGGGGatagga                                                                                                                                         <  1:6365882/150‑4 (MQ=255)
              tgctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGt                                                                                                                                        <  1:102247/150‑1 (MQ=255)
              gcgtgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGt                                                                                                                                        >  1:2751120/4‑150 (MQ=255)
               gctgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCATATCCCCCCCAGGGGGataggatn                                                                                                                                       <  1:2250996/150‑6 (MQ=255)
                 tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaa                                                                                                                                     >  1:4120713/1‑150 (MQ=255)
                 tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaa                                                                                                                                     >  2:5713007/1‑150 (MQ=255)
                 tgcGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaa                                                                                                                                     >  1:6326747/1‑150 (MQ=255)
                  gcgcAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCaaa                                                                                                                                    >  2:6863551/1‑150 (MQ=255)
                      aGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCaa                                                                                                                                <  1:2389815/150‑1 (MQ=255)
                        ggccaccTGCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACt                                                                                                                              >  1:5992093/6‑150 (MQ=255)
                           gCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTAcc                                                                                                                           <  1:3496125/150‑1 (MQ=255)
                             ccTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTAcccc                                                                                                                         >  2:2047288/1‑150 (MQ=255)
                             ccTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCATGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTAcccc                                                                                                                         >  1:4394710/1‑150 (MQ=255)
                              cTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCt                                                                                                                        >  1:6779682/1‑150 (MQ=255)
                                 cGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCa                                                                                                                     <  2:1745137/150‑1 (MQ=255)
                                    aCGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGaa                                                                                                                  >  1:4605517/1‑150 (MQ=255)
                                        ttCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGaacaac                                                                                                              >  2:2270297/1‑150 (MQ=255)
                                        ttCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGaacaac                                                                                                              >  2:871297/1‑150 (MQ=255)
                                           aaCCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAATCCAACTACCCCTTCAGAACAACCAt                                                                                                           >  2:1058621/1‑150 (MQ=255)
                                                  tgatgaTATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTa                                                                                                    <  2:2683705/150‑1 (MQ=255)
                                                     tgatATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGc                                                                                                 <  2:5619408/150‑1 (MQ=255)
                                                      gatATCGCTATGACTTTGGAGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCa                                                                                                >  2:4347881/1‑150 (MQ=255)
                                                         atCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCa                                                                                             >  2:4947050/1‑150 (MQ=255)
                                                          tCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaa                                                                                            <  2:4454518/150‑1 (MQ=255)
                                                          tCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCGGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaa                                                                                            >  1:5461689/1‑150 (MQ=255)
                                                           cGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaa                                                                                           <  2:3233102/150‑1 (MQ=255)
                                                           cGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaa                                                                                           <  1:5788217/150‑1 (MQ=255)
                                                           cGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaa                                                                                           <  2:904140/150‑1 (MQ=255)
                                                           cGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                          <  1:3331569/150‑1 (MQ=255)
                                                            gCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                          <  2:4776415/150‑1 (MQ=255)
                                                            gCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCaaaa                                                                                          >  1:3702779/1‑150 (MQ=255)
                                                             cTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAg                                                                                         >  1:6842598/1‑150 (MQ=255)
                                                             cTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAg                                                                                         <  1:4490575/150‑1 (MQ=255)
                                                              tATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGa                                                                                        >  2:3359785/1‑150 (MQ=255)
                                                                   cTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAacac                                                                                   >  2:3726792/1‑150 (MQ=255)
                                                                     ttGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGa                                                                                 <  2:1458082/150‑1 (MQ=255)
                                                                         tGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATg                                                                             <  2:1396907/150‑1 (MQ=255)
                                                                          gCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGc                                                                            >  2:1173030/1‑150 (MQ=255)
                                                                             gATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGTGTGTTATGTTCAAACCAACTACCCCTTAAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCa                                                                         <  1:1875865/150‑1 (MQ=255)
                                                                             gATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCTTCAAAAGAAACACGACATGCTCa                                                                         <  2:2521971/150‑1 (MQ=255)
                                                                             gATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCa                                                                         <  2:2540627/150‑1 (MQ=255)
                                                                                   cGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAAcc                                                                   >  2:3065094/1‑150 (MQ=255)
                                                                                    gtgtTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAAccc                                                                  >  1:5553265/1‑150 (MQ=255)
                                                                                          gTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAg                                                            >  2:2863042/1‑150 (MQ=255)
                                                                                            gATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCa                                                          <  2:5892277/150‑1 (MQ=255)
                                                                                             aTCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCAc                                                         >  1:1248311/1‑150 (MQ=255)
                                                                                             aTCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCAGCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCAc                                                         >  1:3408184/1‑150 (MQ=255)
                                                                                                 cTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAg                                                     <  1:3212497/150‑1 (MQ=255)
                                                                                                  tCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAgg                                                    >  2:5361860/1‑150 (MQ=255)
                                                                                                       ccTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGAtt                                               <  1:2523443/150‑1 (MQ=255)
                                                                                                           gAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTg                                           >  2:3443735/1‑150 (MQ=255)
                                                                                                            aGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACACCCATCCGCTTAGGCATCAAAAGAACCACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGc                                          <  1:1114913/150‑1 (MQ=255)
                                                                                                            aGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGc                                          <  2:6779682/150‑1 (MQ=255)
                                                                                                              tctcATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGccc                                        >  1:4056041/1‑150 (MQ=255)
                                                                                                                tcATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAAc                                                                                                                                   <  2:522228/57‑1 (MQ=255)
                                                                                                                tcATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCtt                                      >  1:4159504/1‑150 (MQ=255)
                                                                                                                 cATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTg                                     <  1:5153375/150‑1 (MQ=255)
                                                                                                                 cATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACAAGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTg                                     >  1:6093090/1‑150 (MQ=255)
                                                                                                                       tCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTgg                               <  2:638983/150‑1 (MQ=255)
                                                                                                                         gATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGcc                             >  1:1565010/1‑150 (MQ=255)
                                                                                                                         gATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACACCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGcc                             >  2:1343393/1‑150 (MQ=255)
                                                                                                                           tGAATTACGCTCGCATCCTATCCCCCTGGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCa                            >  1:3482409/1‑150 (MQ=255)
                                                                                                                           agAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCAt                           >  1:3911008/2‑150 (MQ=255)
                                                                                                                              aTTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCtt                        >  1:2979469/1‑150 (MQ=255)
                                                                                                                                tACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGa                      <  1:1959704/150‑1 (MQ=255)
                                                                                                                                  cGCTCGCATCCTATCCCCCGGGGGGGGATATGGTCAAACCCACTCCCCCTACCGAACACCCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATa                    <  1:3792441/150‑1 (MQ=255)
                                                                                                                                     tCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCAt                                                                                               <  1:1324659/72‑1 (MQ=255)
                                                                                                                                     tCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCAt                                                                                               >  2:1324659/1‑72 (MQ=255)
                                                                                                                                     tCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGcc                 >  2:4002155/1‑150 (MQ=255)
                                                                                                                                     tCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGcc                 >  2:6159921/1‑150 (MQ=255)
                                                                                                                                     tCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGcc                 >  2:53363/1‑150 (MQ=255)
                                                                                                                                     tCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGcc                 >  2:4228894/1‑150 (MQ=255)
                                                                                                                                      cGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCt                >  1:5660009/1‑150 (MQ=255)
                                                                                                                                       gCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTg               >  1:6200175/1‑150 (MQ=255)
                                                                                                                                        cATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTgg              <  2:2929807/150‑1 (MQ=255)
                                                                                                                                         accatATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTggg             >  1:3064624/5‑150 (MQ=255)
                                                                                                                                          ccatATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCCGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCCTCTGGAGAGCGTGGGc            >  2:4689280/4‑150 (MQ=255)
                                                                                                                                          ccatATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGc            >  2:3615442/4‑150 (MQ=255)
                                                                                                                                           catATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGcc           >  2:3848692/3‑150 (MQ=255)
                                                                                                                                           catATCCCCCCCAGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTGTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGcc           >  1:1797359/3‑150 (MQ=255)
                                                                                                                                            cTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCg          >  2:277711/1‑150 (MQ=255)
                                                                                                                                            cTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCg          >  1:858762/1‑150 (MQ=255)
                                                                                                                                             tATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCgg         >  1:6452402/1‑150 (MQ=255)
                                                                                                                                             tATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCgg         >  1:1426416/1‑150 (MQ=255)
                                                                                                                                              aTCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCggg        >  2:5093912/1‑150 (MQ=255)
                                                                                                                                                  cccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGcaca    >  1:2829568/1‑150 (MQ=255)
                                                                                                                                                  cccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGcaca    >  2:5130818/1‑150 (MQ=255)
                                                                                                                                                  cccTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGcaca    >  1:1370390/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCATGCTGACGATGCTGCGCAGGTGGCCCTCCGCACGCTTCAACCTTTTGATGATATCGCTATGACTTTGGTGCTGATGGCCGTGTTCGTGATCGCTCATGCCTTGAGTCTCATGCATCGATGAATTACGCTCGCATCCTATCCCCCTGGGGGGGATATGGTCAAACCAACTACCCCTTCAGAACAACCATCCGCTTAGGCATCAAAAGAAACACGACATGCTCAGCAACCCCATCAGCACTAAGGTGATTTTTGCCCTTGCTCTGGCCATCTTGATAGCCTGGGCCGGGCACACC  >  NC_002947/3376963‑3377260

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