Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,765,162 C→T 5.4% intergenic (‑40/+60) ydiH ← / ← menI protein YdiH/1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,765,1620CT 5.4% 292.7 / 3.5 92intergenic (‑40/+60)ydiH/menIprotein YdiH/1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base C (30/56);  new base T (3/2);  total (34/58)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAAC  >  NC_000913/1765014‑1765310
                                                                                                                                                    |                                                                                                                                                    
cAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATaccacc                                                                                                                                                     >  1:134787/1‑150 (MQ=255)
  aGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGt                                                                                                                                                   >  2:664913/1‑150 (MQ=255)
       tAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACt                                                                                                                                              <  2:1119877/150‑1 (MQ=255)
           tgttTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGc                                                                                                                                          >  2:1255073/1‑150 (MQ=255)
            gttTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCa                                                                                                                                         <  1:451354/150‑1 (MQ=255)
            gttTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCa                                                                                                                                         >  2:941101/1‑150 (MQ=255)
             tttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAg                                                                                                                                        <  2:254004/150‑1 (MQ=255)
             tttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAg                                                                                                                                        <  1:941101/150‑1 (MQ=255)
              ttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGa                                                                                                                                       >  2:1037543/1‑150 (MQ=255)
              ttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAgg                                                                                                                                       <  1:848044/150‑2 (MQ=255)
               tCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTTGGTGGTATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGAt                                                                                                                                      <  2:1025487/150‑1 (MQ=255)
                cAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGTTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATa                                                                                                                                     <  2:861164/150‑1 (MQ=255)
                cAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATa                                                                                                                                     >  2:932555/1‑150 (MQ=255)
                 cgACGTTTTAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGTCGTAAAAATTCAATGGCCAGCTGGTTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATaa                                                                                                                                    <  2:733401/149‑1 (MQ=255)
                 aGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATaa                                                                                                                                    <  1:299353/150‑1 (MQ=255)
                  gACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAAt                                                                                                                                   <  2:1505769/150‑1 (MQ=255)
                    cGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGa                                                                                                                                 >  1:160812/1‑150 (MQ=255)
                      tttGAGATATTGGGCAGGTGAAAGATTGCTCTTGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGAcc                                                                                                                               <  2:1252835/150‑1 (MQ=255)
                      tttGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCGAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGAcc                                                                                                                               >  1:407589/1‑150 (MQ=255)
                        tGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATAACACCGTTTACTGTGCAGATAATGACCaa                                                                                                                             <  2:430660/150‑1 (MQ=255)
                           gaTATTGGGCAGGGGAAAGATTGCGCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTGGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAg                                                                                                                          <  2:1076058/150‑1 (MQ=255)
                            atatTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGc                                                                                                                         <  1:448851/150‑1 (MQ=255)
                             tatTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCa                                                                                                                        <  2:1482747/150‑1 (MQ=255)
                              atTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCaa                                                                                                                       <  2:92483/150‑1 (MQ=255)
                                tGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAata                                                                                                                     <  2:134787/150‑1 (MQ=255)
                                      ggTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTc                                                                                                               <  2:967289/150‑1 (MQ=255)
                                       gTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTca                                                                                                              >  1:334952/1‑150 (MQ=255)
                                         gAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGTTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTcaca                                                                                                            >  2:1039900/1‑150 (MQ=255)
                                          aaaGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCacac                                                                                                           >  1:448264/1‑150 (MQ=255)
                                           aaGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACt                                                                                                          >  1:1250112/1‑150 (MQ=255)
                                            aGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACtt                                                                                                         >  1:777508/1‑150 (MQ=255)
                                                 gCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTg                                                                                                    >  2:710327/1‑150 (MQ=255)
                                                  ctctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTgg                                                                                                   >  1:610641/1‑150 (MQ=255)
                                                    ttGGTTACGGCGTAAAAATTCAATGGCAAGCTGGGTGGGATAAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTGACACTTTTCTGGTg                                                                                                 <  2:448264/149‑1 (MQ=255)
                                                    ctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTg                                                                                                 <  2:386700/150‑1 (MQ=255)
                                                    ctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTg                                                                                                 <  2:466201/150‑1 (MQ=255)
                                                    ctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTg                                                                                                 <  2:627005/150‑1 (MQ=255)
                                                    ctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGGTGAGTAGACATAGCATCCTCTCTTTTTGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTg                                                                                                 <  2:1513159/150‑1 (MQ=255)
                                                    ctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGAGCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTg                                                                                                 <  2:364443/150‑1 (MQ=255)
                                                     tGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGa                                                                                                >  1:1123920/1‑150 (MQ=255)
                                                       gTCTCGGCTTAAAAATTCAATGGCGAGGTGGGTGGGATCATGTTTAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTGACACTTTTCTGGTGACa                                                                                              <  1:252500/150‑1 (MQ=255)
                                                       gTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACa                                                                                              <  2:275234/150‑1 (MQ=255)
                                                       gTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACa                                                                                              <  1:51473/150‑1 (MQ=255)
                                                       gTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACa                                                                                              >  2:920924/1‑150 (MQ=255)
                                                              cGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTcac                                                                                       >  2:1388024/1‑150 (MQ=255)
                                                              cGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTcac                                                                                       <  1:664913/150‑1 (MQ=255)
                                                               gTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTcaca                                                                                      <  2:191089/150‑1 (MQ=255)
                                                                     aTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATgg                                                                                <  2:1003804/150‑1 (MQ=255)
                                                                     aTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATgg                                                                                <  2:949327/150‑1 (MQ=255)
                                                                      ttCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGc                                                                               <  1:692841/150‑1 (MQ=255)
                                                                      ttCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGc                                                                               <  2:854480/150‑1 (MQ=255)
                                                                      ttCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGc                                                                               <  1:1431193/150‑1 (MQ=255)
                                                                       tCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCg                                                                              <  2:373204/150‑1 (MQ=255)
                                                                       tCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCg                                                                              <  2:1004801/150‑1 (MQ=255)
                                                                           tGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGgtcg                                                                          <  2:610641/150‑1 (MQ=255)
                                                                            ggCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTATTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGAAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGAGgtagt                                                                         >  1:730578/1‑150 (MQ=255)
                                                                                      gtgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGAc                                                               <  1:783742/150‑1 (MQ=255)
                                                                                       tgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACg                                                              <  2:1123920/150‑1 (MQ=255)
                                                                                         ggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGaa                                                            <  1:433204/150‑1 (MQ=255)
                                                                                          gATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAAc                                                           >  2:292077/1‑150 (MQ=255)
                                                                                           aTCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACa                                                          <  2:1222883/150‑1 (MQ=255)
                                                                                             gtaGTTGTGTTGACATAGCATCCTCGCTTTTAGACTTGACCTGCACCGTATACCACCGTTTAGTGTGCAGATAATGACCACAAGCAATATTCGTCACACTTTTGTGGTGACAAGGTCACAAAATGGCGGTCGTCAATCGTGTGGAACAGc                                                        <  2:542069/148‑1 (MQ=255)
                                                                                               aGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGcac                                                      >  2:707485/1‑150 (MQ=255)
                                                                                                         aCATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCtt                                            <  2:891802/150‑1 (MQ=255)
                                                                                                           aTAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTtc                                          >  1:1329740/1‑150 (MQ=255)
                                                                                                            tAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTtct                                         <  2:1503406/150‑1 (MQ=255)
                                                                                                             aGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTtctc                                        >  2:203133/1‑150 (MQ=255)
                                                                                                              gCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCa                                       <  2:298323/150‑1 (MQ=255)
                                                                                                               cATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCAt                                      <  1:932555/150‑1 (MQ=255)
                                                                                                                aTCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATc                                     >  2:1465960/1‑150 (MQ=255)
                                                                                                                   cTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGaa                                  <  2:493429/150‑1 (MQ=255)
                                                                                                                    tCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAg                                 <  1:581018/150‑1 (MQ=255)
                                                                                                                    tCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAg                                 <  2:334952/150‑1 (MQ=255)
                                                                                                                       cTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGAtt                              >  2:595040/1‑150 (MQ=255)
                                                                                                                        ttttAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGAttt                             >  2:1032682/1‑150 (MQ=255)
                                                                                                                        ttttAGACAAGACCTGCACAGTATACCAACGTTTACTGTGCAGAGAATGACCAAAACCAATATGCGTCCCAGTTTTCTGGTGACCACGTCACAAAATGGCGGTCGTCAAGCGTCACGCACAGCACAAACGCCCAATCTCATCGAAGAttt                             >  2:1323752/1‑150 (MQ=255)
                                                                                                                         tttAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTc                            <  1:710327/150‑1 (MQ=255)
                                                                                                                         tttAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTc                            <  1:595040/150‑1 (MQ=255)
                                                                                                                          ttAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCa                           <  1:1079277/150‑1 (MQ=255)
                                                                                                                           tAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGACCAGCACAAACGCCCTTTCTCATCGAAGATTTCaa                          >  2:624335/1‑150 (MQ=255)
                                                                                                                             gACAGGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATc                        >  2:884027/1‑150 (MQ=255)
                                                                                                                             gACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTAACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATc                        >  2:99802/1‑150 (MQ=255)
                                                                                                                                 aGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGcc                    <  1:460508/150‑1 (MQ=255)
                                                                                                                                  gACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCa                   >  2:1133411/1‑150 (MQ=255)
                                                                                                                                   atcTGCACAGTAAACGGTGGTATACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAg                  >  1:1272704/3‑150 (MQ=255)
                                                                                                                                     cTGCACAGTAAACGGTGGTATACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAc                >  1:201043/1‑150 (MQ=255)
                                                                                                                                        cacaGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTg             <  1:168409/150‑1 (MQ=255)
                                                                                                                                        cacaGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTg             <  1:1255073/150‑1 (MQ=255)
                                                                                                                                         acaGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTgg            <  1:89055/150‑1 (MQ=255)
                                                                                                                                         acaGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTgg            <  2:1250112/150‑1 (MQ=255)
                                                                                                                                              ataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAcc               >  2:629940/1‑142 (MQ=255)
                                                                                                                                              ataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAcc               <  1:629940/142‑1 (MQ=255)
                                                                                                                                               taCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcgc      >  1:1169634/1‑150 (MQ=255)
                                                                                                                                                accaccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcgcg     >  1:526147/1‑150 (MQ=255)
                                                                                                                                                  caccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGaa   <  2:31934/150‑1 (MQ=255)
                                                                                                                                                   accGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAAc  <  1:203133/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAAC  >  NC_000913/1765014‑1765310

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