Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,017 C→T 10.6% intergenic (+204/+28) 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,0170CT10.6% 370.5 / 24.1 134intergenic (+204/+28)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base C (66/52);  new base T (4/10);  total (72/62)
Fisher's exact test for biased strand distribution p-value = 8.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATC  >  NC_000913/2917868‑2918165
                                                                                                                                                     |                                                                                                                                                    
gACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTcc                                                                                                                                                      <  2:1021704/150‑1 (MQ=255)
  cAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACACCAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAg                                                                                                                                                    >  2:689925/1‑150 (MQ=255)
       tGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATaga                                                                                                                                               >  1:1419112/1‑150 (MQ=255)
         ccTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGc                                                                                                                                             <  1:752917/150‑1 (MQ=255)
         ccTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGCAAATTTTGAGTTTAGAATCAGTATAAAACCGGCTCCAGATAGAGc                                                                                                                                             >  2:292117/1‑150 (MQ=38)
          cTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGcc                                                                                                                                            >  2:875456/1‑150 (MQ=255)
           tgtgTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCg                                                                                                                                           >  1:531162/1‑150 (MQ=255)
              ttttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtt                                                                                                                                        <  1:1661270/149‑1 (MQ=255)
                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                      <  1:664293/150‑2 (MQ=255)
                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                      <  1:1467764/150‑2 (MQ=255)
                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                      <  1:946187/150‑2 (MQ=255)
                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAAATAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttg                                                                                                                                      <  2:1197355/150‑2 (MQ=38)
                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttgag                                                                                                                                    <  1:1600395/150‑4 (MQ=38)
                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                                    <  1:1665144/150‑1 (MQ=255)
                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                                    <  1:660552/150‑1 (MQ=255)
                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                                    <  1:1477958/150‑1 (MQ=255)
                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttgaga                                                                                                                                   <  1:30055/150‑5 (MQ=38)
                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggt                                                                                                                                   >  2:1417091/1‑150 (MQ=255)
                   ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGAATCAGTCTCAAACCGGCTCTATCTGGAGCCggttttggt                                                                                                                                   >  1:1549449/2‑150 (MQ=255)
                     cctgATGCGGCGTTAACGCCTTATACGGCCTACAATACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTTTTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                 <  1:851301/147‑1 (MQ=255)
                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttgagact                                                                                                                                 <  1:1213814/150‑7 (MQ=38)
                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                 <  2:1493684/150‑1 (MQ=255)
                      cggttaaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGCCCTGATACATAGCACAACAGGCTATTTTGCGGTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                >  1:1244643/8‑150 (MQ=255)
                      cGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                >  1:984527/1‑150 (MQ=255)
                      cGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCACTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                >  2:1597961/1‑150 (MQ=255)
                       ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                               <  1:1025606/150‑1 (MQ=255)
                       ggATGCGGCGTTAACGCCTTACCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                               <  1:358495/150‑1 (MQ=255)
                             ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCtt                                                                                                                         >  2:647608/1‑150 (MQ=255)
                              gcgTTAACGCCTTATCCGGCCTACAAAACCAATTTAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACATGCTATTTTGCGTTTATCATCAGTCTCAAACCGGCACGATAAAGAGCCGGTTTTGGTTTTGTGTCttt                                                                                                                        >  2:1083403/1‑149 (MQ=38)
                               cgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACATGATACATAGAACAACAGGCTATTTTGCGTTTAGAATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGTTTTTATGTCTTaa                                                                                                                       >  2:1045443/1‑150 (MQ=255)
                                   aaCGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                   >  2:42384/1‑150 (MQ=255)
                                   aaCGCCTTATCCGGCCTACAAAACAAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACACCAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTTTCTGTCTTAACgaa                                                                                                                   >  2:273373/1‑148 (MQ=255)
                                     cGCCTTATCCGGCCTACAAAACCTATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAcc                                                                                                                 >  2:171327/1‑150 (MQ=255)
                                      gcattATCCGGCCTACAAAACCAATTAAATTGAATGAATTAAAAATTGATGTAGAACTGATACATAGCACAACAGGCTATTTTGTGTTTAGCATCAGTCTCAAACCGGCTCCCGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCa                                                                                                                <  2:855005/147‑1 (MQ=255)
                                          tATGCGGCCTACAATACCAATTAAATTCAATGAATTAAAAAATGTTGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCa                                                                                                            <  2:190792/150‑1 (MQ=255)
                                          tATCAGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCa                                                                                                            <  1:248965/150‑1 (MQ=255)
                                            tCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACg                                                                                                          >  1:1028526/1‑150 (MQ=255)
                                                  cTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGAGCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCGGGTTTTGGTTTTCTGTCTTAACGCACCCTGCAcgcgcgat                                                                                                    >  2:93597/1‑148 (MQ=255)
                                                  cTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCt                                                                                                    <  1:1388382/150‑1 (MQ=255)
                                                      aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGtt                                                                                                >  1:1292957/1‑150 (MQ=255)
                                                      aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGtt                                                                                                >  1:1344895/1‑150 (MQ=255)
                                                       aaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                               >  2:1521379/1‑150 (MQ=255)
                                                          caaaTTAAATTCAATTAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCg                                                                                            <  2:666264/148‑1 (MQ=255)
                                                           cAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTGTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGa                                                                                           >  2:242292/1‑150 (MQ=255)
                                                             aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCTGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                         >  1:208017/1‑150 (MQ=255)
                                                              ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACg                                                                                        <  2:234583/150‑1 (MQ=255)
                                                              ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACg                                                                                        <  1:1694182/150‑1 (MQ=255)
                                                              ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACg                                                                                        <  2:766184/150‑1 (MQ=255)
                                                                 aaTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACTTcc                                                                                     <  2:873024/150‑1 (MQ=255)
                                                                  aTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcca                                                                                    >  2:1330019/1‑150 (MQ=255)
                                                                     cAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                 <  1:1411718/150‑1 (MQ=255)
                                                                       aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGGTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                               <  2:1594899/150‑1 (MQ=255)
                                                                         gAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGCCGGTTTTGTTTTTCTGTATTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCAAGGa                                                                             >  1:1280628/1‑150 (MQ=255)
                                                                         gAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGa                                                                             >  2:626845/1‑150 (MQ=255)
                                                                          aaTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGat                                                                            >  1:1217169/1‑149 (MQ=255)
                                                                          aaTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaa                                                                            <  2:862090/150‑1 (MQ=255)
                                                                            ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAta                                                                          >  1:154177/1‑149 (MQ=255)
                                                                            ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                          >  2:1086435/1‑150 (MQ=255)
                                                                             tAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCa                                                                         <  1:292117/150‑1 (MQ=255)
                                                                                 aaaTGATGTAGACCTGATACATAGCACAACACTCTATTTTGCGTTTAGCATCCGTCTCACACCGGCTCCAGATAGCGCCGGGTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATa                                                                     <  1:853588/150‑1 (MQ=255)
                                                                                   atgatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACt                                                                   <  1:42384/150‑1 (MQ=255)
                                                                                   atgatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACt                                                                   >  1:595846/1‑150 (MQ=255)
                                                                                     gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGGCTCACAGCGGCTCCAGATAGATCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGc                                                                 <  1:277732/150‑1 (MQ=255)
                                                                                      atgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCa                                                                >  1:1448332/1‑150 (MQ=255)
                                                                                         tAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGAACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTc                                                             >  1:851934/1‑150 (MQ=255)
                                                                                           gACCTGATACATAGCACAACAGGCTATTTTGCGATTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCCCCATGCACGGGCGCTTGTTATAGAACGTCCAGCCAGGAATCAGATACTGCATAccc                                                           >  2:1660693/1‑150 (MQ=255)
                                                                                              cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatc                                                        <  1:1381203/150‑1 (MQ=255)
                                                                                                gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCcatcgc                                                      >  2:1204080/1‑150 (MQ=255)
                                                                                                  tACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGCTACTGCATTCCcatcgccg                                                    >  2:1635810/1‑148 (MQ=255)
                                                                                                  aaCATAGCACAACAGGCTATTTTGTGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcat                                                    <  1:1178943/149‑1 (MQ=255)
                                                                                                    cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatcg                                                  <  1:1086435/150‑1 (MQ=255)
                                                                                                    cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatcg                                                  >  1:1114592/1‑150 (MQ=255)
                                                                                                    cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGGTTTGGTTTTCTGTCTTAACGCACCATTCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCcatcgcatag                                                  >  2:734019/1‑148 (MQ=255)
                                                                                                     aTAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGt                                                 >  1:342811/1‑150 (MQ=255)
                                                                                                     aTAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTACAGCCAGGAATCAGATACTGCATTCCCATCGCATCGt                                                 >  1:1443390/1‑150 (MQ=255)
                                                                                                       aGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCa                                               <  1:176949/150‑1 (MQ=255)
                                                                                                       aGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCa                                               <  1:171327/150‑1 (MQ=255)
                                                                                                        gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAc                                              <  2:572869/150‑1 (MQ=255)
                                                                                                           caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATAGCATCGTCAcgcg                                           >  2:783172/1‑150 (MQ=255)
                                                                                                                      tttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCg                                <  1:1324687/150‑1 (MQ=255)
                                                                                                                            tttAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGGTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                          >  1:770631/1‑150 (MQ=255)
                                                                                                                             ttAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                         >  2:1533950/1‑150 (MQ=255)
                                                                                                                               aGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                       <  1:1381394/150‑1 (MQ=255)
                                                                                                                               aGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGGTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTGCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGAGCGCCACGCCCGTGTGTCTgg                       >  1:329518/1‑149 (MQ=255)
                                                                                                                                gCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                      <  1:1136712/150‑1 (MQ=255)
                                                                                                                                gCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                      <  2:1712823/150‑1 (MQ=255)
                                                                                                                                   tCAGTCTCAAACCGGCTCCAGATAGAGCCGGGTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCAGTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACa                   >  1:1168870/1‑150 (MQ=255)
                                                                                                                                    cAGTCTCAAACCCGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAg                  <  2:1549449/150‑1 (MQ=255)
                                                                                                                                       accaaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATTCTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTc               >  1:1586409/5‑150 (MQ=255)
                                                                                                                                         tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt             <  2:748673/150‑1 (MQ=255)
                                                                                                                                         caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATAGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt             >  2:719100/3‑150 (MQ=255)
                                                                                                                                          cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCTCGTGTTTCTGATACAGCTCATg            >  2:970585/1‑150 (MQ=255)
                                                                                                                                            aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTAATggg          >  2:609766/1‑150 (MQ=255)
                                                                                                                                                ggCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCCAGCCCGTGTTTCTGATCCAGCTCATGGGCttt      >  2:1077486/1‑150 (MQ=255)
                                                                                                                                                  cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAGAGCTCATGGGCTTtca    <  1:1606004/150‑1 (MQ=255)
                                                                                                                                                    ccAGATAGAGCCGTTTTTTGTTTTTTTTCTTATCGCTCCATGAACGGGAGCTTGTTATCTAACGTCCATAAAGGAATCAGATTCTGAATTCTAATCGGATAGTCACGCGCGCGAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcatc  <  2:308505/150‑1 (MQ=255)
                                                                                                                                                    ccAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGTATCAGATACTGCATTCCCATCGCATCCTCACTAGCGCCAAGCCCGTGTTTGTGATACCGCTCATGGGCttagatc  >  1:334412/1‑145 (MQ=255)
                                                                                                                                                    ccAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCCCGCGCGCCAAGCCCGTGTTGCTGATACAGCTCATGGGCTTtcatc  >  2:151177/1‑150 (MQ=255)
                                                                                                                                                    ccAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAc                    >  2:339072/1‑132 (MQ=255)
                                                                                                                                                    ccAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcatc  >  2:1617109/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATC  >  NC_000913/2917868‑2918165

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: