Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,020 A→C 6.4% intergenic (+207/+25) 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,0200AC 6.4% 338.2 / 11.3 120intergenic (+207/+25)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base A (61/51);  new base C (3/5);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

Only 100 of 120 total aligned reads displayed.

CTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAC  >  NC_000913/2917878‑2918167
                                                                                                                                               |                                                                                                                                                   
cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGcc                                                                                                                                              <  2:1492978/150‑1 (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)
      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                        <  2:1781656/150‑2 (MQ=255)
      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                        <  1:230647/150‑2 (MQ=255)
      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                        >  1:1499337/1‑150 (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)
           ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                   <  2:881034/150‑1 (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)
                         aaCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                     <  1:901398/150‑1 (MQ=255)
                              cTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                                <  1:1326738/150‑1 (MQ=255)
                              cTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                                <  2:1398068/150‑1 (MQ=255)
                                tATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAACACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGATCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCa                                                                                                              >  2:543591/1‑150 (MQ=255)
                                  tCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACg                                                                                                            >  2:630602/1‑150 (MQ=255)
                                   ccGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                           >  2:1476940/1‑150 (MQ=255)
                                   ccGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                           <  1:1263876/150‑1 (MQ=255)
                                       ccTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                       >  2:794828/1‑150 (MQ=255)
                                        cTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCt                                                                                                      >  1:764250/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:20618/1‑150 (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)
                                                      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)
                                                           cAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                   >  1:682670/1‑150 (MQ=255)
                                                            aaTGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                  >  2:805058/1‑150 (MQ=255)
                                                             aTGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                 <  1:1114697/150‑1 (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                                                                         >  2:467925/1‑150 (MQ=255)
                                                                      aaaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                         <  1:962261/150‑1 (MQ=255)
                                                                       aaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCCCTTGTTATCGAACGTCCAGCCATGAATCAGAt                                                                        >  1:1151236/1‑150 (MQ=255)
                                                                          atgatgTAGACCTGATACATCGCTCAACCGGCTATTTTGCGTTTAGCATGAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACt                                                                     <  2:1218341/150‑1 (MQ=255)
                                                                            gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGc                                                                   <  1:1368309/150‑1 (MQ=255)
                                                                            gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGc                                                                   <  2:28858/150‑1 (MQ=255)
                                                                                tAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTc                                                               >  1:980085/1‑150 (MQ=255)
                                                                                 aGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATTGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTcc                                                              >  2:586963/1‑150 (MQ=255)
                                                                                 aGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTcc                                                              >  2:1484036/1‑150 (MQ=255)
                                                                                   aCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATCGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCca                                                            <  1:51476/150‑1 (MQ=255)
                                                                                   aCCTGATACATAGCACAAAAGGATATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCca                                                            <  2:59268/150‑1 (MQ=255)
                                                                                       gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGTATTCCcatcgc                                                        <  1:530012/150‑1 (MQ=255)
                                                                                       gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgc                                                        <  2:1673627/150‑1 (MQ=255)
                                                                                           cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatcg                                                    >  1:1369219/1‑150 (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)
                                                                                                       ggCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTTTCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCa                                        >  2:71638/1‑150 (MQ=255)
                                                                                                       ggCTATTTTGCGTTTAGCATCAGTCACAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCATAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCACAACGCATAGTCACGCGCGCCa                                        >  2:346487/1‑150 (MQ=255)
                                                                                                          tATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCCGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                     >  2:1562950/1‑150 (MQ=255)
                                                                                                           aTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                    <  1:827605/150‑1 (MQ=255)
                                                                                                           aTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                    <  1:1712183/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)
                                                                                                                       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)
                                                                                                                                tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt               <  1:324706/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)
                                                                                                                                    aCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGAACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGCGACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGc           >  2:1720613/1‑150 (MQ=255)
                                                                                                                                           ccAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTACAGCCAGGAATCAGATACTTCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcatc    >  1:1165307/1‑150 (MQ=255)
                                                                                                                                             aGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAc  >  1:442795/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAA‑TTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAC  >  NC_000913/2917878‑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: ‑