Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,273,229 A→T 8.5% intergenic (+27/+17) chaC → / ← ychN cation transport regulator/putative sulfur relay protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,273,2290AT 8.5% 187.4 / 7.2 71intergenic (+27/+17)chaC/ychNcation transport regulator/putative sulfur relay protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base T (3/3);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

CTGGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACC  >  NC_000913/1273080‑1273378
                                                                                                                                                     |                                                                                                                                                     
cTGTAACAGGAGCTTATCAAACTGGGAATGCAGGACGATTGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCa                                                                                                                                                       <  1:900900/150‑1 (MQ=255)
cTGGAACAGGAGCTTATCAAACTGGGAATCCAGGACGATGGGCTGAATGATTTGCGGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCa                                                                                                                                                       <  2:1492898/150‑1 (MQ=255)
  ggAACAGGAGATTATCAAACTGGGTATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGAGTTACGTCCGGGATTCGCCTTAGAAAACTTCCAGCATATATGGGCGTCAGa                                                                                                                                                     >  2:156869/1‑150 (MQ=255)
   gAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGAc                                                                                                                                                    >  2:610665/1‑150 (MQ=255)
   gAACAGCACCTTATCAAACTGGGAATCAAAGACGAAGCGCTTAATGATTTGCTGTTATCGGTAAAAAAACTGCTGGCGGAGAATTTCAGGATTGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACg                                                                                                                                                   >  2:284044/1‑150 (MQ=255)
    aaCAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATTGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACg                                                                                                                                                   <  2:353645/150‑1 (MQ=255)
     aCAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGc                                                                                                                                                  <  1:1642835/150‑1 (MQ=255)
          aGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGGTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCt                                                                                                                                             <  1:244300/150‑1 (MQ=255)
           gCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTc                                                                                                                                            >  1:1220499/1‑150 (MQ=255)
           gCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCGGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccact                                                                                                                                            <  2:456551/150‑4 (MQ=255)
            cTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccacta                                                                                                                                           <  2:784904/150‑5 (MQ=255)
            cTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCa                                                                                                                                           >  2:1479881/1‑150 (MQ=255)
                tCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACa                                                                                                                                       <  2:362361/150‑1 (MQ=255)
                   aaCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGATAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCTCAAACATTaa                                                                                                                                   <  1:1353635/150‑1 (MQ=255)
                   aaCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTACGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCACAAACATTa                                                                                                                                    >  1:1339026/1‑150 (MQ=255)
                    aCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTaa                                                                                                                                   >  1:971141/1‑150 (MQ=255)
                        ggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtg                                                                                                                               >  1:1440597/1‑150 (MQ=255)
                        ggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtg                                                                                                                               <  2:1440597/150‑1 (MQ=255)
                         gAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                              <  2:408543/150‑1 (MQ=255)
                           aTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCGGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGa                                                                                                                            >  1:1644954/1‑150 (MQ=255)
                            tGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAg                                                                                                                           <  2:1301266/150‑1 (MQ=255)
                              cAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTaaa                                                                                                                                  <  2:10731/141‑1 (MQ=255)
                              cAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTaaa                                                                                                                                  >  1:10731/1‑141 (MQ=255)
                              cAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCa                                                                                                                         >  2:55258/1‑150 (MQ=255)
                              cAGGACGATGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCAc                                                                                                                        >  2:1174492/1‑150 (MQ=255)
                                    gATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTAt                                                                                                                   <  2:626340/150‑1 (MQ=255)
                                       ttgcggtaTGATTTGCTGGTTTCGGTAAAAAAAGTGCTGGCGCAGAAATTTCCGGATGGTTTGTTACGTCCGGTATTAGCCTGATTAAACTTCCCGCATTGTGGGGCGTCAGACGCCCCCTCAAACCTTAAAATGTGAGGACTTTATAgg                                                                                                                <  2:28823/143‑1 (MQ=255)
                                        ggCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGc                                                                                                               >  1:866846/1‑150 (MQ=255)
                                        ggCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGc                                                                                                               <  1:415338/150‑1 (MQ=255)
                                         gCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                              <  2:1540561/150‑1 (MQ=255)
                                         gCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                              <  2:1702444/150‑1 (MQ=255)
                                            gAATGATATGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTCAGTAAACTTACCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGGACTTTATctgaagaa                                                                                                           >  1:564676/1‑143 (MQ=255)
                                              aTGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAg                                                                                                         >  1:1691698/1‑150 (MQ=255)
                                                   ttGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTa                                                                                                                    <  1:1040522/134‑1 (MQ=255)
                                                   ttGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTa                                                                                                                    >  2:1040522/1‑134 (MQ=255)
                                                       tGGTCTCGGTAACAAACCTGCTGGCGGAGTATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGGAATCTTCCCGCTTAGTGGGGCGTCATACGCCCCCTCAAACAATTAATTGTGAGCACTTTATCGGCGGACAGCGTCCTTTg                                                                                                <  1:1160374/150‑1 (MQ=255)
                                                             cGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACCGCGTCCATTGCGCCAg                                                                                          >  1:390593/1‑150 (MQ=255)
                                                             cGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAg                                                                                          <  1:1174492/150‑1 (MQ=255)
                                                             cGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAg                                                                                          <  1:420934/150‑1 (MQ=255)
                                                               gTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGACTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCACACATAAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGtt                                                                                        >  2:1017247/1‑150 (MQ=255)
                                                                tAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTc                                                                                       >  1:1599855/1‑150 (MQ=255)
                                                                  aaaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCca                                                                                     <  2:1713594/150‑1 (MQ=255)
                                                                   taaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCcac                                                                                    >  1:656314/2‑150 (MQ=255)
                                                                   aaaaaCTGCCGGCGGAGAATTTTCCGTAAGGTGTGTTACGACCGGGATTCGCCTGAGTAATCTTCCCGCATAGTGGGGCGTCAGACGCCCCCGCAAACATTAAAATGGGACCACTTTTTCGGCTGACAGCGTACATTGCGCCAGTTCcac                                                                                    <  1:162454/150‑1 (MQ=255)
                                                                        ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAagag                                                                               <  2:631263/150‑1 (MQ=255)
                                                                          gctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTa                                                                             <  2:1691698/150‑1 (MQ=255)
                                                                              gCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGa                                                                         <  2:927059/150‑1 (MQ=255)
                                                                                   gaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTcc                                                                    >  1:1357216/1‑150 (MQ=255)
                                                                                    aaTTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCa                                                                   <  2:1220499/150‑1 (MQ=255)
                                                                                                ggTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatcaatcc                                                       >  2:199269/1‑149 (MQ=255)
                                                                                                  tgtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCaga                                                     <  1:610665/150‑1 (MQ=255)
                                                                                                            ccGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGCGGCCCCCTGAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGGGCATTGAGCCAGTTCCCCACGTGTCCCGCTGTGCCCCCCAGCAATCCGAGGAAGTGGAc                                           >  2:308570/1‑150 (MQ=255)
                                                                                                             cGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACt                                          <  2:419741/150‑1 (MQ=255)
                                                                                                              gtgATTCGCCTGAGTAAACTACCCGCATTGTTGGTCTTCAGACGCACACTCAAACATTAAAATGTGAGCTCTTTTTCGGCTGACTGCGTCCATTGCGCCAGTTCCACAAGAGTACCTATTTCAACCCCATCAATCAGAGGAAGTGTACTa                                         <  1:156869/148‑1 (MQ=255)
                                                                                                               ggATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTaa                                        >  1:408291/1‑150 (MQ=255)
                                                                                                               ggATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTaa                                        <  2:1434728/150‑1 (MQ=255)
                                                                                                                gATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAAt                                       <  2:1079747/150‑1 (MQ=255)
                                                                                                                 aTTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATc                                      >  1:87254/1‑150 (MQ=255)
                                                                                                                 aTTCGCCGGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCATCAATCAGAGGAAGTGTACTAATcc                                     <  2:1644954/150‑1 (MQ=255)
                                                                                                                        tGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGcc                               <  2:971141/150‑1 (MQ=255)
                                                                                                                        tGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGcc                               >  2:1168053/1‑150 (MQ=255)
                                                                                                                        tGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGcc                               <  1:55258/150‑1 (MQ=255)
                                                                                                                        tGAGTAAACTTCCCGAATAGTGGGGCGTCAGACGCCCACTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATGAATCAGTGGAAGTGTACTAATCCCGCGcc                               <  1:284044/150‑1 (MQ=255)
                                                                                                                                   cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAg                    >  2:13394/1‑150 (MQ=255)
                                                                                                                                   cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACACCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAg                    <  2:1172492/150‑1 (MQ=255)
                                                                                                                                     cGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGt                  >  2:546747/1‑150 (MQ=255)
                                                                                                                                        ttgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGtttt               >  2:1263720/6‑150 (MQ=255)
                                                                                                                                         tgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTg              >  2:1695038/5‑150 (MQ=255)
                                                                                                                                         tAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGAGCCCGTCGGTACAGGTTTTg              <  2:1339026/150‑1 (MQ=255)
                                                                                                                                          aGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATATCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGc             <  2:259918/150‑1 (MQ=255)
                                                                                                                                           agggggCGTCTGACGCCCCCTCAAACATTAAAATGTGCGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCCCCCCCTCACTCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGcc            >  1:1059178/3‑149 (MQ=255)
                                                                                                                                             ggggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGcaca          >  2:684231/1‑150 (MQ=255)
                                                                                                                                                gCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTATAGGTTTTGCACAAtt       <  2:87254/150‑1 (MQ=255)
                                                                                                                                                     aGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAcc  <  2:866846/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACC  >  NC_000913/1273080‑1273378

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