Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,019 G→T 5.8% intergenic (+206/+26) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,918,0190GT 5.8% 347.0 / 7.5 120intergenic (+206/+26)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base G (62/51);  new base T (2/5);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 122 total aligned reads displayed.

CAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAC  >  NC_000913/2917870‑2918167
                                                                                                                                                      |                                                                                                                                                    
cAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAg                                                                                                                                                      >  1:1018917/1‑150 (MQ=255)
        cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGcc                                                                                                                                              <  1:802451/150‑1 (MQ=255)
           tttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggt                                                                                                                                           >  1:443609/8‑150 (MQ=39)
           tgtTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggt                                                                                                                                           <  1:1125865/150‑1 (MQ=38)
            ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACACCAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtt                                                                                                                                          >  2:121230/7‑150 (MQ=39)
            gtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtt                                                                                                                                          >  2:587266/1‑150 (MQ=255)
              ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                        <  1:230647/150‑2 (MQ=255)
               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttg                                                                                                                                       <  1:802703/150‑1 (MQ=255)
               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggatttg                                                                                                                                       >  2:1663242/1‑150 (MQ=255)
                aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttgag                                                                                                                                      <  2:1492587/150‑4 (MQ=38)
                  tccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGCAGCCggtttgagac                                                                                                                                    <  1:1686080/149‑6 (MQ=38)
                  gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAAGGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtt                                                                                                                                    >  1:450374/1‑150 (MQ=255)
                  gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtt                                                                                                                                    >  1:428547/1‑150 (MQ=255)
                     ggataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGCCGGTTTTGGTTTTc                                                                                                                                 >  1:366452/7‑150 (MQ=255)
                      gATGCGGAGTTAACGCCTTATCCGGCCTACAAAAACAATTAAATTCAATGAATTTAAAAAATGATTTAGACCTGATACATAGGACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                 <  2:1539607/150‑1 (MQ=38)
                       aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCATTGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAAGAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                               <  2:126199/150‑1 (MQ=255)
                       aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                               <  1:746967/150‑1 (MQ=255)
                       aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                               <  1:1717997/150‑1 (MQ=255)
                         tcggcgTTAACTCGTTATCCGGCCTACATAACCACTTAAATTCAATTAA‑TTAAAAAATGTTGTAGACCTGATACATCGCACAACAGGCTATTTTGCGTCTAGCATCAGTCTCAACCCGGCTCCAGATTGTGACGGTTTTGGTTTTCTGTc                                                                                                                             <  1:950412/149‑1 (MQ=255)
                         gcggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTc                                                                                                                             <  2:176055/150‑1 (MQ=255)
                          cggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                            >  2:962261/1‑150 (MQ=255)
                            gcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTa                                                                                                                          >  2:1622191/1‑150 (MQ=255)
                               ttAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACg                                                                                                                       >  2:1733755/1‑150 (MQ=255)
                                 aaCGCCTTATCCGGCCTACAAACCCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                     >  2:1143581/1‑150 (MQ=255)
                                 aaCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                     >  1:684265/1‑150 (MQ=255)
                                      cTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                                <  2:1398068/150‑1 (MQ=255)
                                      cTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                                <  1:1326738/150‑1 (MQ=255)
                                       ttATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTATCTGTCTTAACGCACCAtcc                                                                                                               >  2:866339/1‑148 (MQ=255)
                                        tATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAACACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGATCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCa                                                                                                              >  2:543591/1‑150 (MQ=255)
                                          tCCGGCCTACAAAACCAATTAACTTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACg                                                                                                            <  2:1510666/150‑1 (MQ=255)
                                           ccGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                           <  1:1263876/150‑1 (MQ=255)
                                           ccGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                           >  2:1476940/1‑150 (MQ=255)
                                               ccTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                       >  2:794828/1‑150 (MQ=255)
                                                 tACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACAATGCACGGGCGCtt                                                                                                     >  2:78147/1‑150 (MQ=255)
                                                 tACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGCCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATTCACGGGCGCtt                                                                                                     >  2:1323562/1‑150 (MQ=255)
                                                  aCAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  1:587266/150‑1 (MQ=255)
                                                  aCAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  2:635113/150‑1 (MQ=255)
                                                  aCAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  1:1492480/150‑1 (MQ=255)
                                                  aCAAAAACAATTAAATTCACTGAA‑TTAAAAAATGATGTAGACCGCCTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  2:1154745/150‑1 (MQ=255)
                                                    aaaaCCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGtt                                                                                                  >  2:1398798/1‑150 (MQ=255)
                                                     aaaCCAATTAAATTCAATTAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGGCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                 <  2:248991/150‑1 (MQ=255)
                                                     aaaCCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                 <  1:625681/150‑1 (MQ=255)
                                                      aaCCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGTTACATAGCACAACAGGATATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCTGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCTCGGGCGCTTGTTAt                                                                                                <  1:826095/150‑1 (MQ=255)
                                                      aaCCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                <  2:1455475/150‑1 (MQ=255)
                                                          aaTTAAATTCAATGAA‑TTAAAAAATGATGTAGCCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTGTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAGCGac                                                                                            >  2:709984/1‑149 (MQ=255)
                                                              aaaTTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACACGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTc                                                                                        <  1:47906/150‑1 (MQ=255)
                                                                 ttCAATGAA‑TTAAAAAATGATGTAGCCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccag                                                                                     >  1:1137850/1‑150 (MQ=255)
                                                                 ttCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGATCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTAGTTATCCAACGTccac                                                                                     >  2:1345519/1‑149 (MQ=255)
                                                                 ttCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccag                                                                                     <  1:440556/150‑1 (MQ=255)
                                                                    aaTGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                  >  2:805058/1‑150 (MQ=255)
                                                                     aTGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                 <  1:1114697/150‑1 (MQ=255)
                                                                      tGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTATTAACGAACCATGCACGGGCGATTGTTATGGAACGTCCAGCCAgg                                                                                >  1:239069/1‑150 (MQ=255)
                                                                        aa‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaa                                                                              >  1:659180/1‑150 (MQ=255)
                                                                           ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCATGTTATCGAACGTCCAGCCAGGAATc                                                                            >  1:1449994/1‑150 (MQ=255)
                                                                           ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGGTTTGGTTTTCTGGCTTAACGCACCATGCACGTCGGCTGGTTATCGAACGTCAAGCAAGTAATc                                                                            >  2:275071/1‑150 (MQ=255)
                                                                             aaaaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAg                                                                          <  2:224967/150‑1 (MQ=255)
                                                                              aaaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                         <  1:962261/150‑1 (MQ=255)
                                                                              aaaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                         >  2:467925/1‑150 (MQ=255)
                                                                               aaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCCCTTGTTATCGAACGTCCAGCCATGAATCAGAt                                                                        >  1:1151236/1‑150 (MQ=255)
                                                                                  atgatgTAGACCTGATACATCGCTCAACCGGCTATTTTGCGTTTAGCATGAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACt                                                                     <  2:1218341/150‑1 (MQ=255)
                                                                                  atgatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTAATAACGCACCAAGAACGGGCGCATGTTATCCAAAGACAAGACAGGAATCATATTCt                                                                     >  1:1618157/1‑150 (MQ=255)
                                                                                    gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGc                                                                   <  2:28858/150‑1 (MQ=255)
                                                                                    gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGc                                                                   <  1:1368309/150‑1 (MQ=255)
                                                                                       gTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCAGTCTTACCGCCCCATGCACGGGCGCTTGTTATCGAACGGCCCGCCAGGAATCAGATACTGCAtt                                                                >  2:1719455/1‑150 (MQ=255)
                                                                                        tAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTc                                                               >  1:980085/1‑150 (MQ=255)
                                                                                         aGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATTGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTcc                                                              >  2:586963/1‑150 (MQ=255)
                                                                                         aGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTcc                                                              >  2:1484036/1‑150 (MQ=255)
                                                                                           aCCTGATACATAGCACAAAAGGATATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCca                                                            <  2:59268/150‑1 (MQ=255)
                                                                                             cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatc                                                          >  1:474909/1‑150 (MQ=255)
                                                                                               gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGTATTCCcatcgc                                                        <  1:530012/150‑1 (MQ=255)
                                                                                                     tAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGCTAGAGCCGGTTTTGGTTTTCTGTCTTAACACACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTc                                                  >  1:733122/1‑150 (MQ=255)
                                                                                                       gCACAACAGGCTATTTTGCGTTTAGCATCAGACTCAAACCTGCTCCAGATAGAGCCGGTTTTGGTTATCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGTAATCAGATACTGCATTCCCATCGCATCGTCAc                                                <  1:283511/150‑1 (MQ=255)
                                                                                                          caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCCGCCCGGAATCGGTAACTGCCTTCCCATCGCATCGTcgagcg                                             >  2:1110186/1‑145 (MQ=255)
                                                                                                            acaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCTCCATGCACGGGCGCTTGTTATCGACCGTCCCGCCAGGCATCAGATACTGCATTCCCATCGCATCGTCAcgcgag                                           >  1:273444/1‑148 (MQ=255)
                                                                                                              aGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGcc                                         <  2:684265/150‑1 (MQ=255)
                                                                                                               ggCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCa                                        <  1:1376143/150‑1 (MQ=255)
                                                                                                                  tATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCCGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                     >  2:1562950/1‑150 (MQ=255)
                                                                                                                   aTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                    <  1:1712183/150‑1 (MQ=255)
                                                                                                                   aTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                    <  1:830325/150‑1 (MQ=255)
                                                                                                                   aTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                    <  1:827605/150‑1 (MQ=255)
                                                                                                                    ttttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGccc                                   <  2:428547/150‑1 (MQ=255)
                                                                                                                         cGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGtt                              >  1:378551/1‑150 (MQ=255)
                                                                                                                            ttAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCCCGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                           >  2:984049/1‑150 (MQ=255)
                                                                                                                             tAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGGTTTGGTTTTGTGTCTTAACGCACCATGCACGGGCGCGTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCACATCGCAGCGTCCCGCGCGCCAAGCCCGTGTTTCTg                          >  1:1439468/1‑150 (MQ=255)
                                                                                                                               gCATCAGTCTCAAAGCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                        <  2:246345/150‑1 (MQ=255)
                                                                                                                               gCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                        <  1:382043/150‑1 (MQ=255)
                                                                                                                               gCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                        <  2:764250/150‑1 (MQ=255)
                                                                                                                                     aaccaaaaCTGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCt                  >  2:964072/6‑150 (MQ=255)
                                                                                                                                       ctcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                >  1:158664/1‑150 (MQ=255)
                                                                                                                                       ctcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAAGTCCAGCCAGGAATCAGATACTGCAATACCAGCGCATCGTCACGCGCGCCAAGGCCGTGGTTCTGATACAGCtaa                >  1:1355555/1‑148 (MQ=255)
                                                                                                                                        tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCATATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt               <  1:76044/150‑1 (MQ=255)
                                                                                                                                        tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt               >  2:575219/1‑150 (MQ=255)
                                                                                                                                         cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCATGCCCGTGTTTCTGATACAGCTCagg              >  2:69630/1‑148 (MQ=255)
                                                                                                                                          aaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGCTACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGAGTTTCTGATACAGCTCATgg             >  2:91807/1‑150 (MQ=255)
                                                                                                                                          aaaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATgg             <  2:450374/150‑1 (MQ=255)
                                                                                                                                           aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTTTCGAACGTCCCGCCAGGAATCAGATACTGCAGTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATggg            >  1:999280/1‑150 (MQ=255)
                                                                                                                                           aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATggg            >  2:448248/1‑150 (MQ=255)
                                                                                                                                            aCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGAACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGCGACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGc           >  2:1720613/1‑150 (MQ=255)
                                                                                                                                                 cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCCCCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca      <  1:1233524/150‑1 (MQ=255)
                                                                                                                                                 cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca      <  1:630602/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAC  >  NC_000913/2917870‑2918167

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: ‑