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

AGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCG  >  NC_000913/1765016‑1765312
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                                 <  1:471028/150‑1 (MQ=255)
                                    ggTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTc                                                                                                                 >  2:680245/1‑150 (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                                                                                    <  1:346854/150‑1 (MQ=255)
                                                                 aaaTTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAAt                                                                                    >  2:38137/1‑150 (MQ=255)
                                                                    ttCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGc                                                                                 <  1:1381791/150‑1 (MQ=255)
                                                                                cTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCg                                                                     <  2:554040/150‑1 (MQ=255)
                                                                                cTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCg                                                                     <  2:1341163/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:1270091/150‑1 (MQ=255)
                                                                                     tgggATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACg                                                                >  2:1474296/1‑150 (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                                     >  1:337620/1‑150 (MQ=255)
                                                                                                                ccTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGa                                     >  2:315964/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)
                                                                                                                                                  ccGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGGCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAAcc   <  1:618686/150‑1 (MQ=255)
                                                                                                                                                   cGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGATGCGAACCg  >  2:1151999/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTTCTAGTTGTTTCAGACGTTTGAGATATTGGGCAGGTGAAAGATTGCTCTGGTCACGGCGTAAAAATTCAATGGCCAGCTGGGTGGGATCAAGTTGAGTAGACATAGCATCCTCGCTTTTAGACAAGACCTGCACAGTATACCACCGTTTACTGTGCAGATAATGACCAAAAGCAATATGCGTCACACTTTTCTGGTGACAACGTCACAAAATGGCGGTCGTCAATCGTGACGAACAGCACAAACGCCCTTTCTCATCGAAGATTTCAATCTGCCAGACCTGGTGACGCGAACCG  >  NC_000913/1765016‑1765312

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: