Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG46.4% 184.4 / 30.0 110D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (8/51);  new base G (51/0);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 3.92e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGG  >  NC_000913/827357‑827656
                                                                                                                                                      |                                                                                                                                                      
gTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGt                                                                                                                                                        <  2:14949/151‑1 (MQ=255)
  gCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTc                                                                                                                                                       <  2:283883/150‑1 (MQ=255)
   ccGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCaa                                                                                                                                                     <  2:1159540/151‑1 (MQ=255)
    cGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAAc                                                                                                                                                    <  1:1445999/151‑1 (MQ=255)
    cGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAAc                                                                                                                                                    <  1:858452/151‑1 (MQ=255)
    cGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCaa                                                                                                                                                     >  1:1157927/1‑150 (MQ=255)
        ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGc                                                                                                                                                  <  2:874356/149‑1 (MQ=255)
        ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCg                                                                                                                                                <  2:1357284/151‑1 (MQ=255)
            aTTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAg                                                                                                                                            <  1:1596812/151‑1 (MQ=255)
            aTTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAg                                                                                                                                            <  1:922395/151‑1 (MQ=255)
              tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGt                                                                                                                                           <  1:1421931/150‑1 (MQ=255)
              tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGt                                                                                                                                           <  1:962604/150‑1 (MQ=255)
               aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTc                                                                                                                                          <  2:723183/150‑1 (MQ=255)
               aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCg                                                                                                                                         <  1:1363552/151‑1 (MQ=255)
               aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCg                                                                                                                                         >  1:1171300/1‑151 (MQ=255)
                aTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCgg                                                                                                                                        <  2:1022890/151‑1 (MQ=255)
                    cccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcg                                                                                                                                    <  1:13426/151‑1 (MQ=255)
                    cccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcg                                                                                                                                    <  1:1095354/151‑1 (MQ=255)
                     ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGtcg                                                                                                                                    >  2:939489/1‑150 (MQ=255)
                      cgcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtc                                                                                                                                  <  1:1570733/151‑1 (MQ=255)
                      cgcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtc                                                                                                                                  <  1:1339535/151‑1 (MQ=255)
                        cgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCga                                                                                                                                <  1:1499235/151‑1 (MQ=255)
                           tACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTTGTCgag                                                                                                                               >  1:1352482/1‑149 (MQ=255)
                            acacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaaa                                                                                                                            <  1:1735140/151‑1 (MQ=255)
                              acCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGaaa                                                                                                                            >  2:1638706/1‑149 (MQ=255)
                               ccAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAg                                                                                                                         <  1:863401/151‑1 (MQ=255)
                                 aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGTCGGGTCGTCGAGAAACa                                                                                                                          >  2:808677/1‑148 (MQ=255)
                                  ggCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAAAAGaa                                                                                                                       >  1:1454246/1‑150 (MQ=255)
                                  ggCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTTGTCGAGAAACAGaa                                                                                                                       >  1:661076/1‑150 (MQ=255)
                                    ccGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGTTTGTCGAGAAAAAAAATGt                                                                                                                    >  1:912318/1‑151 (MQ=255)
                                     cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTc                                                                                                                   >  1:311021/1‑151 (MQ=255)
                                     cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAAAAGAATg                                                                                                                     >  1:1163888/1‑149 (MQ=255)
                                      gATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGGAGTCGGGTTGTCGAGAAACAGAATGTcc                                                                                                                  >  1:816955/1‑151 (MQ=255)
                                               gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGc                                                                                                         >  1:1180836/1‑151 (MQ=255)
                                                    atTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGGAAAAGAATGGCCGGGTcct                                                                                                           >  2:600645/1‑142 (MQ=255)
                                                    atTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGGCGGGTTGTCGAGAAACAGAAAGTCCGGGTTATg                                                                                                          >  1:559865/1‑145 (MQ=255)
                                                         gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                                                                               <  2:347766/151‑1 (MQ=255)
                                                         gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTTGAGAAACAGAAAGTCCGGGTTATGCATCCGCGaa                                                                                                >  2:1746769/1‑150 (MQ=255)
                                                         gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGGAAAAGAAAGTCCGGGTTATGCATCAGCGaa                                                                                                >  1:257409/1‑150 (MQ=255)
                                                           ttCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAACg                                                                                                  >  2:326955/1‑146 (MQ=255)
                                                            tCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAgg                                                                                            >  1:1006092/1‑151 (MQ=255)
                                                                 cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGGGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGgcc                                                                                          >  1:388180/1‑148 (MQ=255)
                                                                  aTAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCTTGCATCAGCgcaaaagcccg                                                                                        >  1:950562/1‑139 (MQ=255)
                                                                    aaaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcc                                                                                    <  2:1697862/151‑1 (MQ=255)
                                                                     aaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGTGGGGGGGCACACCGGAAATAGGTTCCTCGAAAAACAAAAATTTCGGGTCCTGCCTCAGGCTAAAGgcccgcgc                                                                                     >  2:1019175/1‑149 (MQ=255)
                                                                       gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgccagc                                                                                 <  2:1181759/151‑1 (MQ=255)
                                                                         gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCg                                                                               <  1:622072/151‑1 (MQ=255)
                                                                         gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGGTCGTTGAGAAAAAGAATGGCCGGGTCATGCATCAACGAACAGGCCCGCGCCAGcc                                                                                >  2:127825/1‑150 (MQ=255)
                                                                             ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCCACCCCCGGAGGCGGTTCCTCGAGAAACAGAATGTCCGGGTCATGCCTCAGCGGACAGGCCAGCCCCcgccgc                                                                             >  2:1775706/1‑148 (MQ=255)
                                                                             ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGt                                                                          <  2:1157855/151‑1 (MQ=255)
                                                                              gtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTTGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGGACAGGCCCGcgcccgcc                                                                                >  1:1072171/1‑145 (MQ=255)
                                                                                gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTa                                                                       <  2:846184/151‑1 (MQ=255)
                                                                                gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCCTCAGCGAACAGGCCAGCGCCAGCCGCTGttt                                                                        >  1:1161085/1‑150 (MQ=255)
                                                                                 aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGttt                                                                        >  2:620922/1‑149 (MQ=255)
                                                                                 aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa                                                                      <  2:277593/151‑1 (MQ=255)
                                                                                 aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGCGAAACAGAATGTCCCGGTCATGCATCAGCGAACAGGCCCGcgcccgcccccggtt                                                                         >  2:859866/1‑142 (MQ=255)
                                                                                   cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAAGCCAGCGCCAGCCGCTGTTTtaa                                                                     >  1:1426352/1‑150 (MQ=255)
                                                                                   cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGGCCCGTTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCCGTTTaaa                                                                     >  1:135139/1‑150 (MQ=255)
                                                                                   cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTTGGTTCGTCGAGAAACAAAATGGCCGGTTCCCGCATCTGCGAACACGCCAACGCCAAGCGCTGTTTaaa                                                                     >  1:272979/1‑150 (MQ=255)
                                                                                    aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTTATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTaaaa                                                                    >  1:710888/1‑150 (MQ=255)
                                                                                     tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTTGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGGACAGGGCAGCGCCCGCCGCTGTTTaaaa                                                                    >  1:1228508/1‑149 (MQ=255)
                                                                                        ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCCGAAGTCGGTTCGTCGAGAAACAAAATGTCCGGTTCATGCATCAGCGAACAGGCCCGCGCCCGCCGCTGTTTaaaaac                                                                  >  1:382754/1‑146 (MQ=255)
                                                                                         cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCCGCGCCCGCCCCCGTTTTAAAACtca                                                               >  2:216575/1‑148 (MQ=255)
                                                                                          gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAAAGCCGGGAGTGGGGTCGTCGAGAAAAAGAATGTTCGGGTTATGGAACAACGAAAAAGCCCGCGCCAGCCGCCGGTTTAAAAccaa                                                                >  2:866857/1‑147 (MQ=255)
                                                                                            gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                           <  2:381078/151‑1 (MQ=255)
                                                                                             aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAATGGGCGGGTCATGCATCAGCGGAAAGGGCCGCGCCCGCCGCTGTTTTAAAACTAATg                                                             >  2:50222/1‑148 (MQ=255)
                                                                                             aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCa                                                          <  1:298394/151‑1 (MQ=255)
                                                                                             aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAAGGGGGGCACGCCGGGAGTCGGGTCGGCCAGAAACAGAGTGTGCCGGTCGTGTATCAGCGGAAAGGCCAGGGCCAGCCGCTCGTTGAAACCTaa                                                               >  2:807392/1‑146 (MQ=255)
                                                                                              cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGCGGGGGGCCACGCCGGAAGTCCGGTTGTCGCGAAAAAGAATGTCCGGGTCATGGCTCCGCGAAAAAGCCCGCGCCCACCGCCGTTTTAAAACTAAttg                                                            >  2:750800/1‑146 (MQ=255)
                                                                                               gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGAAGGCGGTTCGTCGAGAAACAAAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCCGTTTAAAACCTAATgg                                                            >  1:198728/1‑147 (MQ=255)
                                                                                                ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGtt                                                       <  1:181066/151‑1 (MQ=255)
                                                                                                  tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGGCCGCGCCAGCCGCTGTTTTAAACCTAATGGccggt                                                       >  1:217168/1‑145 (MQ=255)
                                                                                                       ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTTGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCg                                                  <  2:469177/149‑1 (MQ=255)
                                                                                                       ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCg                                                  <  2:57048/149‑1 (MQ=255)
                                                                                                         aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTgg                                              <  1:513794/151‑1 (MQ=255)
                                                                                                          ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAACGAACAAGCCAGCGCCCGCCGCTGTTTAAAACCTAATGGCAGTTCATCGgggg                                              >  1:1230649/1‑150 (MQ=255)
                                                                                                           cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCg                                            <  1:1587902/151‑1 (MQ=255)
                                                                                                           cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAACCTAATGGCAGTTCATCGGTgg                                              >  2:620978/1‑149 (MQ=255)
                                                                                                            aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAACTAATGGCAGTTCATCGGGGGCg                                            >  1:1066952/1‑150 (MQ=255)
                                                                                                              gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAAGGCCAGCGCCAGCCGCCGTTTAAAAACTAATGGGAGTTCATCGGGGGCgggg                                         >  1:41073/1‑151 (MQ=255)
                                                                                                               cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGGCGGGTCGTCGCGAAAAAGAAAGTCCCGGTCATGCCTCAGCGAACAGGCCAGCGCCCGCCGCTGTTTAAAACCTAATGGCAGCTCCTCGgggggggggg                                        >  1:1212966/1‑151 (MQ=255)
                                                                                                                   tGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCCCTGTTTTAAAACTAATGGCCGGTCATCGGTGGCgggggggg                                     >  1:844939/1‑150 (MQ=255)
                                                                                                                   tGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTCCGGTTCATGCCTCAGCGAACAGGCCCGCGCCAGCCGCTGTTTAAAAACCAATGGCAGTTCATCGGGGGCgggggggg                                     >  2:944399/1‑150 (MQ=255)
                                                                                                                       gtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAt                                 <  2:1610161/150‑1 (MQ=255)
                                                                                                                                  aTTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCAACAGCGGACAGGCCAGCGCCAGCCGCTGTTTTAAAACTAATGGCAGGTCATCGGGGGCGGGGGGGGCGGTAATTTTCCGGcc                     >  2:1036409/1‑151 (MQ=255)
                                                                                                                                     cACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGaa                  <  2:112814/151‑1 (MQ=255)
                                                                                                                                     cACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGaa                  <  1:199913/151‑1 (MQ=255)
                                                                                                                                      aCGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCCTGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAAACTAATGGCAGGTCATCCGTGGCGGGGGGGGGGATACTTTTCCGGcccga                  >  2:173927/1‑147 (MQ=255)
                                                                                                                                        ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGc               >  1:751552/1‑151 (MQ=255)
                                                                                                                                          cggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCt             >  2:1765008/1‑151 (MQ=255)
                                                                                                                                           gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTc            <  2:231344/151‑1 (MQ=255)
                                                                                                                                           gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTc            >  1:875617/1‑151 (MQ=255)
                                                                                                                                           gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTc            <  1:1313570/151‑1 (MQ=255)
                                                                                                                                            ggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCg           <  2:627637/151‑1 (MQ=255)
                                                                                                                                               gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCt         <  1:1514499/150‑1 (MQ=255)
                                                                                                                                               gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTc        <  1:310994/151‑1 (MQ=255)
                                                                                                                                                aGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCa       <  2:1738457/151‑1 (MQ=255)
                                                                                                                                                  ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg     <  1:1246841/151‑1 (MQ=255)
                                                                                                                                                   gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGc    >  1:272549/1‑151 (MQ=255)
                                                                                                                                                     gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCgg  <  2:822029/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGG  >  NC_000913/827357‑827656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: