Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,765,157 T→A 7.9% intergenic (‑35/+65) ydiH ← / ← menI uncharacterized protein/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,1570TA 7.9% 240.1 / 8.2 89intergenic (‑35/+65)ydiH/menIuncharacterized protein/1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base A (2/5);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

AAATTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGC  >  NC_000913/1765009‑1765306
                                                                                                                                                    |                                                                                                                                                     
aaaTTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGtata                                                                                                                                                      <  1:1640462/150‑1 (MQ=255)
 aaTTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATac                                                                                                                                                     <  1:414744/150‑1 (MQ=255)
  aTTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATacc                                                                                                                                                    <  2:168441/150‑1 (MQ=255)
  aTTCAAGTTCTAGTTGTTTCAGACGTTTGAGATAATGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATacc                                                                                                                                                    <  2:234398/150‑1 (MQ=255)
    tCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATaccac                                                                                                                                                  >  1:409845/1‑150 (MQ=255)
    acAATTTATAGTTGTTTCAGACGTATGAGATATTGGGCAGGTTAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTTGGATCAAGTTTAGTAGACATAGCATCCTCGCTATTAGACAAGACCTGCACAGTATaccac                                                                                                                                                  <  1:320504/149‑1 (MQ=255)
       aGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGt                                                                                                                                               >  1:685604/1‑150 (MQ=255)
       aGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCATCGCTTTTAGACAAGACCTGCACAGTATACCACCGt                                                                                                                                               >  1:696638/1‑150 (MQ=255)
             tgttgttTCAGACGTTTGATATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAATTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACtg                                                                                                                                         <  1:1464973/149‑1 (MQ=255)
              gttgttTCAGATTATTTAGATAATGTGCTGTGGAAATATTGTTCTTGTAACGGAGTAAAAATTCAATGGCAAGCTGGGTGTTATCAATTTGAGTATACATAGCATCCTCGATTTTAGACAAGACCTGCACAGTATACCACCGTTTACtgt                                                                                                                                        <  1:115272/150‑1 (MQ=255)
                tgttTCAGACGTTTGAGATTTTGGGCAGGTGATAGATTGCTCTGGTCACTGGGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATCCCACCGTTTACTGTGc                                                                                                                                      <  1:79002/150‑1 (MQ=255)
                tgttTCAGACGTTTGAGATTTTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGc                                                                                                                                      <  2:1424706/150‑1 (MQ=255)
                tggttCCGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGc                                                                                                                                      <  1:1545513/147‑1 (MQ=255)
                 gttTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCa                                                                                                                                     >  2:1662355/1‑150 (MQ=255)
                  tttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAg                                                                                                                                    >  2:1710364/1‑150 (MQ=255)
                  tttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAg                                                                                                                                    <  2:1562983/150‑1 (MQ=255)
                  ttgcAGACGTTTGAGATATTGGTCAGGTGAAAGATTGCTCTGGTCACGGAGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAg                                                                                                                                    <  2:75565/147‑1 (MQ=255)
                   ttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAgg                                                                                                                                   <  1:537476/150‑2 (MQ=255)
                   ttCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCGGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAgg                                                                                                                                   <  2:1044626/150‑2 (MQ=255)
                    tCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTAAACGGTGGTATACTGTGCAggt                                                                                                                                  <  2:751969/150‑3 (MQ=255)
                      aGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGaaaa                                                                                                                                >  2:346854/1‑150 (MQ=255)
                       gACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAAt                                                                                                                               <  2:1669606/150‑1 (MQ=255)
                        aCGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTCGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATg                                                                                                                              <  1:1472372/150‑1 (MQ=255)
                          gTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACAACCGTTTACTGTGCAGATAATGAc                                                                                                                            >  2:259886/1‑150 (MQ=255)
                           ttttAGATAATCGTCAGGTGAAAGATTGCTATGGTCACGGCCTAAAAATTCAATGGACAGCTTGGTGGGATCAAGTTGAGTAGTCATAGCATCCTCGCTTTTAGACAAGACCTGAACAGTATAGCACCGTATACTGTGCAGATAATGAcc                                                                                                                           <  1:98101/150‑1 (MQ=255)
                           tttGAGATATTGGGCAGGTGAAGGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGAcc                                                                                                                           <  1:498740/150‑1 (MQ=255)
                             tGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCaa                                                                                                                         <  2:1567780/150‑1 (MQ=255)
                                 atatTGCGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAAATCCATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACCTAGCATCCTCGCTTTTAGACAAGACCAGCACAGTATACCACCGTGTACTGTGCAGATAAAGACCaaaacc                                                                                                                     <  2:78325/150‑3 (MQ=255)
                                     tGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAata                                                                                                                 <  2:641083/150‑1 (MQ=255)
                                       ggCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTTCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATg                                                                                                               <  1:1506648/150‑1 (MQ=255)
                                       ggCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATg                                                                                                               >  1:536661/1‑150 (MQ=255)
                                        gCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATAAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGc                                                                                                              <  1:1520523/150‑1 (MQ=255)
                                           ggTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTc                                                                                                           >  2:680245/1‑150 (MQ=255)
                                           ggTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTc                                                                                                           <  1:471028/150‑1 (MQ=255)
                                            gTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTca                                                                                                          >  2:277441/1‑150 (MQ=255)
                                             tGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTcac                                                                                                         >  1:1005361/1‑150 (MQ=255)
                                              gAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTcaca                                                                                                        >  2:1381791/1‑150 (MQ=255)
                                                aaGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACt                                                                                                      >  1:993127/1‑150 (MQ=255)
                                                  gATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACttt                                                                                                    >  2:1002368/1‑150 (MQ=255)
                                                       ctctGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTgg                                                                                               <  1:1648966/150‑1 (MQ=255)
                                                            gTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACAACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACa                                                                                          >  1:1513624/1‑150 (MQ=255)
                                                             taaCGGAGTAAAAATAAAATCGCCAGCAGGGTCGGATCAAGTTCAGTATACATACCATCCTCCCTTTTAGACAAGACCTCCACATTATAACAACGATAACTCTGCAGATAATGACCAAAAGCAAAAACCTTCACACTATTCTGGtcaaaa                                                                                         >  2:726067/3‑145 (MQ=255)
                                                             tCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACaa                                                                                         >  1:1420060/1‑150 (MQ=255)
                                                              cACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAAc                                                                                        <  1:125748/150‑1 (MQ=255)
                                                               aCGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACg                                                                                       <  2:176058/150‑1 (MQ=255)
                                                                  gCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTca                                                                                    >  1:1447719/1‑150 (MQ=255)
                                                                  gCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTca                                                                                    >  1:1033173/1‑150 (MQ=255)
                                                                   cGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACACCGTcac                                                                                   >  2:985011/1‑150 (MQ=255)
                                                                      aaaaaTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACaaa                                                                                >  1:1232320/1‑150 (MQ=255)
                                                                       aaaaTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACaaaa                                                                               <  2:1703421/150‑1 (MQ=255)
                                                                        aaaTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAAt                                                                              >  2:38137/1‑150 (MQ=255)
                                                                        aaaTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAAt                                                                              <  1:346854/150‑1 (MQ=255)
                                                                           ttCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGc                                                                           <  1:1381791/150‑1 (MQ=255)
                                                                                       cTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCg                                                               <  2:1341163/150‑1 (MQ=255)
                                                                                       cTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCg                                                               <  2:554040/150‑1 (MQ=255)
                                                                                         gggtgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTg                                                             <  1:1345378/150‑1 (MQ=255)
                                                                                         gggtgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTg                                                             >  2:420587/1‑150 (MQ=255)
                                                                                          ggtgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGa                                                            <  2:385318/150‑1 (MQ=255)
                                                                                           gtgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGAc                                                           <  1:1662355/150‑1 (MQ=255)
                                                                                            tgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACg                                                          >  2:1474296/1‑150 (MQ=255)
                                                                                            tgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACg                                                          <  2:1270091/150‑1 (MQ=255)
                                                                                               gATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAAc                                                       >  2:678126/1‑150 (MQ=255)
                                                                                                aTCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACa                                                      <  1:277441/150‑1 (MQ=255)
                                                                                                    aGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTTACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGcac                                                  >  2:388355/1‑150 (MQ=255)
                                                                                                    aGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTATACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGcac                                                  <  2:1518596/150‑1 (MQ=255)
                                                                                                         aGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACg                                             <  1:680245/150‑1 (MQ=255)
                                                                                                           tAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGcc                                           >  1:1165325/1‑150 (MQ=255)
                                                                                                             gACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCt                                         <  2:1033173/150‑1 (MQ=255)
                                                                                                               cATAGCATCCTCGCTTTTAGCCAAGACCTGCACAGTATACCACCGTTAACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCttt                                       >  2:82181/1‑150 (MQ=255)
                                                                                                                 tAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTtct                                     <  2:1447719/150‑1 (MQ=255)
                                                                                                                 tAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTtct                                     >  2:440839/1‑150 (MQ=255)
                                                                                                                  aGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTtctc                                    <  1:420587/150‑1 (MQ=255)
                                                                                                                     aTCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATc                                 >  2:622793/1‑150 (MQ=255)
                                                                                                                       ccTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGa                               >  2:315964/1‑150 (MQ=255)
                                                                                                                       ccTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGa                               >  1:337620/1‑150 (MQ=255)
                                                                                                                          cGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGa                            <  1:1590654/150‑1 (MQ=255)
                                                                                                                            cTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGAtt                          >  1:288041/1‑150 (MQ=255)
                                                                                                                              tttAGACAAGACCTTCACAGTATACCACCGTTTACTCTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGAttta                        <  1:259886/150‑2 (MQ=255)
                                                                                                                                tAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCaa                      <  2:729268/150‑1 (MQ=255)
                                                                                                                                tAGACAAGACCTGCACAGGTTTCCACCGTTTACTGTGGAGATAATGACCTAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCACTCGTGGCGAGCAGGACAAGGGCGCTTTCTCAGCGGAGAtgtgat                      >  2:80111/1‑145 (MQ=255)
                                                                                                                                  gACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTTACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATc                    >  2:320535/1‑150 (MQ=255)
                                                                                                                                     aaGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGc                 <  2:1598189/150‑1 (MQ=255)
                                                                                                                                        atcTGCACAGTAAACGGTGGTATACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAg              >  2:1561693/3‑150 (MQ=255)
                                                                                                                                        aCCTGCACAGTATACCACCTTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAg              <  2:1263762/150‑1 (MQ=255)
                                                                                                                                          cTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGCTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAgcc            >  2:42415/1‑148 (MQ=255)
                                                                                                                                          cTGCACAGTAAACGGTGGTATACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAc            >  2:370940/1‑150 (MQ=255)
                                                                                                                                           tGCCCAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACCGCACCCACGCCCTTTCTCATCGAAGATTTCAAGCTGCCAGAcc           >  2:1491168/1‑150 (MQ=255)
                                                                                                                                           tGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAcc           >  1:26655/1‑150 (MQ=255)
                                                                                                                                           tGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAcc           <  1:320535/150‑1 (MQ=255)
                                                                                                                                                 gTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGa     <  2:406040/150‑1 (MQ=255)
                                                                                                                                                  tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAc    <  1:315964/150‑1 (MQ=255)
                                                                                                                                                   ataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcg   >  2:903955/1‑150 (MQ=255)
                                                                                                                                                   ataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcg   <  1:1085155/150‑1 (MQ=255)
                                                                                                                                                    taCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcgc  >  1:1328666/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAATTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGC  >  NC_000913/1765009‑1765306

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