Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,375,737 G→A A259T (GCC→ACC)  ycjQ → D‑guloside 3‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,375,7370GA94.5% 156.4 / 0.3 58A259T (GCC→ACC) ycjQD‑guloside 3‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (21/31);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

AACGCGCTGGCGCTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGACCA  >  NC_000913/1375495‑1375975
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
aaCGCGCTGGCGCTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGACACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCCCCATCTCCTATGTCGCGTTTACCAAGCCg                                                                                                                                                                                                                                        <  2:77829/251‑1 (MQ=255)
               cTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGtttt                                                                                                                                                                                                                         >  1:472422/1‑251 (MQ=255)
                   ggtgATTGGCGTCGCTCCTATCGCCCATCGCTGTGTTATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAAcc                                                                                                                                                                                                                     <  2:119422/251‑1 (MQ=255)
                          ggCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGAc                                                                                                                                                                                                               <  1:224883/250‑1 (MQ=255)
                                      ctCGCCCATCGCTGTGATAGTGCCCGGCGCCACGGCGCGGATTTCTGCCTTAACCCCAGCGTCACTGATGTATGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTCTGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGGTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCAt                                                                                                                                                                                                    <  2:387748/248‑1 (MQ=255)
                                       tcgaccATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCAttt                                                                                                                                                                                                  <  1:237983/246‑1 (MQ=255)
                                                gCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGGTTTAACCTCGGACGCGAAGCGCATTTCAATAACGcc                                                                                                                                                                                        >  1:326007/1‑251 (MQ=255)
                                                                           cgGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGAGCCGATGTGATTATCGACACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAAACTCGGACGCGAAGCGCAATTCAATAACGCCAAAATTGTCTTTTCTCGCGCGTACAGc                                                                                                                                                             >  1:292300/1‑251 (MQ=255)
                                                                                ttCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGcgaacc                                                                                                                                                        >  1:406752/1‑251 (MQ=255)
                                                                                 tCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGTCGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGACGcgaaccg                                                                                                                                                       >  1:458956/1‑251 (MQ=255)
                                                                                 gcttccTTAACCCCAACGGCACTGATGTAGGTAAAGAGATCAAAACGCGGACCGGCAAGCAGGGTGCGGAGGGGATTATCGAAACCAGCGGCTACGCCTACGCGCTGCAATCGGCGCTCCGCGGACTGGCTTATGGCGGCACCATCTCCTAGGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCGCGGACGCGAAGCTCATTTCAATAGCGCCAAAAATGTCTTCTCTCGCGCGTGCAGctaaccg                                                                                                                                                       <  2:453082/247‑1 (MQ=255)
                                                                                   tGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAAc                                                                                                                                                                                           <  1:22226/213‑1 (MQ=255)
                                                                                     ccTTAACCCCAGCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGAGGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCgaaccgaacc                                                                                                                                                   <  2:244249/251‑1 (MQ=255)
                                                                                               aTCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGCACCGAACCCGGATTATcc                                                                                                                                         >  1:5431/1‑251 (MQ=255)
                                                                                                           gTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCg                                                                                                                             <  2:113978/251‑1 (MQ=255)
                                                                                                                 aaaGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCCTCTCCTATGTCGCGGTTACCAAGACGTTTGCCGAAGGTTTTAACCTCTGACGCGCAGCGCATTTCCATAACGCCAACATTGTCTTCTCTCGCGCGTGCAGCGCACCGACACCGGATTTTCCCCCCTTGAGCCCCAAGCg                                                                                                                       >  2:411442/1‑251 (MQ=255)
                                                                                                                    gagaTCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGACGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTAt                                                                                                                    >  1:324433/1‑251 (MQ=255)
                                                                                                                     agagcAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCGCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGa                                                                                                                                 <  2:240934/233‑1 (MQ=255)
                                                                                                                     agaTCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCCCCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGCAGCGCATTTCAAAAACGCCAAAATTATCTTCTCTCGCGCGTGCAGCTCACCGAAACCGGATTATCCACGCTGGAGCCGCAAGCGTAtt                                                                                                                   >  1:146564/1‑251 (MQ=255)
                                                                                                                     agaTCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGa                                                                                                                                 >  1:240934/1‑237 (MQ=255)
                                                                                                                        tCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCCTCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGt                                                                                                                                                                  >  2:475800/1‑201 (MQ=255)
                                                                                                                        tCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGt                                                                                                                                                                  <  1:475800/201‑1 (MQ=255)
                                                                                                                          aaaaCGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGGTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATACGGCCTGGAGCCGCAAGCGTATTga                                                                                                                 >  1:399094/1‑248 (MQ=255)
                                                                                                                          aaaaCGCTGACCGGCAAGCAGGGTGCCGATGAGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCCCGTTTAACAAGCCGTTTGCCGAAGGTTTTAACCTCGGGAGCAAAGCGCATTTCAATAACGCCACAATTGTCTTCTCTCGCGCGTGCAGCGCACCGCACCCGGGTTATCCGCGCTGGGGCCGcccgcggcttgg                                                                                                                 >  2:290265/1‑237 (MQ=255)
                                                                                                                              cGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCTTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAAcc                                                                                                          <  1:31436/251‑1 (MQ=255)
                                                                                                                               gCTGACCGGCAAGCAGGTTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCTCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCt                                                                                                         <  1:462668/251‑1 (MQ=255)
                                                                                                                               gCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCt                                                                                                         <  1:166814/251‑1 (MQ=255)
                                                                                                                                      ggCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGACGGTTTTAAACTCGGACGCGAAGCCCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAAGAACCTGTTgggg                                                                                                  >  2:380219/1‑250 (MQ=255)
                                                                                                                                      ggCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGa                                                                                                  <  1:127359/251‑1 (MQ=255)
                                                                                                                                         aaGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAAct                                                                                               >  1:126816/1‑251 (MQ=255)
                                                                                                                                                  gCCGATGGGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCGCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGGCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCGCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGATCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGaa                                                                                      <  2:413929/251‑1 (MQ=255)
                                                                                                                                                   ccGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCGTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAAc                                                                                     <  2:408685/251‑1 (MQ=255)
                                                                                                                                                   ccGAGGTGATTATCGAAACCAGCGGCCACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCGCCATCTCCTATGTCGCGTTTTCCAAGCCGGTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCtgatgcac                                                                                     <  2:9673/251‑4 (MQ=255)
                                                                                                                                                    cGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCCTCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCCAGCGTATTGAAGAAACCTGTTGGGAACTGTTGATGAACg                                                                                    >  2:477768/1‑251 (MQ=255)
                                                                                                                                                      aTGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGt                                                                                  <  1:111798/251‑1 (MQ=255)
                                                                                                                                                        gtgATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTa                                                                                <  1:372780/251‑1 (MQ=255)
                                                                                                                                                         tgATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCCCCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGCACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAAAGGTTAt                                                                               >  1:28307/1‑251 (MQ=255)
                                                                                                                                                                     ccAGCGGCTACGGCGACTCGCTGCAATCGGCGCTCCGCGGTCTGGCTGATGGCGGCACCATCTCCTATGTCGCGTTTTCCAAGCCGTGTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCGTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGTGCCGCAAGCGGATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGac                                                                   <  2:146564/251‑2 (MQ=255)
                                                                                                                                                                      cAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAGCCTCGGACGCGAAGCGCATTTCAATAACGCCaaa                                                                                                                                                                                     <  1:127536/136‑1 (MQ=255)
                                                                                                                                                                       aGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGa                                                                 <  1:477768/251‑1 (MQ=255)
                                                                                                                                                                          ggCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAActg                                                                                              <  1:95687/219‑1 (MQ=255)
                                                                                                                                                                          ggCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGGTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTCTTGAAGAAACCTGTTGGGAActg                                                                                              >  2:95687/1‑219 (MQ=255)
                                                                                                                                                                                 ccGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTTCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAGCGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGCAGATTTAATCgcc                                                       >  1:367213/1‑249 (MQ=255)
                                                                                                                                                                                 ccGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGAc                                                       >  1:176622/1‑251 (MQ=255)
                                                                                                                                                                                     cgcgCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTCTCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGGTGAACGGTTATCTCACTTGCGAAGATTTAATCGACCCgg                                                   >  1:201580/1‑251 (MQ=255)
                                                                                                                                                                                          ggCAATCGGCGCTCCGGGGTTTGGCTTATGGCGGCACCATCGCCTAGTGCGCGTTGGCCAAGCCGTTGGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAGCGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTGGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGgtggtg                                              <  2:367213/250‑1 (MQ=255)
                                                                                                                                                                                           gCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGgtggtg                                              <  2:176622/250‑1 (MQ=255)
                                                                                                                                                                                                           ggTCTGGCTTATGGCGGCACCATCTCCTCTGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGc                                                                      >  1:395193/1‑210 (MQ=255)
                                                                                                                                                                                                           ggTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGTGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTCTGCCAATAGccc                             <  1:395511/251‑1 (MQ=255)
                                                                                                                                                                                                           ggTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGc                                                                      <  2:395193/210‑1 (MQ=255)
                                                                                                                                                                                                            gTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCg                            <  2:28307/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggcggcACCATCTCCTATGTCGCGTTTACCAAGCCGTTGGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCtata                  <  2:5431/250‑1 (MQ=255)
                                                                                                                                                                                                                            cACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGTATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGt            <  1:485664/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGt        <  2:326007/247‑1 (MQ=255)
                                                                                                                                                                                                                                      gaTGTCGTGTTTACAAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAAGTGCCTTCTCTCGCGCGTGCAGAGAAACGAACCCGGATGAGCCGGGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGACCGGTTATCTCAATTGCGAAGATTTAATCGCCCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGTAGTATGTCGTCCa  <  2:376008/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGAcc   <  1:290265/249‑1 (MQ=255)
                                                                                                                                                                                                                                               gTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATcc                                                                                                                                         <  2:170596/108‑1 (MQ=255)
                                                                                                                                                                                                                                               gTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATcc                                                                                                                                         >  1:170596/1‑108 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
AACGCGCTGGCGCTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGACCA  >  NC_000913/1375495‑1375975

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