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

CGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAA  >  NC_000913/2236753‑2237045
                                                                                                                                                     |                                                                                                                                                
cgcgcgTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTcc                                                                                                                                                <  1:287704/151‑1 (MQ=255)
           cAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACAccgccg                                                                                                                                     >  2:845008/1‑151 (MQ=255)
           cAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACAccgcc                                                                                                                                      <  2:148056/150‑1 (MQ=255)
              tGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTg                                                                                                                                  >  2:698217/1‑151 (MQ=255)
                 aTCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAAc                                                                                                                               >  1:783443/1‑151 (MQ=255)
                   ccAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACg                                                                                                                              >  1:860666/1‑150 (MQ=255)
                     aGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGAt                                                                                                                            >  2:38316/1‑150 (MQ=255)
                      gcgtTACGGCGAAGATAATAATCGCCCCTTTGATGATTTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGAt                                                                                                                            <  1:435165/149‑1 (MQ=255)
                               cGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTAc                                                                                                                 >  1:328545/1‑151 (MQ=255)
                                 aaGATAATACTCGCCCCTTTGATGTGGTACTGCCAGTATGGTTTTACGCCAATATACGTCAGCCCATAGTTGATGAGGGT‑AAAAATAATTCCCCCGTCCCCCCCCCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTAc                                                                                                                 <  2:624493/149‑1 (MQ=255)
                                 aaGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTAcc                                                                                                                >  1:915999/1‑150 (MQ=255)
                                   gATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACcacg                                                                                                             >  2:43299/1‑151 (MQ=255)
                                    ataataATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACcacgc                                                                                                            >  2:1091609/1‑151 (MQ=255)
                                     taataaTCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACcacgc                                                                                                            >  2:81139/1‑150 (MQ=255)
                                         aatCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTCCCCCGGTCACCACGCCAATCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgcc                                                                                                       <  2:1021741/151‑1 (MQ=255)
                                          atCGCCCCTTTAAGGATGTCCGCCCAATTTGGGTTTCCCCCGATATCCGCCAGGCCATAGTTGAGGCCGGA‑AAAAATAATCCCCCCGGTCCCCCCCCCAACCACCGTCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCAcgccg                                                                                                      <  1:844975/151‑1 (MQ=255)
                                                  tttGATAATGTACTGCCAGTATGGTTTTCCGCCAATATACGTCAGGCCATATTTGATGACGGT‑AAAAAAATTTCCCCCGGCCCCCACGCCAACCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGc                                                                                               <  1:183905/150‑1 (MQ=255)
                                                    tgatgaTGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCTTATTTGATGACGGT‑AAAAATAATTCCCCCGGCCCCCACGCCAACCCCCGTCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGc                                                                                               <  2:386745/148‑1 (MQ=255)
                                                     gatgatGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATcc                                                                                           >  2:984914/1‑151 (MQ=255)
                                                      atgatgTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATcc                                                                                           >  2:124729/1‑150 (MQ=255)
                                                       taagttccTGCAAGTATGGGTTTACGCCGATATACGTCGGGCCATATTTGATGACGGT‑AAAAAAAATTCCCCCGGCCCCCACCCCAATCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                                                                          <  1:816749/143‑1 (MQ=255)
                                                         atgTACTGCCATTATGGTTTTACGCCAATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                                                                          <  2:1191692/148‑1 (MQ=255)
                                                          tgTACTGCCAGTAGGGTTTTACGCCGATATACGTCGGGCCATAGTTGATGACGGT‑AAAAATAATTCCCCCGGTCCCCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                                                                          <  1:1055829/147‑1 (MQ=255)
                                                           gTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTc                                                                                      >  1:1049190/1‑150 (MQ=255)
                                                            tACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcat                                                                                    >  2:1143972/1‑151 (MQ=255)
                                                               tGCCAGTAGGGGTTTACGCCGATATCCGTCAGGCCATAGTTGATGACTGT‑AAAAAAAATTCCCCCGCCCCCCCCGCCAATCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGc                                                                                        <  2:744016/144‑1 (MQ=255)
                                                                 ccAGTATGGTTTTACGCCGATATCCGCCAGGCCATAGTTGATGACGGT‑AAAAATAATTCCCCCGGTCCCCACCCCAATCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacata                                                                               <  1:552990/151‑1 (MQ=255)
                                                                 ccAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacat                                                                                >  2:916524/1‑150 (MQ=255)
                                                                  cAGTATGGGTTTCGGCCGATATACGTCAGGCAATGGTGGTTGCCGGA‑AAAAATAATCCCCCCGGTCACCGCGCCAATCACCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacata                                                                               <  2:1022175/150‑1 (MQ=255)
                                                                   aGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaa                                                                              >  1:214567/1‑150 (MQ=255)
                                                                     tATGGGTTTCCGCCATTATACGTCAGGCCTTATTTGAGGACGTT‑AAAAATAATTACCCCGTTCCCCCCGCCAACCACCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaaa                                                                             <  1:439229/149‑1 (MQ=255)
                                                                        gggTTTCCGCCAAAATACGCCGGCCCATACTTGAGGACGGT‑AAAAAAAATTCCCCCGGCCCCCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAg                                                                            <  2:92855/147‑1 (MQ=255)
                                                                        gggTTTACGCCAATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTCCCCCGGTCACCACGCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCg                                                                         <  1:584397/150‑1 (MQ=255)
                                                                         ggTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGcc                                                                          >  1:272525/1‑148 (MQ=255)
                                                                         ggTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCg                                                                         >  2:76438/1‑149 (MQ=255)
                                                                           tttACGCCGATATACGTAAGGCCATAGTTGATGACGGT‑AAAAATAATTCCCCCGGTCCCCACCCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAgg                                                                      <  2:1122280/150‑1 (MQ=255)
                                                                               cGCCGATATACGTCGGGCCATAGTGGATGACGGT‑AAAAAAAATTCCCCCGCTCCCCCCGCCAACCCCCGTCCCCCCCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgttg                                                                   <  1:191250/149‑1 (MQ=255)
                                                                                  cGATCTACGTGAGCCCATAGTGGAAAGCGTA‑AAAAATAATCCCCCCGGCCACAACGCCAACCCCCCCCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgttgt                                                                  <  1:1074782/147‑1 (MQ=255)
                                                                                  cGATATCCGTCGGCCCATAGTTATTACCGGA‑AAAAATATTTCCCCCGGTCCCCACCCCAACCCCCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGttgttg                                                                <  2:865024/149‑1 (MQ=255)
                                                                                   gATATACGCCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCACGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggt                                                              >  1:703982/1‑150 (MQ=255)
                                                                                     tatcCGTCAGCCAATATTTGA‑GAGGGTAAAAAAAAATTCCCCCGGCCCCCACCCCAATCCCCGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggtgg                                                            <  1:750332/150‑1 (MQ=255)
                                                                                        aCGTCAGTCAATAGTTAATACCAGT‑AAAAAAAATTCCCCCGCTGACCCCGCAAGTCCCCCCCCCCACCCCGCCCCTGAACGATACCCCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCag                                                         <  1:19127/150‑1 (MQ=255)
                                                                                         cGTAAGGCCATAGTTGATGCCGGT‑AAAAAAAATTCCCCCGGCCCCCCCGCCAATCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCaga                                                        <  1:197744/150‑1 (MQ=255)
                                                                                          gTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGc                                                      >  1:795148/1‑151 (MQ=255)
                                                                                              ggCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGa                                                   >  1:919595/1‑150 (MQ=255)
                                                                                              gcccATAGTTGATGACGGA‑AAAAATAATTACCCCGGTCCCCCCGCCAACCCCCGCCCCCACCCCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGc                                                      <  1:321448/145‑1 (MQ=255)
                                                                                                caaTAGTTAATGACGGT‑AAAAATAATTACCCGGATCCCCACGCAAATCCCGGCCCCCCCACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATAc                                                <  2:382850/149‑1 (MQ=255)
                                                                                                  aTAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGt                                              >  1:929159/1‑151 (MQ=255)
                                                                                                   tAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACATCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTc                                             >  1:548661/1‑151 (MQ=255)
                                                                                                      ttGATGACGGT‑AAAAATAATTCCCCCGGTCCCCACGCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCg                                           <  2:718363/150‑1 (MQ=255)
                                                                                                         ataccGGT‑AAAAATAATTCCCCGGGCCCCCCCGCCAACCACGGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGc                                         <  1:398415/145‑1 (MQ=255)
                                                                                                           gaCGGT‑AAAAAAAATTCCCCCGGTCCCCACGCCATTCCCCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCt                                        <  2:660525/148‑1 (MQ=255)
                                                                                                            aCGGT‑AAAAATAATTCCCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCt                                     <  2:516520/150‑1 (MQ=255)
                                                                                                              gat‑AAAAATAATTCCCCCGCCCCCCCCGCCAACCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTaa                                   <  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)
                                                                                                                                                     |                                                                                                                                                
CGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTGATGATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGT‑AAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAA  >  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: ‑