Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,022,725 C→T F470F (TTC→TTT xdhD → fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,022,7250CT98.6% 239.8 / ‑4.3 73F470F (TTC→TTTxdhDfused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (40/32);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTCCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAGCAACATTTG  >  NC_000913/3022483‑3022967
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
cagTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTacac                                                                                                                                                                                                                                              <  1:252046/251‑1 (MQ=255)
     cccGACTGAAACACAGATCTGCTCTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCGCCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTg                                                                                                                                                                                                                                         <  1:79810/251‑1 (MQ=255)
        gACTGAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTgaag                                                                                                                                                                                                                                      <  1:88507/251‑1 (MQ=255)
         aCTGAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAATAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTgaaga                                                                                                                                                                                                                                     <  1:291789/251‑1 (MQ=255)
            gAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGGGTGCGCCTGGGCCACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACCCCCGTggaga                                                                                                                                                                                                                                     >  2:88434/1‑248 (MQ=255)
            gAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGGTTTGGTTCCAAAAAGGACATCCTGCTGGAAGAAGGGTGCGCCTGGGCAACCTGCGTGACCGGGGGTCCGGGCCGGTGTCGCTACACCCGTGaagaagg                                                                                                                                                                                                                                  >  1:23924/1‑250 (MQ=255)
             aaaCACAGATCCTCTTTACGCGTATGGAGGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTCCCATGGCACTTACGCCGCCAGGTCGCGTGCGTCGTGGGCCTGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGTAAGAAGTGTGCGCCTGGGCTATCTGCGTGACCGGGCGTCCGGTGCTGTTTCGCTTCTCCCGAGAAGAagag                                                                                                                                                                                                                                 <  2:464980/251‑1 (MQ=255)
                 acaTATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTa                                                                                                                                                                                                                             <  2:176014/251‑1 (MQ=255)
                 acaTATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGATTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGGGCGCCTGGGCAACCTGCGGGAACGGGCGTCCGGTACTGTTTTGCTACACCCGTGAAGAagaggttt                                                                                                                                                                                                                             >  1:391067/1‑247 (MQ=255)
                      tCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAAAAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGGTTCGCTACACCCGTGAAGAAGAGTTTTTTGCt                                                                                                                                                                                                                        >  1:503817/1‑251 (MQ=255)
                           tttACTCGTATGGACGGCGATCGTCTGGTTCTCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAacac                                                                                                                                                                                                                   >  1:198659/1‑251 (MQ=255)
                            ttACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATTAAGTGCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCACCCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTTAAGAAGAGTTTATTGCTAacac                                                                                                                                                                                                                   <  2:247422/250‑1 (MQ=255)
                            ttACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGATTTTGTTCCAAACAGGACATCCTGCTGGAAGAAGTGCGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCTCTACACCCGTGAAGGAGAGTGTATTGCTAACAcc                                                                                                                                                                                                                  >  2:320852/1‑251 (MQ=255)
                              aCTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGGTCCCAACAGGACATCCTGGTGGGAGGAGTGTGGGCCTTGGCAACCTGAGTGACCGGGCGTGCGGTACGGTTTCGCTACACCCGGGAAGAAGAGGTTATTGCTAACAActc                                                                                                                                                                                                                >  2:129420/1‑251 (MQ=255)
                                 cGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCtctc                                                                                                                                                                                                              <  1:278471/249‑1 (MQ=255)
                                 cGTATGGACGGAGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACctct                                                                                                                                                                                                               <  1:326576/249‑1 (MQ=255)
                                   tATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGCCAACCTGCTGGGAGAAGGGTGCGCCTGGGCAACCTGCGGGACCGGGGGTCCGGGACTGTTTCGCTACACCCGTGACTAAGAGGTTATTGCTAACACCTCTCggc                                                                                                                                                                                                           >  1:336369/1‑249 (MQ=255)
                                   tATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTc                                                                                                                                                                                                           >  1:366948/1‑251 (MQ=255)
                                   tATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGGGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGGACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTc                                                                                                                                                                                                           >  2:202710/1‑251 (MQ=255)
                                      ggACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACg                                                                                                                                                                                                        <  1:88434/251‑1 (MQ=255)
                                      ggACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGAGGGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGGTTATTGCTacccactctcgcnacc                                                                                                                                                                                                        >  2:595142/1‑236 (MQ=255)
                                      ggACGGCGATCGTCGGGTTATCCACGCCTCCGCCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACTCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACg                                                                                                                                                                                                        <  2:10442/251‑1 (MQ=255)
                                            cGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGGTTTGGTTCCAAACAGGAAATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGAACCGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTtgcg                                                                                                                                                                                                   >  2:432986/1‑250 (MQ=255)
                                                cGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGt                                                                                                                                                                                              <  1:147464/251‑1 (MQ=255)
                                                  tCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGGGTGCGCCTGGGCCACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGGTTATTGCTAACACCTCTCGTCCCGTCCCGAAAGTc                                                                                                                                                                                             >  2:72188/1‑250 (MQ=255)
                                                            cACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACt                                                                                                                                                                                 >  1:203460/1‑251 (MQ=255)
                                                               gCCTCCACCCAGGTACCATGGAACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGATTTGGTTCCAAACAAGACATTCTGCTTGAAGAAGTGTGCGCCTGGGGAAACTGCGTGAACGGGCGTCCGGTACTGTTTCGCTACACACGTAAAGAAGAGCTTATTGCTAATACCCATCGTCCAGTCGCGAAAGTCACCGTCAAACTggg                                                                                                                                                                              >  2:562139/1‑251 (MQ=255)
                                                                ccTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGCTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTgggg                                                                                                                                                                             >  2:545919/1‑250 (MQ=255)
                                                                  tCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGCAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGc                                                                                                                                                                           >  1:631666/1‑251 (MQ=255)
                                                                        cAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAACGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGGGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGGTTATTGCTAACACCTCTCGTCACGTCGCGAAAGACACCGTCCAACTGGGAGCGanaa                                                                                                                                                                      >  2:594414/1‑250 (MQ=255)
                                                                         aGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGaaaaa                                                                                                                                                                     >  1:617064/1‑250 (MQ=255)
                                                                           gTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTACAAATGGGGGCGaaaaa                                                                                                                                                                     >  1:399741/1‑248 (MQ=255)
                                                                                aTGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTc                                                                                                                                                             <  2:366948/251‑1 (MQ=255)
                                                                                aTGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTc                                                                                                                                                             >  1:317730/1‑251 (MQ=255)
                                                                                aTGGCACTTAATCCTCCATGTCGCGCGCCGCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGTTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGTGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGATTTTATTGCTAACATCTCTCGTCACGTCGCGAAAGTCACCGTCAAACGTGGAGCTAAAAAAGATGGTc                                                                                                                                                             <  2:646273/251‑1 (MQ=255)
                                                                                           gcggccAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTTGCCGTCAAACTGGGAGCGAAAAAAGATGGTCGTCTTACGGCa                                                                                                                                                 <  2:618136/248‑1 (MQ=255)
                                                                                              gccAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGGACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTg                                                                                                                                              >  1:148503/1‑251 (MQ=255)
                                                                                               ccAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTg                                                                                                                                              >  1:388505/1‑250 (MQ=255)
                                                                                                              cGTGGGCATGCAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTAGGTCCAAACAGGACATACCGCTGGGAGAAGTGTGCGCCTGGGCAAACTGCGTGACCGGGCGTCCGGGACTGTTTCGCTCAACACGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCTCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATTGGCGCCTAAAGGCAGTGAAGATGGATTTcnccc                                                                                                                              >  2:598156/1‑247 (MQ=255)
                                                                                                              cGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCcg                                                                                                                                <  1:537907/249‑1 (MQ=255)
                                                                                                              cGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCACGCTGGAAGACGTGTGCGCCTGGGCAAACTGCGTGAACGGGCGTCCGGTAATGTTTCGCTACACCCGTGAAGAAGAGTTTCTTTCTAACACCGCTCGTCCAGTCGCGAAAATCACCGTCAAACTGGGAGTGGAAAAAGATGGGCGCCCTGTGGCAGTGCAGATGGATTTCcgtg                                                                                                                              >  2:435343/1‑249 (MQ=255)
                                                                                                                    cATGACACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGAACGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCCCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCCGTGAAGATGGATTTCCGCACCAac                                                                                                                         >  2:372294/1‑250 (MQ=255)
                                                                                                                       gAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTg                                                                                                                     <  2:148472/251‑1 (MQ=255)
                                                                                                                          aCAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGttt                                                                                                                                                                                                                              <  2:451509/145‑1 (MQ=255)
                                                                                                                          aCAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGttt                                                                                                                                                                                                                              >  1:451509/1‑145 (MQ=255)
                                                                                                                           cAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGGGAAGATGGATTTCCGCGCCAACACTGGccc                                                                                                                 >  2:432276/1‑251 (MQ=255)
                                                                                                                              cATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTa                                                                                                              <  2:503817/251‑1 (MQ=255)
                                                                                                                               aTAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTAc                                                                                                             <  1:511054/251‑1 (MQ=255)
                                                                                                                                             aTTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACcactcact                                                                                               >  1:259270/1‑251 (MQ=255)
                                                                                                                                                aaaGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACtcactcac                                                                                            <  1:392047/251‑1 (MQ=255)
                                                                                                                                                aaaGAGCGAGTTGGCGGCGGTTGTGGTTCCAAACAGTACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTTACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACtcactcac                                                                                            <  1:42209/251‑1 (MQ=255)
                                                                                                                                                    agCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCaccgta                                                                                        <  1:320852/251‑1 (MQ=255)
                                                                                                                                                             gcggcgGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTaat                                                                               <  2:617064/251‑2 (MQ=255)
                                                                                                                                                               ggcggTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTaatgg                                                                             <  1:482646/251‑4 (MQ=255)
                                                                                                                                                                 cggTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGAc                                                                           <  1:371954/251‑1 (MQ=255)
                                                                                                                                                                         gTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTg                                                                   >  2:92745/1‑251 (MQ=255)
                                                                                                                                                                                  aGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCg                                                          <  1:72188/251‑1 (MQ=255)
                                                                                                                                                                                     aCATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGtt                                                        >  1:473170/1‑250 (MQ=255)
                                                                                                                                                                                        tCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTatat                                                    >  1:90188/1‑251 (MQ=255)
                                                                                                                                                                                         ccTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTatat                                                    <  2:317730/250‑1 (MQ=255)
                                                                                                                                                                                         ccTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATc                                                   >  2:546942/1‑251 (MQ=255)
                                                                                                                                                                                             ctgGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGCAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTg                                               >  2:450999/1‑251 (MQ=255)
                                                                                                                                                                                                gaagaaGTGTGCGCCTGGGAAACCTTCGTGACCGGGCGTCCGGTACTGTTTCGCTAAACCCGTGAAGCAGAGTTTAGTGCTAACACCTCTCGTCACGTCGCGACAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGAAGTGAAGATGGATTTCCGCGCCAACACTAGCCCTTACGGCAACCACGCACTCACCGTACCGTGTAATGGACCGGCGCTGGCGCTGCCGTTATATCCGTGCg                                             >  1:584961/1‑250 (MQ=255)
                                                                                                                                                                                                 aagaagTGTGCGCCTGGGCATCCTGCGTGACCGGGCGTCCGTTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTAACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACTGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATTGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGag                                           <  2:103258/251‑2 (MQ=255)
                                                                                                                                                                                                    aagTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAAc                                        <  1:441362/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccTGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTcacc                      >  2:571152/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGCGTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTcaccac                    >  2:547100/1‑251 (MQ=255)
                                                                                                                                                                                                                           gTGACCGGGCGTCCGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACct                  >  2:341469/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gggCGTCCGGTACTGTTTCGCTACGCCCGGGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTTACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAg           <  2:253073/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAGCAACa      <  1:594414/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAGCAACAtt    >  2:211051/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAGCAACAttt   <  1:212978/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTACTGTTTCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAATGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAGCAACATTTg  >  2:548662/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
CAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTCCGCTACACCCGTGAAGAAGAGTTTATTGCTAACACCTCTCGTCACGTCGCGAAAGTCACCGTCAAACTGGGAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCGCTGCCGTTATATCCGTGCGATAACGTCGATTTCCAGGTCACCACCTACTACAGCAACATTTG  >  NC_000913/3022483‑3022967

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