Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,765,160 C→G 7.9% intergenic (‑38/+62) 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,1600CG 7.9% 251.7 / 11.4 88intergenic (‑38/+62)ydiH/menIuncharacterized protein/1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base C (37/42);  new base G (2/5);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

ATTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAA  >  NC_000913/1765011‑1765309
                                                                                                                                                     |                                                                                                                                                     
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:320535/150‑1 (MQ=255)
                                                                                                                                         tGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGAcc              >  1:26655/1‑150 (MQ=255)
                                                                                                                                               gTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGa        <  2:406040/150‑1 (MQ=255)
                                                                                                                                                tataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAc       <  1:315964/150‑1 (MQ=255)
                                                                                                                                                 ataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcg      <  1:1085155/150‑1 (MQ=255)
                                                                                                                                                 ataCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcg      >  2:903955/1‑150 (MQ=255)
                                                                                                                                                  taCCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcgc     >  1:1328666/1‑150 (MQ=255)
                                                                                                                                                   accaccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTTACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGAcgcg    <  1:388355/150‑1 (MQ=255)
                                                                                                                                                    ccaccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGa   >  1:1551541/1‑150 (MQ=255)
                                                                                                                                                     caccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGaa  >  2:1088191/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCAAGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAA  >  NC_000913/1765011‑1765309

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: