Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,830 A→G 6.9% intergenic (+22/+18) 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,8300AG 6.9% 512.6 / 18.3 159intergenic (+22/+18)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (76/72);  new base G (6/5);  total (82/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 160 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                                                                                                                                                                                                                                                                   >  5:711312/1‑149 (MQ=255)
                                                                                                                            gCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                                                                                                                                                                                                                                   <  6:282898/149‑1 (MQ=255)
                                                                                                                            gCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCATGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                                                                                                                                                                                                                                   >  6:849624/1‑149 (MQ=255)
                                                                                                                            gCGAACATGTGTGAGCGTGTGGTGAAAAATGCGATCCGCAGATGGCGTTATGGGCCGGGTAAGCCAGGCATTGGGATTGTGGTGAATAACCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtt                                                                                                                                                                                                                                                          <  4:17480/158‑1 (MQ=255)
                                                                                                                             cGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                                                                                                                                                                                                                                  <  8:678049/149‑1 (MQ=255)
                                                                                                                              gAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                                                                                                                                                                                 >  8:1114960/1‑149 (MQ=255)
                                                                                                                                aCATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                                                                                                                                                                               <  8:970258/149‑1 (MQ=255)
                                                                                                                                     tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtt                                                                                                                                                                                                                                                          >  7:991971/1‑149 (MQ=255)
                                                                                                                                      ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCttt                                                                                                                                                                                                                                                         <  5:584797/149‑1 (MQ=255)
                                                                                                                                      ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                                                                                                                                                                                         <  6:1379450/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  7:998931/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  5:1261893/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  5:1271900/149‑1 (MQ=255)
                                                                                                                                       tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  7:287557/149‑1 (MQ=255)
                                                                                                                                            gTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTa                                                                                                                                                                                                                                                   >  7:62045/1‑149 (MQ=255)
                                                                                                                                             gggggggAAAAAGGCGGTGGGAAGATGGNTNTGTGGGCCGGGTAAGCCAGGCAGTGGGATTGTGGTAAATATCCTTNTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACNATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGTGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAata                                                                                                                    <  3:30777/269‑1 (MQ=255)
                                                                                                                                                 gTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACt                                                                                                                                                                                                                                              >  8:477025/1‑149 (MQ=255)
                                                                                                                                                  tGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                                                                                                                                      <  8:28344/140‑1 (MQ=255)
                                                                                                                                                  tGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                                                                                                                                      >  7:28344/1‑140 (MQ=255)
                                                                                                                                                         tgcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGAc                                                                                                         <  2:62954/274‑1 (MQ=255)
                                                                                                                                                           cgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAGCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTTCCCAATTGGTTCAGAAGCTACTTTTTTAACCCAAACTTCCATTTTATTGCTGACGGa                                                                                                      >  2:305785/1‑275 (MQ=255)
                                                                                                                                                            gatgcgCAGATGGCTTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTAc                                                                                                                                                                                                                                               <  5:478483/137‑1 (MQ=255)
                                                                                                                                                            gatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAg                                                                                                                                                                                                                                   >  6:519863/1‑149 (MQ=255)
                                                                                                                                                            gatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTAAAAAGCCTCCGACCGGAGGCTTTTGACTATTAc                                                                                                                                                                                                                                               >  6:478483/1‑137 (MQ=255)
                                                                                                                                                                   aGATGGGGTTATGGGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTcc                                                                                              <  2:212855/275‑1 (MQ=255)
                                                                                                                                                                   aGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            <  6:809919/149‑1 (MQ=255)
                                                                                                                                                                      tGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTCAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            <  8:997751/146‑1 (MQ=255)
                                                                                                                                                                        gCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                                                                                                                                                                                       >  5:1477689/1‑149 (MQ=255)
                                                                                                                                                                         cGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCATGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                                                                                                                                                                                      <  5:849624/149‑1 (MQ=255)
                                                                                                                                                                          gTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                                                                                                                                         <  5:521819/129‑1 (MQ=255)
                                                                                                                                                                           ttATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                                                                                                                                                                                                    >  8:1354216/1‑149 (MQ=255)
                                                                                                                                                                             aTGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTc                                                                                                                                                                                                                  <  8:798037/149‑1 (MQ=255)
                                                                                                                                                                               gAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                                                                                                                                                                                <  5:276551/149‑1 (MQ=255)
                                                                                                                                                                                 gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              >  5:77887/1‑149 (MQ=255)
                                                                                                                                                                                 gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              <  8:661532/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)
                                                                                                                                                                                         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                                                                                                                                                                                      >  7:300392/1‑149 (MQ=255)
                                                                                                                                                                                            cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAgg                                                                                                                                                                                                                                                             >  5:352358/1‑91 (MQ=255)
                                                                                                                                                                                            cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                                                                                                                                        >  8:1163565/1‑128 (MQ=255)
                                                                                                                                                                                             aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                                                                                                                                                                                  >  7:1396138/1‑149 (MQ=255)
                                                                                                                                                                                                     ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          >  7:873194/1‑149 (MQ=255)
                                                                                                                                                                                                       aTTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaata                                                                                                                                                                                        >  8:750159/1‑149 (MQ=255)
                                                                                                                                                                                                         tgtgGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                                                                                                                                                                                      <  7:794963/149‑1 (MQ=255)
                                                                                                                                                                                                           tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              <  8:1439828/123‑1 (MQ=255)
                                                                                                                                                                                                           tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              >  7:1439828/1‑123 (MQ=255)
                                                                                                                                                                                                           tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGc                                                                                                                                                                                    >  7:1108699/1‑149 (MQ=255)
                                                                                                                                                                                                              tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                                                                                                                                      >  7:1053140/1‑80 (MQ=255)
                                                                                                                                                                                                              tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                                                                                                                                      <  8:1053140/80‑1 (MQ=255)
                                                                                                                                                                                                                aaTATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAAAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGCCATACTTAAATAATGCTTCTg                                                                                                                                                                               >  7:507561/1‑149 (MQ=255)
                                                                                                                                                                                                                    tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                                                                                                                                                                                >  6:1458482/1‑144 (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)
                                                                                                                                                                                                                        gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATAc                                                                                                                                                                                               >  7:7376/1‑125 (MQ=255)
                                                                                                                                                                                                                        gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCCCACCCCATCGCCGACCGCAACCCGc                                         >  4:260780/1‑275 (MQ=255)
                                                                                                                                                                                                                           tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgc                                                                                                                                                                                                                                 <  8:365759/88‑1 (MQ=255)
                                                                                                                                                                                                                              aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGt                                                                                                                                                                 <  8:189962/149‑1 (MQ=255)
                                                                                                                                                                                                                              aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGt                                                                                                                                                                 >  7:685448/1‑149 (MQ=255)
                                                                                                                                                                                                                              aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTAGTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAg                                   <  4:112248/275‑1 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  7:41813/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  8:1417772/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  5:303655/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  6:772886/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCg                                                                              >  2:33233/1‑227 (MQ=255)
                                                                                                                                                                                                                                     ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCaa                                                                                                                                                          <  6:769603/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAAt                                                                                                                                                         <  8:1060395/149‑1 (MQ=255)
                                                                                                                                                                                                                                       caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCa                                                                          <  4:33233/227‑1 (MQ=255)
                                                                                                                                                                                                                                       caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGGCGCTTTATAGCGTTGCCCAAtt                                                                                                                                                        <  5:900032/149‑1 (MQ=255)
                                                                                                                                                                                                                                           accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                                                                                                                                        >  7:486903/1‑81 (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)
                                                                                                                                                                                                                                                  ttCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGc                                                                                                                                             >  7:1305449/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTAc                                                                                                                                          <  7:796561/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACt                                                                                                                                         >  6:564332/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTc                                                                                                                                   <  6:1477689/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGc                                                                                                                                                              >  6:1407255/1‑121 (MQ=255)
                                                                                                                                                                                                                                                             gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGc                                                                                                                                                              <  5:1407255/121‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaGTCAATGCCCTCCCACCTGGGGCTTTTTGCTATTACTCAACAGGTAAGGCGGGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcg  <  4:21150/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          <  8:1193448/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                                                                                                                                                                                          >  7:1193448/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCTCCCacac                                                                                                                           >  7:372714/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACt                                                                                                                          >  5:1002154/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt                                                                                                                         >  5:1424228/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)
                                                                                                                                                                                                                                                                             cGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATAtt                                                                                                                  >  8:828967/1‑149 (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: ‑