Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,831 C→T 7.0% intergenic (+23/+17) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8310CT 7.0% 506.8 / 15.4 156intergenic (+23/+17)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (74/71);  new base T (6/5);  total (80/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 161 total aligned reads displayed.

CGTAATCAGCCGCAGTATCCGGCACGAGCACA‑GGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACG  >  NC_000913/1311560‑1312088
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                 
cGTAATCAGCCGCAGTATCCGGCACGGGCACA‑GGCATTGCGCATTGAAGGGCAGGGTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGCGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                                                                                                                                                                                <  2:103283/275‑1 (MQ=255)
cGTAATCAGCCGAAGGTGCCGGCAAG‑GGACATGGCATTGGGCAGTGAAGGGAAGTGTAAAGTTAAATTTGACGGTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCGCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                                                                                                                                                                                <  4:335766/275‑1 (MQ=255)
   aaTCAGCCGCCGTATCCGGCACGAGCACA‑GGCATTGCGCATTGAATGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGCTAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGCGAGGTGAAAATTGCGATGCGCAGATGGCGTTATGAGACGCGTAAGCCACGCAGTGGGATTGTGGTGACTATCCTGTTTAAAATTAACGCCACCACCGAACTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACAGGa                                                                                                                                                                                                                                                               >  1:67796/1‑273 (MQ=255)
                                                                                  taTGGTGGGGTGGATAACGTACAAATCCTCTCAGCCAGGCCTGCGGACAGGTTTGAGCGTGAGGTGAAAATTGCGATNCGGAGGTGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGAGTTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTg                                                                                                                                                                               <  3:67796/274‑1 (MQ=255)
                                                                                        cgcgTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTCTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCtt                                                                                                                                                                         >  2:112248/1‑275 (MQ=255)
                                                                                                                        gCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTCTGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCACAAGCCTCCGACCggcg                                                                                                                                                                                                                                                              >  2:17480/1‑156 (MQ=255)
                                                                                                                            gCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                                                                                                                                                                                                                                   <  6:282898/149‑1 (MQ=255)
                                                                                                                            gCGAACATGTGTGAGCGTGTGGTGAAAAATGCGATCCGCAGATGGCGTTATGGGCCGGGTAAGCCAGGCATTGGGATTGTGGTGAATAACCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtt                                                                                                                                                                                                                                                          <  4:17480/158‑1 (MQ=255)
                                                                                                                             cGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                                                                                                                                                                                                                                  <  8:678049/149‑1 (MQ=255)
                                                                                                                                aCATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                                                                                                                                                                               <  5:728111/149‑1 (MQ=255)
                                                                                                                                aCATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                                                                                                                                                                               <  8:970258/149‑1 (MQ=255)
                                                                                                                                     tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAAATTACGGCACCACCGAAATTCCGTTAGCCAAAAATTAAAAGCCTTCCACCCGAAGCtt                                                                                                                                                                                                                                                          >  8:796561/1‑149 (MQ=255)
                                                                                                                                      ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCttt                                                                                                                                                                                                                                                         <  5:584797/149‑1 (MQ=255)
                                                                                                                                      ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                                                                                                                                                                                         <  6:1379450/149‑1 (MQ=255)
                                                                                                                                      ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGCAGGCttt                                                                                                                                                                                                                                                         <  6:809332/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                                                                                                                                        <  8:1099485/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  7:998931/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  6:569560/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  5:1261893/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  5:1271900/149‑1 (MQ=255)
                                                                                                                                          gCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                                                                                                                                                                                     >  8:1327031/1‑149 (MQ=255)
                                                                                                                                             gggggggAAAAAGGCGGTGGGAAGATGGNTNTGTGGGCCGGGTAAGCCAGGCAGTGGGATTGTGGTAAATATCCTTNTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACNATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGTGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAata                                                                                                                    <  3:30777/269‑1 (MQ=255)
                                                                                                                                                  tGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                                                                                                                                      >  7:28344/1‑140 (MQ=255)
                                                                                                                                                         tgcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGAc                                                                                                         <  2:62954/274‑1 (MQ=255)
                                                                                                                                                           cgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAGCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTTCCCAATTGGTTCAGAAGCTACTTTTTTAACCCAAACTTCCATTTTATTGCTGACGGa                                                                                                      >  2:305785/1‑275 (MQ=255)
                                                                                                                                                            gatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTAAAAAGCCTCCGACCGGAGGCTTTTGACTATTAc                                                                                                                                                                                                                                               >  6:478483/1‑137 (MQ=255)
                                                                                                                                                                   aGATGGGGTTATGGGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTcc                                                                                              <  2:212855/275‑1 (MQ=255)
                                                                                                                                                                   aGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            <  6:809919/149‑1 (MQ=255)
                                                                                                                                                                   aGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            >  7:661532/1‑149 (MQ=255)
                                                                                                                                                                     aTGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGt                                                                                                                                                                                                                          >  6:994971/1‑149 (MQ=255)
                                                                                                                                                                      tGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTCAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            <  8:997751/146‑1 (MQ=255)
                                                                                                                                                                      tGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGATTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            >  7:997751/1‑146 (MQ=255)
                                                                                                                                                                        gCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                                                                                                                                                                                       >  5:1477689/1‑149 (MQ=255)
                                                                                                                                                                         cGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                                                                                                                                                                                      <  8:1150082/149‑1 (MQ=255)
                                                                                                                                                                         cGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCATGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                                                                                                                                                                                      <  5:849624/149‑1 (MQ=255)
                                                                                                                                                                          gTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                                                                                                                                         <  5:521819/129‑1 (MQ=255)
                                                                                                                                                                          gTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGAGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                                                                                                                                         >  6:521819/1‑129 (MQ=255)
                                                                                                                                                                               gAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                                                                                                                                                                                <  5:276551/149‑1 (MQ=255)
                                                                                                                                                                                 gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              >  5:77887/1‑149 (MQ=255)
                                                                                                                                                                                 gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              <  8:661532/149‑1 (MQ=255)
                                                                                                                                                                                  ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAt                                                                                                                                                                                                             <  5:519863/149‑1 (MQ=255)
                                                                                                                                                                                   cGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATc                                                                                                                                                                                                            >  7:981076/1‑149 (MQ=255)
                                                                                                                                                                                     ggTAAGCCACGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACAGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTg                                                                            <  4:77629/275‑1 (MQ=255)
                                                                                                                                                                                       tAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAcc                                                                                                                                                                                                        <  8:1101518/149‑1 (MQ=255)
                                                                                                                                                                                         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAt                                                                                                                                                                                                             <  2:292160/142‑1 (MQ=255)
                                                                                                                                                                                            cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAgg                                                                                                                                                                                                                                                             >  5:352358/1‑91 (MQ=255)
                                                                                                                                                                                            cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                                                                                                                                        <  7:1163565/128‑1 (MQ=255)
                                                                                                                                                                                             aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                                                                                                                                                                                  >  7:1396138/1‑149 (MQ=255)
                                                                                                                                                                                                     ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          <  6:1408774/149‑1 (MQ=255)
                                                                                                                                                                                                       aTTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaata                                                                                                                                                                                        >  8:750159/1‑149 (MQ=255)
                                                                                                                                                                                                         tgtgGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                                                                                                                                                                                      <  7:794963/149‑1 (MQ=255)
                                                                                                                                                                                                             gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt                                                                                                                                                                                  >  8:28084/1‑149 (MQ=255)
                                                                                                                                                                                                              tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                                                                                                                                      >  7:1053140/1‑80 (MQ=255)
                                                                                                                                                                                                                    tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                                                                                                                                                                                >  6:1458482/1‑144 (MQ=255)
                                                                                                                                                                                                                    tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                                                                                                                                                                                <  5:1458482/144‑1 (MQ=255)
                                                                                                                                                                                                                     ccTGTTTAGAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCt                                                                                                                                                                          <  5:576912/149‑1 (MQ=255)
                                                                                                                                                                                                                     ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCt                                                                                                                                                                          >  7:1349869/1‑149 (MQ=255)
                                                                                                                                                                                                                      cTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCtt                                                                                                                                                                         <  8:271000/149‑1 (MQ=255)
                                                                                                                                                                                                                      cTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCtt                                                                                                                                                                         <  8:981076/149‑1 (MQ=255)
                                                                                                                                                                                                                       tGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt                                                                                                                                                                        <  7:28084/149‑1 (MQ=255)
                                                                                                                                                                                                                        gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCCCACCCCATCGCCGACCGCAACCCGc                                         >  4:260780/1‑275 (MQ=255)
                                                                                                                                                                                                                         tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtat                                                                                                                                                                      >  7:1321599/1‑149 (MQ=255)
                                                                                                                                                                                                                         tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtat                                                                                                                                                                      >  6:1097714/1‑149 (MQ=255)
                                                                                                                                                                                                                         tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtat                                                                                                                                                                      >  7:292050/1‑149 (MQ=255)
                                                                                                                                                                                                                           tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgc                                                                                                                                                                                                                                 <  8:365759/88‑1 (MQ=255)
                                                                                                                                                                                                                           tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgc                                                                                                                                                                                                                                 >  7:365759/1‑88 (MQ=255)
                                                                                                                                                                                                                              aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGt                                                                                                                                                                 <  7:859767/149‑1 (MQ=255)
                                                                                                                                                                                                                              aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGt                                                                                                                                                                 >  7:685448/1‑149 (MQ=255)
                                                                                                                                                                                                                              aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTAGTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAg                                   <  4:112248/275‑1 (MQ=255)
                                                                                                                                                                                                                               aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt                                                                                                                                                                <  7:750159/149‑1 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  8:1417772/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  7:443645/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  5:303655/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCg                                                                              >  2:33233/1‑227 (MQ=255)
                                                                                                                                                                                                                                     ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCaa                                                                                                                                                          <  6:769603/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAAt                                                                                                                                                         <  8:300392/149‑1 (MQ=255)
                                                                                                                                                                                                                                       caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCa                                                                          <  4:33233/227‑1 (MQ=255)
                                                                                                                                                                                                                                           accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                                                                                                                                        <  8:486903/81‑1 (MQ=255)
                                                                                                                                                                                                                                           accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                                                                                                                                        >  7:486903/1‑81 (MQ=255)
                                                                                                                                                                                                                                           accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                                                                                                                                                    >  7:591728/1‑149 (MQ=255)
                                                                                                                                                                                                                                           accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                                                                                                                                                    >  5:1293499/1‑149 (MQ=255)
                                                                                                                                                                                                                                           accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                                                                                                                                                    <  7:1077103/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTc                                                                                                                                                   >  6:1158762/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 aTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACg                                                                                                                                              >  6:530003/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  ttCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGc                                                                                                                                             <  5:747254/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCTCAATTGGTTCAGACGCTACt                                                                                                                                         >  8:1447480/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGc                                                                                                                                                              <  5:1407255/121‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAc                                                                                                                                 <  8:1059976/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaGTCAATGCCCTCCCACCTGGGGCTTTTTGCTATTACTCAACAGGTAAGGCGGGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcg  <  4:21150/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAccc                                                                                                                               >  7:973481/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          >  7:1193448/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          <  8:1193448/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCTCCCacac                                                                                                                           >  7:372714/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACt                                                                                                                          >  5:1002154/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          >  6:578122/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                      aaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          <  5:578122/84‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt                                                                                                                         >  5:807466/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAat                                                                                                                     <  6:238146/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATAtt                                                                                                                  >  8:828967/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               aCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCGACACTTCAATATTaa                                                                                                                <  8:1396138/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                 
CGTAATCAGCCGCAGTATCCGGCACGAGCACA‑GGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACG  >  NC_000913/1311560‑1312088

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