Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,273,229 A→T 8.3% 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.3% 195.3 / 6.4 59intergenic (+27/+17)chaC/ychNcation transport regulator/putative sulfur relay protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base T (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAC  >  NC_000913/1273085‑1273377
                                                                                                                                                 |                                                                                                                                                    
aCAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACg                                                                                                                                                  >  3:120447/1‑149 (MQ=255)
aCAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACg                                                                                                                                                  >  1:341379/1‑149 (MQ=255)
    gAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGGCGTCTGACGccc                                                                                                                                               <  3:214014/149‑1 (MQ=255)
          aTCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCaaa                                                                                                                                        <  3:112385/149‑1 (MQ=255)
          aTCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCaaa                                                                                                                                        >  4:788150/1‑149 (MQ=255)
              aaCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACAtt                                                                                                                                    >  2:574703/1‑149 (MQ=255)
              aaCTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACAtt                                                                                                                                    <  1:602569/149‑1 (MQ=255)
                 tGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTaaa                                                                                                                                 >  3:153834/1‑149 (MQ=255)
                  gggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTaaaa                                                                                                                                <  4:697168/149‑1 (MQ=255)
                    gAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAAtg                                                                                                                              <  1:593234/149‑1 (MQ=255)
                                aTGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTAt                                                                                                                  >  2:323964/1‑149 (MQ=255)
                                 tGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATc                                                                                                                 <  3:360913/149‑1 (MQ=255)
                                 tGGGCTGAATGATTTGCTGGTATCGGT‑AAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCg                                                                                                                <  1:376305/149‑1 (MQ=255)
                                     cTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                             <  2:167753/149‑1 (MQ=255)
                                      tGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTg                                                                                                            >  2:46462/1‑149 (MQ=255)
                                       gAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGa                                                                                                           >  1:488165/1‑149 (MQ=255)
                                        aaTGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑‑GGGCGTCTGACGcc                                                                                                                                                <  2:620778/110‑1 (MQ=255)
                                        aaTGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑‑GGGCGTCTGACGcc                                                                                                                                                >  1:620778/1‑110 (MQ=255)
                                          tGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAg                                                                                                        >  1:6522/1‑149 (MQ=255)
                                              ttGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTc                                                                                                    <  2:470492/149‑1 (MQ=255)
                                                    gTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgc                                                                                              >  2:781465/1‑149 (MQ=255)
                                                    gTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgc                                                                                              <  2:40142/149‑1 (MQ=255)
                                                    gTATCGGTAAAAAAACTGCTGGCGAAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTGCCCGGATAGT‑GGGGGGTCAGACGCCCCCACAAACATTAGAATGTGAGCGCTTTATCGGCTGACAGCGTCAATTgc                                                                                              <  4:359811/149‑1 (MQ=255)
                                                         ggTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAg                                                                                         >  1:542980/1‑149 (MQ=255)
                                                          gTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCgg                                                                                                               >  4:429503/1‑126 (MQ=255)
                                                          gTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCgg                                                                                                               <  3:429503/126‑1 (MQ=255)
                                                          gTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGt                                                                                        >  4:344593/1‑149 (MQ=255)
                                                           tAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGtt                                                                                       >  3:342454/1‑149 (MQ=255)
                                                              aaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCca                                                                                    >  4:659289/1‑149 (MQ=255)
                                                               aaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                             >  2:259797/1‑107 (MQ=255)
                                                               aaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                             <  1:259797/107‑1 (MQ=255)
                                                                 aaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGGCGCCCCCTCAAACATTAAAATGTGAGCAATTTATCGGCTGACAGCGTCCATTGCGCCAGTGCCACaa                                                                                 <  2:170822/149‑1 (MQ=255)
                                                                 aaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACaa                                                                                 <  3:113963/149‑1 (MQ=255)
                                                                   ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAaga                                                                               <  3:671969/149‑1 (MQ=255)
                                                                       tgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTATAATGTGAGCACTTTATCGGCTGATAGCGTCCATTGCGCCAGTTCCACAAGAGTAc                                                                           >  1:253485/1‑149 (MQ=255)
                                                                              gaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTc                                                                    >  2:101046/1‑149 (MQ=255)
                                                                                 ttttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAc                                                                 <  1:46462/149‑1 (MQ=255)
                                                                                  tttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                                <  4:120447/149‑1 (MQ=255)
                                                                                      cGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCat                                                            <  3:287879/149‑1 (MQ=255)
                                                                                        gATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatca                                                          >  3:523485/1‑149 (MQ=255)
                                                                                          tGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatcaat                                                        >  4:716603/1‑149 (MQ=255)
                                                                                             tgtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCag                                                     >  3:563716/1‑149 (MQ=255)
                                                                                              gtgtTACGTCCGGGAGTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCaga                                                    <  3:233505/149‑1 (MQ=255)
                                                                                              cgatctcGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCaga                                                    <  3:538517/143‑1 (MQ=255)
                                                                                                 ttACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTACCACCCCATCAATCAGAGGa                                                 >  3:774608/1‑149 (MQ=255)
                                                                                                  tACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGaa                                                <  2:53663/149‑1 (MQ=255)
                                                                                                      tCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAgtgt                                            >  3:268432/1‑149 (MQ=255)
                                                                                                       ccGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTa                                           <  3:788150/149‑1 (MQ=255)
                                                                                                                gCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCcg                                  >  1:186020/1‑149 (MQ=255)
                                                                                                                  cTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCcgcg                                >  1:535049/1‑149 (MQ=255)
                                                                                                                     aGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGccc                             >  1:280076/1‑149 (MQ=255)
                                                                                                                          aCTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCgg                        <  4:316441/149‑1 (MQ=255)
                                                                                                                            ttCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTa                      <  3:344593/149‑1 (MQ=255)
                                                                                                                               ccGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAg                   <  4:563716/149‑1 (MQ=255)
                                                                                                                                  cATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGtt                <  2:186020/149‑1 (MQ=255)
                                                                                                                                   tttgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAg                   <  2:586184/140‑1 (MQ=255)
                                                                                                                                   tttgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAg                   >  1:586184/7‑146 (MQ=255)
                                                                                                                                         ggggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGcac          >  3:288776/1‑149 (MQ=255)
                                                                                                                                           ggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACaa        <  3:44392/149‑1 (MQ=255)
                                                                                                                                                 aGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAc  <  2:320947/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGT‑GGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCAC  >  NC_000913/1273085‑1273377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑