Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,236,9010TC39.2% 100.3 / 16.6 74T285A (ACG→GCG) mglCmethyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (30/15);  new base C (0/29);  total (30/44)
Fisher's exact test for biased strand distribution p-value = 7.98e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAA  >  NC_000913/2236753‑2237045
                                                                                                                                                     |                                                                                                                                                
cgcgcgTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTcc                                                                                                                                                <  1:287704/151‑1 (MQ=255)
           cAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACAccgccg                                                                                                                                     >  2:845008/1‑151 (MQ=255)
           cAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACAccgcc                                                                                                                                      <  2:148056/150‑1 (MQ=255)
              tGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTg                                                                                                                                  >  2:698217/1‑151 (MQ=255)
                 aTCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAAc                                                                                                                               >  1:783443/1‑151 (MQ=255)
                   ccAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACg                                                                                                                              >  1:860666/1‑150 (MQ=255)
                     aGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGAt                                                                                                                            >  2:38316/1‑150 (MQ=255)
                      gcgtTACGGCGAAGATAATAATCGCCCCTTTGATGATTTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGAt                                                                                                                            <  1:435165/149‑1 (MQ=255)
                               cGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTAc                                                                                                                 >  1:328545/1‑151 (MQ=255)
                                 aaGATAATACTCGCCCCTTTGATGTGGTACTGCCAGTATGGTTTTACGCCAATATACGTCAGCCCATAGTTGATGAGGGTAAAAATAATTCCCCCGTCCCCCCCCCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTAc                                                                                                                 <  2:624493/149‑1 (MQ=255)
                                 aaGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTAcc                                                                                                                >  1:915999/1‑150 (MQ=255)
                                   gATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACcacg                                                                                                             >  2:43299/1‑151 (MQ=255)
                                    ataataATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACcacgc                                                                                                            >  2:1091609/1‑151 (MQ=255)
                                     taataaTCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACcacgc                                                                                                            >  2:81139/1‑150 (MQ=255)
                                         aatCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTCCCCCGGTCACCACGCCAATCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgcc                                                                                                       <  2:1021741/151‑1 (MQ=255)
                                          atCGCCCCTTTAAGGATGTCCGCCCAATTTGGGTTTCCCCCGATATCCGCCAGGCCATAGTTGAGGCCGGAAAAAATAATCCCCCCGGTCCCCCCCCCAACCACCGTCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCAcgccg                                                                                                      <  1:844975/151‑1 (MQ=255)
                                                  tttGATAATGTACTGCCAGTATGGTTTTCCGCCAATATACGTCAGGCCATATTTGATGACGGTAAAAAAATTTCCCCCGGCCCCCACGCCAACCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGc                                                                                               <  1:183905/150‑1 (MQ=255)
                                                    tgatgaTGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCTTATTTGATGACGGTAAAAATAATTCCCCCGGCCCCCACGCCAACCCCCGTCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGc                                                                                               <  2:386745/148‑1 (MQ=255)
                                                     gatgatGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATcc                                                                                           >  2:984914/1‑151 (MQ=255)
                                                      atgatgTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATcc                                                                                           >  2:124729/1‑150 (MQ=255)
                                                       taagttccTGCAAGTATGGGTTTACGCCGATATACGTCGGGCCATATTTGATGACGGTAAAAAAAATTCCCCCGGCCCCCACCCCAATCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                                                                          <  1:816749/143‑1 (MQ=255)
                                                         atgTACTGCCATTATGGTTTTACGCCAATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                                                                          <  2:1191692/148‑1 (MQ=255)
                                                          tgTACTGCCAGTAGGGTTTTACGCCGATATACGTCGGGCCATAGTTGATGACGGTAAAAATAATTCCCCCGGTCCCCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                                                                          <  1:1055829/147‑1 (MQ=255)
                                                           gTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTc                                                                                      >  1:1049190/1‑150 (MQ=255)
                                                            tACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcat                                                                                    >  2:1143972/1‑151 (MQ=255)
                                                               tGCCAGTAGGGGTTTACGCCGATATCCGTCAGGCCATAGTTGATGACTGTAAAAAAAATTCCCCCGCCCCCCCCGCCAATCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGc                                                                                        <  2:744016/144‑1 (MQ=255)
                                                                 ccAGTATGGTTTTACGCCGATATCCGCCAGGCCATAGTTGATGACGGTAAAAATAATTCCCCCGGTCCCCACCCCAATCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacata                                                                               <  1:552990/151‑1 (MQ=255)
                                                                 ccAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacat                                                                                >  2:916524/1‑150 (MQ=255)
                                                                  cAGTATGGGTTTCGGCCGATATACGTCAGGCAATGGTGGTTGCCGGAAAAAATAATCCCCCCGGTCACCGCGCCAATCACCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacata                                                                               <  2:1022175/150‑1 (MQ=255)
                                                                   aGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaa                                                                              >  1:214567/1‑150 (MQ=255)
                                                                     tATGGGTTTCCGCCATTATACGTCAGGCCTTATTTGAGGACGTTAAAAATAATTACCCCGTTCCCCCCGCCAACCACCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaaa                                                                             <  1:439229/149‑1 (MQ=255)
                                                                        gggTTTCCGCCAAAATACGCCGGCCCATACTTGAGGACGGTAAAAAAAATTCCCCCGGCCCCCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAg                                                                            <  2:92855/147‑1 (MQ=255)
                                                                        gggTTTACGCCAATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTCCCCCGGTCACCACGCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCg                                                                         <  1:584397/150‑1 (MQ=255)
                                                                         ggTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGcc                                                                          >  1:272525/1‑148 (MQ=255)
                                                                         ggTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCg                                                                         >  2:76438/1‑149 (MQ=255)
                                                                           tttACGCCGATATACGTAAGGCCATAGTTGATGACGGTAAAAATAATTCCCCCGGTCCCCACCCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAgg                                                                      <  2:1122280/150‑1 (MQ=255)
                                                                               cGCCGATATACGTCGGGCCATAGTGGATGACGGTAAAAAAAATTCCCCCGCTCCCCCCGCCAACCCCCGTCCCCCCCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgttg                                                                   <  1:191250/149‑1 (MQ=255)
                                                                                  cGATCTACGTGAGCCCATAGTGGAAAGCGTAAAAAATAATCCCCCCGGCCACAACGCCAACCCCCCCCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgttgt                                                                  <  1:1074782/147‑1 (MQ=255)
                                                                                  cGATATCCGTCGGCCCATAGTTATTACCGGAAAAAATATTTCCCCCGGTCCCCACCCCAACCCCCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGttgttg                                                                <  2:865024/149‑1 (MQ=255)
                                                                                   gATATACGCCAGGCCATAGTTGATGACGGTAAAAATAATTACCACGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggt                                                              >  1:703982/1‑150 (MQ=255)
                                                                                     tatcCGTCAGCCAATATTTGAGAGGGTAAAAAAAAATTCCCCCGGCCCCCACCCCAATCCCCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggtgg                                                            <  1:750332/150‑1 (MQ=255)
                                                                                        aCGTCAGTCAATAGTTAATACCAGTAAAAAAAATTCCCCCGCTGACCCCGCAAGTCCCCCCCCCCACCCCGCCCCTGAACGATACCCCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCag                                                         <  1:19127/150‑1 (MQ=255)
                                                                                         cGTAAGGCCATAGTTGATGCCGGTAAAAAAAATTCCCCCGGCCCCCCCGCCAATCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCaga                                                        <  1:197744/150‑1 (MQ=255)
                                                                                          gTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGc                                                      >  1:795148/1‑151 (MQ=255)
                                                                                              ggCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGa                                                   >  1:919595/1‑150 (MQ=255)
                                                                                              gcccATAGTTGATGACGGAAAAAATAATTACCCCGGTCCCCCCGCCAACCCCCGCCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGc                                                      <  1:321448/145‑1 (MQ=255)
                                                                                                caaTAGTTAATGACGGTAAAAATAATTACCCGGATCCCCACGCAAATCCCGGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATAc                                                <  2:382850/149‑1 (MQ=255)
                                                                                                  aTAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGt                                              >  1:929159/1‑151 (MQ=255)
                                                                                                   tAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACATCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTc                                             >  1:548661/1‑151 (MQ=255)
                                                                                                      ttGATGACGGTAAAAATAATTCCCCCGGTCCCCACGCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCg                                           <  2:718363/150‑1 (MQ=255)
                                                                                                         ataccGGTAAAAATAATTCCCCGGGCCCCCCCGCCAACCACGGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGc                                         <  1:398415/145‑1 (MQ=255)
                                                                                                           gaCGGTAAAAAAAATTCCCCCGGTCCCCACGCCATTCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCt                                        <  2:660525/148‑1 (MQ=255)
                                                                                                            aCGGTAAAAATAATTCCCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCt                                     <  2:516520/150‑1 (MQ=255)
                                                                                                              gatAAAAATAATTCCCCCGCCCCCCCCGCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTaa                                   <  1:676764/148‑1 (MQ=255)
                                                                                                                  aaaaaTAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATc                               >  1:1134737/1‑151 (MQ=255)
                                                                                                                   aaaaTATTTCCCCCGGCCCCCACGCCAACCCCCGTCCCAACCCCGCCGCTGAACGATGCACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAAc                                  <  2:12931/147‑1 (MQ=255)
                                                                                                                   aaaaTAATTACCCCGTTCCCCACGCCAATCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATcc                              <  1:266952/151‑1 (MQ=255)
                                                                                                                    aaaTAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATccc                             >  2:805864/1‑151 (MQ=255)
                                                                                                                     aataatTCCCCCGGTCCCCCCCCCAATCACCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCg                            <  2:648547/151‑1 (MQ=255)
                                                                                                                        aatTCCCCCGCCCACCACGCCAACCACCGTCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGcc                          <  2:562419/150‑1 (MQ=255)
                                                                                                                            aCCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAgg                     >  2:378392/1‑151 (MQ=255)
                                                                                                                                cggcccccacGCCAACCCCCGCCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATa                 <  1:200638/145‑1 (MQ=255)
                                                                                                                                    ccccccGCCAACCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGa               <  2:987029/144‑1 (MQ=255)
                                                                                                                                    ccccacCCCAACCCCCTCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAac             <  2:685513/149‑1 (MQ=255)
                                                                                                                                    caccacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAac             >  1:372466/1‑151 (MQ=255)
                                                                                                                                      ccacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAacac           >  2:217998/1‑151 (MQ=255)
                                                                                                                                      ccacGCCAACCCCCCCCCCCACACCGCCGCTGAACGATACACCGCCTGCCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGTCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAacac           <  1:300733/151‑1 (MQ=255)
                                                                                                                                       cacCCCAATCCCCTCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAacac           <  1:43294/150‑1 (MQ=255)
                                                                                                                                         cccaactccccGCCCCCCCGCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGcc        <  2:701662/142‑1 (MQ=255)
                                                                                                                                         cGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGcc        >  1:293630/1‑151 (MQ=255)
                                                                                                                                          ccactcccccGCCCCCGCACCGCCGCTGAACGATACACCGCCTACCTCCCACGTCGCGATGGCATCCAGCTCATACATACAGCCGCGTTTGTTGGTGTCAGCGCCGCTTCGTCCTGCTTCTAACTTCCCGCCAAAGTCATAGAACACg          <  1:353968/140‑1 (MQ=255)
                                                                                                                                             aacccccGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCa       <  1:886026/143‑1 (MQ=255)
                                                                                                                                             aaTCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCa       <  1:867134/148‑1 (MQ=255)
                                                                                                                                                cACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACaa  >  1:960746/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAA  >  NC_000913/2236753‑2237045

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