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.

CAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAC  >  NC_000913/2917870‑2918167
                                                                                                                                                      |                                                                                                                                                    
cAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAg                                                                                                                                                      >  1:1018917/1‑150 (MQ=255)
        cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGcc                                                                                                                                              <  1:802451/150‑1 (MQ=255)
           tttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggt                                                                                                                                           >  1:443609/8‑150 (MQ=39)
           tgtTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggt                                                                                                                                           <  1:1125865/150‑1 (MQ=38)
            ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACACCAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtt                                                                                                                                          >  2:121230/7‑150 (MQ=39)
            gtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtt                                                                                                                                          >  2:587266/1‑150 (MQ=255)
              ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                        <  1:230647/150‑2 (MQ=255)
               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttg                                                                                                                                       <  1:802703/150‑1 (MQ=255)
               tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggatttg                                                                                                                                       >  2:1663242/1‑150 (MQ=255)
                aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttgag                                                                                                                                      <  2:1492587/150‑4 (MQ=38)
                  tccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGCAGCCggtttgagac                                                                                                                                    <  1:1686080/149‑6 (MQ=38)
                  gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAAGGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtt                                                                                                                                    >  1:450374/1‑150 (MQ=255)
                  gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtt                                                                                                                                    >  1:428547/1‑150 (MQ=255)
                     ggataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGCCGGTTTTGGTTTTc                                                                                                                                 >  1:366452/7‑150 (MQ=255)
                      gATGCGGAGTTAACGCCTTATCCGGCCTACAAAAACAATTAAATTCAATGAATTTAAAAAATGATTTAGACCTGATACATAGGACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                 <  2:1539607/150‑1 (MQ=38)
                       aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCATTGAATTAAAAAATGATGTAGACCTGATACATAGCACAAGAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                               <  2:126199/150‑1 (MQ=255)
                       aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                               <  1:746967/150‑1 (MQ=255)
                       aTGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                               <  1:1717997/150‑1 (MQ=255)
                         tcggcgTTAACTCGTTATCCGGCCTACATAACCACTTAAATTCAATTAATTAAAAAATGTTGTAGACCTGATACATCGCACAACAGGCTATTTTGCGTCTAGCATCAGTCTCAACCCGGCTCCAGATTGTGACGGTTTTGGTTTTCTGTc                                                                                                                             <  1:950412/149‑1 (MQ=255)
                         gcggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTc                                                                                                                             <  2:176055/150‑1 (MQ=255)
                          cggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                            >  2:962261/1‑150 (MQ=255)
                            gcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTa                                                                                                                          >  2:1622191/1‑150 (MQ=255)
                               ttAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACg                                                                                                                       >  2:1733755/1‑150 (MQ=255)
                                 aaCGCCTTATCCGGCCTACAAACCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                     >  2:1143581/1‑150 (MQ=255)
                                 aaCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                     >  1:684265/1‑150 (MQ=255)
                                      cTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                                <  2:1398068/150‑1 (MQ=255)
                                      cTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                                <  1:1326738/150‑1 (MQ=255)
                                       ttATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTATCTGTCTTAACGCACCAtcc                                                                                                               >  2:866339/1‑148 (MQ=255)
                                        tATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAACACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGATCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCa                                                                                                              >  2:543591/1‑150 (MQ=255)
                                          tCCGGCCTACAAAACCAATTAACTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACg                                                                                                            <  2:1510666/150‑1 (MQ=255)
                                           ccGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                           <  1:1263876/150‑1 (MQ=255)
                                           ccGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                           >  2:1476940/1‑150 (MQ=255)
                                               ccTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                       >  2:794828/1‑150 (MQ=255)
                                                 tACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACAATGCACGGGCGCtt                                                                                                     >  2:78147/1‑150 (MQ=255)
                                                 tACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGCCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATTCACGGGCGCtt                                                                                                     >  2:1323562/1‑150 (MQ=255)
                                                  aCAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  1:587266/150‑1 (MQ=255)
                                                  aCAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  2:635113/150‑1 (MQ=255)
                                                  aCAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  1:1492480/150‑1 (MQ=255)
                                                  aCAAAAACAATTAAATTCACTGAATTAAAAAATGATGTAGACCGCCTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTg                                                                                                    <  2:1154745/150‑1 (MQ=255)
                                                    aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGtt                                                                                                  >  2:1398798/1‑150 (MQ=255)
                                                     aaaCCAATTAAATTCAATTAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGGCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                 <  2:248991/150‑1 (MQ=255)
                                                     aaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                                 <  1:625681/150‑1 (MQ=255)
                                                      aaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGTTACATAGCACAACAGGATATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCTGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCTCGGGCGCTTGTTAt                                                                                                <  1:826095/150‑1 (MQ=255)
                                                      aaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                <  2:1455475/150‑1 (MQ=255)
                                                          aaTTAAATTCAATGAATTAAAAAATGATGTAGCCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTGTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAGCGac                                                                                            >  2:709984/1‑149 (MQ=255)
                                                              aaaTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACACGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTc                                                                                        <  1:47906/150‑1 (MQ=255)
                                                                 ttCAATGAATTAAAAAATGATGTAGCCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccag                                                                                     >  1:1137850/1‑150 (MQ=255)
                                                                 ttCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGATCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTAGTTATCCAACGTccac                                                                                     >  2:1345519/1‑149 (MQ=255)
                                                                 ttCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccag                                                                                     <  1:440556/150‑1 (MQ=255)
                                                                    aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                  >  2:805058/1‑150 (MQ=255)
                                                                     aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                 <  1:1114697/150‑1 (MQ=255)
                                                                      tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTATTAACGAACCATGCACGGGCGATTGTTATGGAACGTCCAGCCAgg                                                                                >  1:239069/1‑150 (MQ=255)
                                                                        aaTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaa                                                                              >  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)
                                                                                                                                                      |                                                                                                                                                    
CAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAC  >  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: