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

CAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAC‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCG  >  NC_000913/2348122‑2348609
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGCGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGTCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTCTGCAGGCAGTCGGTTTGGATCTGCCGTTTCAGAT‑GCGCTc                                                                                                                                                                                                                                               >  1:155011/1‑251 (MQ=255)
   aataggCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCGGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCaa                                                                                                                                                                                                                                            <  1:450974/246‑1 (MQ=255)
      atGCTGAATTTGCTGCTCTGCGGCGCGGACGATCCTGGGATGGCGGAAATTGCCGAAGAGTTTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTAAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAG‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCCCCAAccc                                                                                                                                                                                                                                         <  2:282562/251‑1 (MQ=255)
           gAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATc                                                                                                                                                                                                                                     <  1:311897/250‑1 (MQ=255)
               ctgctgCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGCTCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTg                                                                                                                                                                                                                                >  1:399462/1‑251 (MQ=255)
                tgctgcGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTgg                                                                                                                                                                                                                               >  1:344903/1‑251 (MQ=255)
                  ctgcGCAGCGGCGCGGACGCTCCTGAGCTGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCCGAC‑GCGCTCCCAACCGATCCCGTgg                                                                                                                                                                                                                               >  2:208993/1‑249 (MQ=255)
                        aGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAacac                                                                                                                                                                                                                       <  1:263918/251‑1 (MQ=255)
                           ggcgCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTg                                                                                                                                                                                                                    <  2:41751/251‑1 (MQ=255)
                                       ccTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTg                                                                                                                                                                                                         <  1:504167/250‑1 (MQ=255)
                                                gCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCa                                                                                                                                                                                               <  1:131057/251‑1 (MQ=255)
                                                     aaTTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAGAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTg                                                                                                                                                                                          >  1:322015/1‑251 (MQ=255)
                                                        tGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGCCGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTc                                                                                                                                                                                       >  1:374764/1‑251 (MQ=255)
                                                                agTGGAAGCTGGAGTGCTATGACCTGTTGGTTCAGGCAGCTCAAAAGGAGAAAGACTGGGCTGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCCCCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAggaa                                                                                                                                                                               <  2:323486/251‑1 (MQ=255)
                                                                     aaGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGCTTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAggaagtgg                                                                                                                                                                           >  1:78357/1‑250 (MQ=255)
                                                                                  aTGACCTGTTTGTTCAGGCAGCTCAAAAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCGCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTAt                                                                                                                                                             <  2:320406/251‑1 (MQ=255)
                                                                                             gttCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGc                                                                                                                                                   >  2:119145/1‑250 (MQ=255)
                                                                                                   gCAGCGCAACAGGAAAAAGACTGTGCGGATTTTCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑AGCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGGCAGTTCTTATCTGGTCGGGCAGATTGa                                                                                                                                            <  2:469036/251‑1 (MQ=255)
                                                                                                      gCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTc                                                                                                                                         <  1:200343/251‑1 (MQ=255)
                                                                                                                    aGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACcg                                                                                                                           <  2:430342/251‑1 (MQ=255)
                                                                                                                        tGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCACCCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTTCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCAGGCAGATTGACTCGGAAGTGGAAAccac                                                                                                                          >  2:338636/1‑246 (MQ=255)
                                                                                                                            gggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTg                                                                                                                   <  1:386452/250‑1 (MQ=255)
                                                                                                                                   tCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACt                                                                                                            >  1:191717/1‑251 (MQ=255)
                                                                                                                                         ccGCGACGGTGGGAGGATTGGTGTGAATAAAGGCAATCGCTGCGAGTTCGGTGAATACAACAACAAG‑TGCCGTATGCAGTCAGGCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGAATCGGAAGTGGACTCCGACGAGTTGAGTAACTTccgg                                                                                                       <  2:121643/250‑3 (MQ=255)
                                                                                                                                          cgcgACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGCTCTGACGTTCAAGAA‑GAGCTCCAACCCGATCCCGTGATTCACCACTTGGCTGGTTGCTGATAACGTGCAGAATGCTCAGCAGGAAATGGAAGTCAGTTCTTATCTGTACGGGCAGATTGCCTCGGAAGTGGAAACCAACGATTTGAGCAACGTCCAGct                                                                                                     >  2:376822/1‑251 (MQ=255)
                                                                                                                                            cgACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGctct                                                                                                   >  1:405312/1‑251 (MQ=255)
                                                                                                                                             gACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGCTCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTg                                                                                                  >  2:373355/1‑251 (MQ=255)
                                                                                                                                                ggTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATg                                                                                               >  1:186460/1‑251 (MQ=255)
                                                                                                                                                 gTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGCTCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGACGTGGACAACGACGATTTGAGAAACTTCCAGCTCTGATgg                                                                                              >  1:37327/1‑251 (MQ=255)
                                                                                                                                                         atTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAGCACCAAT‑ATCCGTAGGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGtt                                                                                      <  2:322015/251‑1 (MQ=255)
                                                                                                                                                                 gAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCt                                                                                 <  1:269407/248‑1 (MQ=255)
                                                                                                                                                                 gAATAAAGACATTCTATGCCAGTACGTTGACTACATCACCAAT‑ATCCGTCTGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAC‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACATGCAGGTTGCTCCGCAGGAAGAGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTgcg                                                                              >  1:396162/1‑251 (MQ=255)
                                                                                                                                                                    tAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGCCTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGCTGGCCCGCGTTACCCTGCGCAt                                                                           >  2:150545/1‑251 (MQ=255)
                                                                                                                                                                     aaaGACATTCCCGGCCTGTACGTTTAATACATC‑CTAATTATCCGTATGCAGGCCGTCGGTTTGGAGCTGCCGTTCCAGAC‑GCGCTCCAACCCGATCCCGTGGATCAACACGTGGCTGGTGTCTGATATCATGCAGGGTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTTATGGCCCGCGTTCCCCTGCGCATc                                                                          <  2:396162/251‑1 (MQ=255)
                                                                                                                                                                       aGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCAc                                                                        >  2:48939/1‑251 (MQ=255)
                                                                                                                                                                                  gCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGcacaca                                                               <  1:178165/249‑1 (MQ=255)
                                                                                                                                                                                   ccAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCAc                                                                        <  1:438648/239‑1 (MQ=255)
                                                                                                                                                                                   ccAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAAct                                                            >  1:351733/1‑251 (MQ=255)
                                                                                                                                                                                   ccAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCA‑ACAGCGCTCCAACCCGATCCCGTGGATCACAACTTGGCTGGTGACTGATAACGTGCGGGTTGCTCAGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAGCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATccc                                                                        >  2:438648/1‑237 (MQ=255)
                                                                                                                                                                                     agtacctggAATACCTCACCAAT‑AGCCTTATGCAGGCAGTCGGGTTGGATCGGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACGActgc                                                          <  2:481638/243‑1 (MQ=255)
                                                                                                                                                                                      gTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTTGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAActgct                                                         <  2:144496/251‑1 (MQ=255)
                                                                                                                                                                                           ttGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGcc                                                    >  1:484799/1‑251 (MQ=255)
                                                                                                                                                                                                 aCATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGAt                                               <  1:150545/250‑1 (MQ=255)
                                                                                                                                                                                                    tCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCAGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCATGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAac                                           >  1:218182/1‑251 (MQ=255)
                                                                                                                                                                                                    ccACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAac                                           <  1:444848/250‑1 (MQ=255)
                                                                                                                                                                                                     cACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAaca                                          <  1:48939/251‑1 (MQ=255)
                                                                                                                                                                                                     cACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAaca                                          <  1:373259/251‑1 (MQ=255)
                                                                                                                                                                                                                          cagTCGGTTTGGCTCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGAGGACACCGACGATTTGAGGAACTTCAAGCTCTGATGGCCCGCGTTACCCTGCGCATCCCTAGCACACAACAGCTGTGCCAGGATGCACCCCCTTCCATTCTGGCAgcgc                      >  2:260752/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTCGGTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTg                    >  1:306135/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggTTTGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGAGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCCCTGGCACACAACTGCTGTGCCCGGCTGAACACCCTTCCCTTCTGGCGGCGCTGGa                  >  2:451583/1‑250 (MQ=255)
                                                                                                                                                                                                                                 tttGGATCTGCCGTTCCAGAA‑GCGCTACAACACTATCCCGTAGATCAACCCTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATATGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAAt                >  1:409709/1‑250 (MQ=255)
                                                                                                                                                                                                                                   tGGATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATccc             >  2:449547/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ggATCTGCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAGCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCGAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATccc             <  1:70592/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCGTTCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAAtgt      <  2:405312/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tCCAGAA‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCg  <  1:373355/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATTGCCGAAGAGTGTAAGCAGGAGTGCTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAAT‑ATCCGTATGCAGGCAGTCGGTTTGGATCTGCCGTTCCAGAC‑GCGCTCCAACCCGATCCCGTGGATCAACACTTGGCTGGTGTCTGATAACGTGCAGGTTGCTCCGCAGGAAGTGGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCG  >  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: ‑