Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,348,366 C→A T328K (ACG→AAG)  nrdB → ribonucleoside‑diphosphate reductase 1 subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,348,3660CA92.6% 166.9 / ‑1.8 55T328K (ACG→AAG) nrdBribonucleoside‑diphosphate reductase 1 subunit beta
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (23/27);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCG  >  NC_000913/2348122‑2348609
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGCGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGTCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTCTGCAGGCAGTCGGTTTGGATCTGCCGTTTCAGATGCGCTc                                                                                                                                                                                                                                               >  1:155011/1‑251 (MQ=255)
   aataggCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCGGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCaa                                                                                                                                                                                                                                            <  1:450974/246‑1 (MQ=255)
      atGCTGAATTTGCTGCTCTGCGGCGCGGACGATCCTGGGATGGCGGAAATTGCCGAAGAGTTTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTAAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAGATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCCCCAAccc                                                                                                                                                                                                                                         <  2:282562/251‑1 (MQ=255)
           gAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATc                                                                                                                                                                                                                                     <  1:311897/250‑1 (MQ=255)
               ctgctgCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGCTCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTg                                                                                                                                                                                                                                >  1:399462/1‑251 (MQ=255)
                tgctgcGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTgg                                                                                                                                                                                                                               >  1:344903/1‑251 (MQ=255)
                  ctgcGCAGCGGCGCGGACGCTCCTGAGCTGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCCGACGCGCTCCCAACCGATCCCGTgg                                                                                                                                                                                                                               >  2:208993/1‑249 (MQ=255)
                        aGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAacac                                                                                                                                                                                                                       <  1:263918/251‑1 (MQ=255)
                           ggcgCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTg                                                                                                                                                                                                                    <  2:41751/251‑1 (MQ=255)
                                       ccTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTg                                                                                                                                                                                                         <  1:504167/250‑1 (MQ=255)
                                                gCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCa                                                                                                                                                                                               <  1:131057/251‑1 (MQ=255)
                                                     aaTTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAGAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTg                                                                                                                                                                                          >  1:322015/1‑251 (MQ=255)
                                                        tGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGCCGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTc                                                                                                                                                                                       >  1:374764/1‑251 (MQ=255)
                                                                agTGGAAGCTGGAGTGCTATGACCTGTTGGTTCAGGCAGCTCAAAAGGAGAAAGACTGGGCTGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCCCCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAggaa                                                                                                                                                                               <  2:323486/251‑1 (MQ=255)
                                                                     aaGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGCTTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAggaagtgg                                                                                                                                                                           >  1:78357/1‑250 (MQ=255)
                                                                                  aTGACCTGTTTGTTCAGGCAGCTCAAAAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCGCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTAt                                                                                                                                                             <  2:320406/251‑1 (MQ=255)
                                                                                             gttCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGc                                                                                                                                                   >  2:119145/1‑250 (MQ=255)
                                                                                                   gCAGCGCAACAGGAAAAAGACTGTGCGGATTTTCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATAGCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGGCAGTTCTTATCTGGTCGGGCAGATTGa                                                                                                                                            <  2:469036/251‑1 (MQ=255)
                                                                                                      gCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTc                                                                                                                                         <  1:200343/251‑1 (MQ=255)
                                                                                                                    aGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACcg                                                                                                                           <  2:430342/251‑1 (MQ=255)
                                                                                                                        tGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCACCCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTTCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCAGGCAGATTGACTCGGAAGTGGAAAccac                                                                                                                          >  2:338636/1‑246 (MQ=255)
                                                                                                                            gggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTg                                                                                                                   <  1:386452/250‑1 (MQ=255)
                                                                                                                                   tCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACt                                                                                                            >  1:191717/1‑251 (MQ=255)
                                                                                                                                         ccGCGACGGTGGGAGGATTGGTGTGAATAAAGGCAATCGCTGCGAGTTCGGTGAATACAACAACAAGTGCCGTATGCAGTCAGGCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGAATCGGAAGTGGACTCCGACGAGTTGAGTAACTTccgg                                                                                                       <  2:121643/250‑3 (MQ=255)
                                                                                                                                          cgcgACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGCTCTGACGTTCAAGAAGAGCTCCAACCCGATCCCGTGATTCACCACTTGGCTGGTTGCTGATAACGTGCAGAATGCTCAGCAGGAAATGGAAGTCAGTTCTTATCTGTACGGGCAGATTGCCTCGGAAGTGGAAACCAACGATTTGAGCAACGTCCAGct                                                                                                     >  2:376822/1‑251 (MQ=255)
                                                                                                                                            cgACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGctct                                                                                                   >  1:405312/1‑251 (MQ=255)
                                                                                                                                             gACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGCTCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTg                                                                                                  >  2:373355/1‑251 (MQ=255)
                                                                                                                                                ggTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATg                                                                                               >  1:186460/1‑251 (MQ=255)
                                                                                                                                                 gTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGCTCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGACGTGGACAACGACGATTTGAGAAACTTCCAGCTCTGATgg                                                                                              >  1:37327/1‑251 (MQ=255)
                                                                                                                                                         atTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAGCACCAATATCCGTAGGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGtt                                                                                      <  2:322015/251‑1 (MQ=255)
                                                                                                                                                                 gAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCt                                                                                 <  1:269407/248‑1 (MQ=255)
                                                                                                                                                                 gAATAAAGACATTCTATGCCAGTACGTTGACTACATCACCAATATCCGTCTGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACATGCAGGTTGCTCCGCAGGAAGAGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTgcg                                                                              >  1:396162/1‑251 (MQ=255)
                                                                                                                                                                    tAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGCCTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGCTGGCCCGCGTTACCCTGCGCAt                                                                           >  2:150545/1‑251 (MQ=255)
                                                                                                                                                                     aaaGACATTCCCGGCCTGTACGTTTAATACATCCTAATTATCCGTATGCAGGCCGTCGGTTTGGAGCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACGTGGCTGGTGTCTGATATCATGCAGGGTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTTATGGCCCGCGTTCCCCTGCGCATc                                                                          <  2:396162/251‑1 (MQ=255)
                                                                                                                                                                       aGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCAc                                                                        >  2:48939/1‑251 (MQ=255)
                                                                                                                                                                                  gCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGcacaca                                                               <  1:178165/249‑1 (MQ=255)
                                                                                                                                                                                   ccAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCAc                                                                        <  1:438648/239‑1 (MQ=255)
                                                                                                                                                                                   ccAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAAct                                                            >  1:351733/1‑251 (MQ=255)
                                                                                                                                                                                   ccAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAACAGCGCTCCAACCCGATCCCGTGGATCACAACTTGGCTGGTGACTGATAACGTGCGGGTTGCTCAGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAGCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATccc                                                                        >  2:438648/1‑237 (MQ=255)
                                                                                                                                                                                     agtacctggAATACCTCACCAATAGCCTTATGCAGGCAGTCGGGTTGGATCGGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACGActgc                                                          <  2:481638/243‑1 (MQ=255)
                                                                                                                                                                                      gTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTTGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAActgct                                                         <  2:144496/251‑1 (MQ=255)
                                                                                                                                                                                           ttGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGcc                                                    >  1:484799/1‑251 (MQ=255)
                                                                                                                                                                                                 aCATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGAt                                               <  1:150545/250‑1 (MQ=255)
                                                                                                                                                                                                    tCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCAGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCATGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAac                                           >  1:218182/1‑251 (MQ=255)
                                                                                                                                                                                                    ccACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAac                                           <  1:444848/250‑1 (MQ=255)
                                                                                                                                                                                                     cACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAaca                                          <  1:48939/251‑1 (MQ=255)
                                                                                                                                                                                                     cACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAaca                                          <  1:373259/251‑1 (MQ=255)
                                                                                                                                                                                                                          cagTCGGTTTGGCTCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGAGGACACCGACGATTTGAGGAACTTCAAGCTCTGATGGCCCGCGTTACCCTGCGCATCCCTAGCACACAACAGCTGTGCCAGGATGCACCCCCTTCCATTCTGGCAgcgc                      >  2:260752/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTCGGTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTg                    >  1:306135/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggTTTGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGAGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCCCTGGCACACAACTGCTGTGCCCGGCTGAACACCCTTCCCTTCTGGCGGCGCTGGa                  >  2:451583/1‑250 (MQ=255)
                                                                                                                                                                                                                                 tttGGATCTGCCGTTCCAGAAGCGCTACAACACTATCCCGTAGATCAACCCTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATATGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAAt                >  1:409709/1‑250 (MQ=255)
                                                                                                                                                                                                                                   tGGATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATccc             >  2:449547/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ggATCTGCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAGCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCGAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATccc             <  1:70592/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCGTTCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAAtgt      <  2:405312/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tCCAGAAGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCg  <  1:373355/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGACGCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCG  >  NC_000913/2348122‑2348609

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