Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,748,110 C→A 7.5% L343M (CTG→ATG) ‡ lyxK → L‑xylulose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,748,1100CA 7.5% 233.9 / 8.5 90L343M (CTG→ATG) ‡lyxKL‑xylulose kinase
Reads supporting (aligned to +/- strand):  ref base C (40/43);  new base A (2/5);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

GTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACT  >  NC_000913/3747961‑3748259
                                                                                                                                                     |                                                                                                                                                     
gTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTcc                                                                                                                                                       <  2:1004659/150‑1 (MQ=255)
 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                                                                                                                                             <  2:756638/150‑1 (MQ=255)
          cGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcag                                                                                                                                             <  1:1100146/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                                                                            <  2:300590/150‑1 (MQ=255)
                                                                           tctcGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCAcc                                                                            <  1:804138/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)
                                                                                                                                                     |                                                                                                                                                     
GTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACT  >  NC_000913/3747961‑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: ‑