Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 653,723 G→A 5.8% A497T (GCT→ACT)  dpiB → sensory histidine kinase DpiB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913653,7230GA 5.8% 273.1 / 4.3 86A497T (GCT→ACT) dpiBsensory histidine kinase DpiB
Reads supporting (aligned to +/- strand):  ref base G (40/40);  new base A (2/3);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

TGCGTAGCGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAACCTAATGACAGCTCC  >  NC_000913/653580‑653872
                                                                                                                                               |                                                                                                                                                     
tGCGTAGCGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACg                                                                                                                                                 >  2:635604/1‑150 (MQ=255)
  cGTAGCGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAg                                                                                                                                               <  2:629405/150‑1 (MQ=255)
  cGTAGCGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGCACGCGTACTGAcccc                                                                                                                                               <  1:1040135/150‑4 (MQ=255)
     aGCGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGccc                                                                                                                                            <  2:642927/150‑1 (MQ=255)
       cGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCgg                                                                                                                                          >  1:704262/1‑150 (MQ=255)
       cGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGCACGCGTACTGAccccct                                                                                                                                             >  2:1172456/1‑142 (MQ=255)
       cGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGCACGCGTACTGAccccct                                                                                                                                             <  1:1172456/147‑6 (MQ=255)
       cGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGGCAGTACGCGTGCTGACGAGCCCgg                                                                                                                                          <  1:913776/150‑1 (MQ=255)
        gATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGt                                                                                                                                         >  1:1004764/1‑150 (MQ=255)
         aTGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTACAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTg                                                                                                                                        >  1:1375378/1‑150 (MQ=255)
          tGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGa                                                                                                                                       >  1:1426054/1‑150 (MQ=255)
           gAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGaa                                                                                                                                      >  2:1083405/1‑150 (MQ=255)
           gAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGaa                                                                                                                                      >  2:1088514/1‑150 (MQ=255)
            aaGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAAc                                                                                                                                     <  2:314180/150‑1 (MQ=255)
            aaGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAAc                                                                                                                                     <  2:203320/150‑1 (MQ=255)
                   cAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCAt                                                                                                                              >  2:1263318/1‑150 (MQ=255)
                     aGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTg                                                                                                                            <  2:1078666/150‑1 (MQ=255)
                        tCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGt                                                                                                                         <  1:364442/150‑1 (MQ=255)
                          gTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTg                                                                                                                       <  2:663701/150‑1 (MQ=255)
                          gTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTg                                                                                                                       <  1:583182/150‑1 (MQ=255)
                          gTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTg                                                                                                                       >  1:409114/1‑150 (MQ=255)
                            tGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTgg                                                                                                                           <  1:862158/144‑1 (MQ=255)
                            tGAATTATTCCTCAGAGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTgg                                                                                                                           >  2:862158/1‑144 (MQ=255)
                                ttattCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTg                                                                                                                 >  1:393659/1‑150 (MQ=255)
                                  attCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGAt                                                                                                               >  2:1306232/1‑150 (MQ=255)
                                  attCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGAt                                                                                                               <  1:1376930/150‑1 (MQ=255)
                                     ccTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGc                                                                                                            >  1:729839/1‑150 (MQ=255)
                                     ccTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTACTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGc                                                                                                            >  1:918413/1‑150 (MQ=255)
                                     ccTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTACAGAGTCACTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGACCGGTGAACATGGCATAGGGTTGTACTTTATTGc                                                                                                            >  2:1376264/1‑150 (MQ=255)
                                      cTCAGCTATGAATGCGATGTTGTTTTGATTGAAGTCGCCGTTCTGGGGTGGGGTGTTCCTGATTGTCTACGAGACAAAATATTTGAGCATGGTGTTAGTACGCGTTCTTACGAGCCCGGTGAATATGGGATTGGGTTTTACTTGATTTcc                                                                                                           <  2:1260785/150‑1 (MQ=255)
                                      cTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGcc                                                                                                           <  2:1998/150‑1 (MQ=255)
                                       tCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCa                                                                                                          <  2:1426054/150‑1 (MQ=255)
                                        cAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAg                                                                                                         >  1:205681/1‑150 (MQ=255)
                                           cGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTa                                                                                                      <  1:293050/150‑1 (MQ=255)
                                                aaGGCGATGATGTGGGGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTaa                                                                                                 <  1:1340300/150‑1 (MQ=255)
                                                  ggCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGa                                                                                                                                                >  1:714772/1‑101 (MQ=255)
                                                  ggCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGa                                                                                                                                                >  1:719614/1‑101 (MQ=255)
                                                  ggCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGa                                                                                                                                                <  2:714772/101‑1 (MQ=255)
                                                  ggCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGa                                                                                                                                                <  2:719614/101‑1 (MQ=255)
                                                   gCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACg                                                                                                                                                 >  1:1061725/1‑99 (MQ=255)
                                                   gCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACg                                                                                                                                                 <  2:1061725/99‑1 (MQ=255)
                                                        gtgggggtgATTGAAGTCGCGGATCAGGGCTGCGGCGTTCCAGAGTCTCTTCGATACAAAATTTTTGAGGAGGGGGTCAGTTCGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGc                                                                                         <  2:1186981/145‑1 (MQ=255)
                                                        gatgTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTg                                                                                                                                        <  2:1203263/103‑1 (MQ=255)
                                                        gatgTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTg                                                                                                                                        >  1:1203263/1‑103 (MQ=255)
                                                           gtggtgATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGTAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCggt                                                                                      >  2:113835/1‑150 (MQ=255)
                                                           gtggtgATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCggt                                                                                      >  1:948170/1‑150 (MQ=255)
                                                           gtggtgATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCggt                                                                                      >  2:99613/1‑150 (MQ=255)
                                                           gtggtgATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCggt                                                                                      >  1:1184798/1‑150 (MQ=255)
                                                                  ttGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTa                                                                               >  2:1275165/1‑150 (MQ=255)
                                                                   tGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTAt                                                                              >  1:378710/1‑150 (MQ=255)
                                                                    gAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACTATACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATc                                                                             <  2:930333/150‑1 (MQ=255)
                                                                    gAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATc                                                                             <  2:205681/150‑1 (MQ=255)
                                                                     aaGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCa                                                                            >  2:918201/1‑150 (MQ=255)
                                                                      aGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCTGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCAc                                                                           >  1:647973/1‑150 (MQ=255)
                                                                       gTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCAct                                                                          >  1:443094/1‑150 (MQ=255)
                                                                         cGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCActct                                                                        <  1:635604/150‑1 (MQ=255)
                                                                          gCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACtctc                                                                       <  1:847076/150‑1 (MQ=255)
                                                                                   ggCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAAt                                                              <  2:1004764/150‑1 (MQ=255)
                                                                                    gCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATg                                                             >  1:41767/1‑150 (MQ=255)
                                                                                    gCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATg                                                             >  1:90473/1‑150 (MQ=255)
                                                                                         ggcgTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATccc                                                        <  2:948170/150‑1 (MQ=255)
                                                                                            gTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGc                                                     <  1:1306232/150‑1 (MQ=255)
                                                                                            gTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGc                                                     >  1:421177/1‑150 (MQ=255)
                                                                                              tCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATcc                                                         <  2:906915/144‑1 (MQ=255)
                                                                                              tCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATcc                                                         >  1:906915/1‑144 (MQ=255)
                                                                                               ccAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGaa                                                                    <  1:1328402/132‑1 (MQ=255)
                                                                                               ccAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGaa                                                                    >  2:1328402/1‑132 (MQ=255)
                                                                                                cAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTa                                                 >  2:186406/1‑150 (MQ=255)
                                                                                                   agTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTTGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCt                                              <  2:1375378/150‑1 (MQ=255)
                                                                                                    gTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCtt                                             <  1:315719/150‑1 (MQ=255)
                                                                                                    gTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCtt                                             >  1:12691/1‑150 (MQ=255)
                                                                                                         tACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGATCATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTAtttt                                        >  1:248640/1‑150 (MQ=255)
                                                                                                                aaaaaTATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCta                                 <  1:1403448/149‑1 (MQ=255)
                                                                                                                  aaaTATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCtata                               >  1:1307344/1‑150 (MQ=255)
                                                                                                                        ttGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGa                         >  1:658987/1‑150 (MQ=255)
                                                                                                                         tGAGCAGGGGGTCAGTACGCGTACTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGaa                        <  2:918413/150‑1 (MQ=255)
                                                                                                                          gAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGaaa                       <  2:253729/150‑1 (MQ=255)
                                                                                                                             cAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTg                    <  2:729839/150‑1 (MQ=255)
                                                                                                                             cAGGGGGGCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTg                    <  2:919715/150‑1 (MQ=255)
                                                                                                                                 gggTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTTGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAAc                <  1:1088514/150‑1 (MQ=255)
                                                                                                                                 gggTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAAc                <  2:90473/150‑1 (MQ=255)
                                                                                                                                 gggTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAAc                <  1:1083405/150‑1 (MQ=255)
                                                                                                                                 gggTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAAc                <  2:41767/150‑1 (MQ=255)
                                                                                                                                    tCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAACCTa             >  1:1165279/1‑150 (MQ=255)
                                                                                                                                      aGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAACCTAAt           <  1:1275165/150‑1 (MQ=255)
                                                                                                                                         aCGCGTGCTGCCGCGCCCGGTGAACATGGCTTTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTGTCCATCTATATTCCGAAAGTGAAACCTAATGAc        >  1:1081360/1‑150 (MQ=255)
                                                                                                                                               gCTGACGAGCCCGGTGAACATGGCATTGGGGTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCCCTCTCGAAGATAATGATCCCTGCTGTACCTTATTTTCAATCTATATTCCGAAAGTGCAACCTAATGACAGCTcc  >  2:796971/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TGCGTAGCGATGAAGGAAACAAGATCGTTGAATTATTCCTCAGCGATGAAGGCGATGATGTGGTGATTGAAGTCGCCGATCAGGGCTGCGGCGTTCCAGAGTCTCTACGAGACAAAATATTTGAGCAGGGGGTCAGTACGCGTGCTGACGAGCCCGGTGAACATGGCATTGGGTTGTACTTGATTGCCAGCTACGTAACGCGCTGCGGTGGTGTTATCACTCTCGAAGATAATGATCCCTGCGGTACCTTATTTTCAATCTATATTCCGAAAGTGAAACCTAATGACAGCTCC  >  NC_000913/653580‑653872

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑