Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,748,111 T→G 7.7% L343R (CTG→CGG) ‡ lyxK → L‑xylulose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,748,1110TG 7.7% 248.5 / 10.3 91L343R (CTG→CGG) ‡lyxKL‑xylulose kinase
Reads supporting (aligned to +/- strand):  ref base T (40/44);  new base G (2/5);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

TGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACT  >  NC_000913/3747962‑3748259
                                                                                                                                                     |                                                                                                                                                    
tGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCt                                                                                                                                                      >  1:305011/1‑150 (MQ=255)
  aaattATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCATGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGt                                                                                                                                                    <  2:559712/147‑1 (MQ=255)
  aaTTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGt                                                                                                                                                    >  1:406871/1‑150 (MQ=255)
       aTCGTTCACGTAGCCAGCCCTACCTCTTCCGGCATCCTCGATTGGTTTACCGCACAGTGGGGAGAATTCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgc                                                                                                                                               <  1:1134317/150‑1 (MQ=255)
         cGTTCACGACGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATGTCGTTTGATGAGATCAATCAGGACGTTGCCAGCTGGGCGCGGGCAGGGGGCGGGCGCGGTGTCCTGCCGTTACTGTCCGgcag                                                                                                                                             >  2:60933/1‑150 (MQ=255)
         cGTTCACGAAGCCAGCCCTCCCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcag                                                                                                                                             <  1:859876/150‑1 (MQ=255)
         cGTTCACGAAGCCAGCCCTACCTCTTCCGTCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcag                                                                                                                                             <  1:864041/150‑1 (MQ=255)
         cGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcag                                                                                                                                             >  2:422798/1‑150 (MQ=255)
         cGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcag                                                                                                                                             <  1:1100146/150‑1 (MQ=255)
         cGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcag                                                                                                                                             <  2:756638/150‑1 (MQ=255)
          gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagc                                                                                                                                            <  2:766500/150‑1 (MQ=255)
              aCGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGGCAGGaaaa                                                                                                                                        <  2:752889/150‑3 (MQ=255)
                gAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGcc                                                                                                                                      <  1:5674/150‑1 (MQ=255)
                  agccagccCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCgg                                                                                                                                    <  1:298953/150‑1 (MQ=255)
                    ccagccCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGAc                                                                                                                                  >  2:780457/1‑150 (MQ=255)
                    ccagccCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGAc                                                                                                                                  <  1:1012797/150‑1 (MQ=255)
                       gccCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCg                                                                                                                               >  1:459188/1‑150 (MQ=255)
                         ccTACCTCTTCCGGCAAACTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCCGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCgag                                                                                                                             >  1:11087/1‑150 (MQ=255)
                          cTACCTCTTCCGGCAACCTCGAATGGGTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCgaga                                                                                                                            <  2:960551/150‑1 (MQ=255)
                              ctctTCCGGCAACCTCGAATGG‑TTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGAcc                                                                                                                       <  2:1096977/150‑1 (MQ=255)
                               tctTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCATTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGAcc                                                                                                                       <  1:845075/150‑1 (MQ=255)
                                ctTCCGGCAACCTCGATTGGTTTACCGCCCAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGACGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCa                                                                                                                      >  2:1284706/1‑150 (MQ=255)
                                  tCCGGCAACCTCGATTGGTTTACCGCACAGTGGGGAGAAATCTCGCTTGATGAGATCAATCAGGCCGTTGCCAGCATGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGGAACGACGGACTCGAGATGACCAGt                                                                                                                    >  1:275855/1‑150 (MQ=255)
                                  tCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTCCGGCAGCAACGCCGGACTCGAGATGACCAGt                                                                                                                    >  2:110294/1‑150 (MQ=255)
                                   ccGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTg                                                                                                                   >  2:150908/1‑150 (MQ=255)
                                    cGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTgg                                                                                                                  >  1:271132/1‑150 (MQ=255)
                                     ggCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCGACGCCGGACTCGAGATGACCAGTGGt                                                                                                                 >  2:1183149/1‑150 (MQ=255)
                                     ggCAACCTCGAATGGTTTACCGAACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGt                                                                                                                 >  2:272334/1‑150 (MQ=255)
                                        aaCCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTc                                                                                                              <  2:957548/150‑1 (MQ=255)
                                          ccTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCGTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGCCTCGAGATGACCAGTGGTTTCTa                                                                                                            >  2:248952/1‑150 (MQ=255)
                                           cTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTAc                                                                                                           >  1:417259/1‑150 (MQ=255)
                                            tCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACg                                                                                                          <  1:1025431/150‑1 (MQ=255)
                                             cGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGAGCTCTTTTGCATGCCGTGCCTGTACGGCCGCACCGCCGGCGTCGAGATGACAGGTGGTTAATACgg                                                                                                         >  2:1079279/1‑150 (MQ=255)
                                              gAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACggg                                                                                                        >  2:923340/1‑150 (MQ=255)
                                               aaTGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGa                                                                                                       <  1:408474/150‑1 (MQ=255)
                                               aaTGGTTGACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGa                                                                                                       <  1:835374/150‑1 (MQ=255)
                                                  ggTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGc                                                                                                    <  1:31696/150‑1 (MQ=255)
                                                  ggTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGc                                                                                                    <  1:977481/150‑1 (MQ=255)
                                                   gTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCAGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGTCAGCAACGCCGGAATCGAGATGACCAGTGGTTTCAACGGGATGCa                                                                                                   >  2:367463/1‑150 (MQ=255)
                                                         cGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCAt                                                                                             >  2:760286/1‑150 (MQ=255)
                                                          gCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCAtt                                                                                            <  2:181763/150‑1 (MQ=255)
                                                           cacaGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTc                                                                                           <  1:736442/150‑1 (MQ=255)
                                                           cacaGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTc                                                                                           >  1:605421/1‑150 (MQ=255)
                                                              aGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTcaca                                                                                        <  2:994494/150‑1 (MQ=255)
                                                               gTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCATGCCTAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGACGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCAGTCccac                                                                                       >  1:852449/1‑150 (MQ=255)
                                                                   ggAGATATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTCTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGCGATGCAGGCCAGTCACACCcgc                                                                                   >  2:719021/1‑150 (MQ=255)
                                                                   ggAGAAATCTCGTTTGATGAGATCAATCAGGCGGTATCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCTTACCTGTACGGCAGCAACGCCGGACTCGAGATGCCGAGTGGTTTCTCCGGGATGCAGGCCATTCACACCcgc                                                                                   >  2:520388/1‑150 (MQ=255)
                                                                   ggAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCTTTCCTGTACGGCAGCAACGCCGGACTCGACATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCcgg                                                                                   >  2:1211419/1‑149 (MQ=255)
                                                                        aaTCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCAGAAGGATGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACCCCcgcgcgaa                                                                              >  1:1210253/1‑148 (MQ=255)
                                                                         aTCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCAc                                                                             <  2:965459/150‑1 (MQ=255)
                                                                          tctcGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCAGACCCGCGCGCAcc                                                                            <  2:305011/150‑1 (MQ=255)
                                                                          tctcGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCAcc                                                                            <  1:804138/150‑1 (MQ=255)
                                                                          tctcGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCAcc                                                                            <  2:300590/150‑1 (MQ=255)
                                                                              gTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGtt                                                                        <  2:892328/150‑1 (MQ=255)
                                                                               gttGCTGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTg                                                                       <  2:25696/149‑1 (MQ=255)
                                                                                    ggagaTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGc                                                                  <  1:1241405/149‑1 (MQ=255)
                                                                                      atatcactctGGCCGTTTCCATCTTGCCGAAGGCTTGGGGGGATCGCTTTTTCCTGCCTTTCCTGTAAGGAAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTTCGGGATGCAGGCCATTCACACCCGCGCGCCCCTGTTGCAGGCCa                                                                <  1:619864/148‑1 (MQ=255)
                                                                                      agaTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCa                                                                >  2:514060/1‑150 (MQ=255)
                                                                                       gaTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCAt                                                               <  2:962255/150‑1 (MQ=255)
                                                                                        atcaatcaGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATc                                                              <  2:361393/150‑1 (MQ=255)
                                                                                          caatgaGGCCGTTGCCTGCTTGGCTAAGGCTGGGGGCTATCTCATTATCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCAGACCCGCGCGCACCTGTTGCAGGCCATCTa                                                            <  1:579988/150‑1 (MQ=255)
                                                                                            atcaGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGGACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATg                                                          <  2:129704/150‑1 (MQ=255)
                                                                                               aGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAg                                                       <  1:1034275/150‑1 (MQ=255)
                                                                                                ggCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGATGCAGGCCATCTATGAAgg                                                      <  1:272334/150‑1 (MQ=255)
                                                                                                ggCCGGTGCCAGCTGGCCGAAGGCGGGGGGCGAGCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAgg                                                      <  1:58220/150‑1 (MQ=255)
                                                                                                   cGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGggt                                                   >  2:19632/1‑150 (MQ=255)
                                                                                                    gTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCAGACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGggtg                                                  <  1:923340/150‑1 (MQ=255)
                                                                                                     ttGCCAGCTTGCCGAAGGCTGGGGGCGAGCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGggtgg                                                 <  2:490629/150‑1 (MQ=255)
                                                                                                       gCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGgtggtg                                               >  1:1220050/1‑150 (MQ=255)
                                                                                                        ccAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGgtgt                                              >  1:980004/1‑150 (MQ=255)
                                                                                                         cAGCTGGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGtt                                             <  1:444885/150‑1 (MQ=255)
                                                                                                            cTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAg                                          <  2:616817/150‑1 (MQ=255)
                                                                                                            cTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAg                                          >  1:105041/1‑150 (MQ=255)
                                                                                                              tGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGCCCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCAGCTATGAAGGGGTGGTGTTCAGcc                                        >  1:975064/1‑150 (MQ=255)
                                                                                                              tGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGcc                                        >  1:1046861/1‑150 (MQ=255)
                                                                                                                ccgccggCGGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCCGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCat                                      <  2:217043/145‑1 (MQ=255)
                                                                                                                ccGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCCCACCCGCGCGCACCTGTTGCAGGCCCTCTATGAAGGGGTGGTGTTCAGCCat                                      >  1:418902/1‑150 (MQ=255)
                                                                                                                   aaGTCTGGGGGCGATCTCTTTTTCCTGCCTTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCTCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATg                                   <  1:482493/150‑1 (MQ=255)
                                                                                                                   aaGGCTGGGTGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCTGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATg                                   <  2:681004/150‑1 (MQ=255)
                                                                                                                    aGGCTGGGGGCGATCTCTTTTTCCTGCCGTACCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCTCACCAGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGACATATGa                                  >  1:965040/1‑150 (MQ=255)
                                                                                                                     ggCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGAc                                 <  1:482837/150‑1 (MQ=255)
                                                                                                                       cTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGAccc                               >  1:832090/1‑150 (MQ=255)
                                                                                                                        tGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGGGGTGTTCAGCCATATGACCCa                              <  2:406085/150‑1 (MQ=255)
                                                                                                                          ggggCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCAcc                            <  2:1096748/150‑1 (MQ=255)
                                                                                                                              cGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCaa                        >  2:1290816/1‑150 (MQ=255)
                                                                                                                                    ctTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAAt                  >  2:1203046/1‑150 (MQ=255)
                                                                                                                                          gcTGCCGTACAGGAACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCCCCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGCCCCACCTCAACCGAATGCGCGa            >  2:896203/2‑150 (MQ=255)
                                                                                                                                          gcTGCCGTACAGGAACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGa            >  2:739484/2‑150 (MQ=255)
                                                                                                                                            tGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAAc          >  2:1211884/1‑150 (MQ=255)
                                                                                                                                             gCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACg         >  1:1167042/1‑150 (MQ=255)
                                                                                                                                                   ccTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTAc   <  2:406871/150‑1 (MQ=255)
                                                                                                                                                    cTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACt  >  1:938553/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACT  >  NC_000913/3747962‑3748259

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