Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,795,562 C→T D795N (GAT→AAT)  pheT ← phenylalanine‑‑tRNA ligase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,795,5620CT98.2% 185.1 / NA 55D795N (GAT→AAT) pheTphenylalanine‑‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/19);  minor base G (0/1);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCTGC  >  NC_000913/1795322‑1795802
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
gcgcCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAatgat                                                                                                                                                                                                                                        <  2:140366/251‑1 (MQ=255)
    cGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCt                                                                                                                                                                                                                                    >  1:173386/1‑251 (MQ=255)
    cGTGCTGTAATGGGAATATCCTAGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCt                                                                                                                                                                                                                                    <  1:378344/251‑1 (MQ=255)
          gTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAatgat                                                                                                                                                                                                                                        >  1:463144/1‑241 (MQ=255)
          gTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAatgat                                                                                                                                                                                                                                        <  2:463144/241‑1 (MQ=255)
                                     gggTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCAACTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCACGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGCTGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTgg                                                                                                                                                                                                   >  2:5790/1‑251 (MQ=255)
                                            gTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGt                                                                                                                                                                                            >  1:18118/1‑251 (MQ=255)
                                             tCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGt                                                                                                                                                                                            >  2:152933/1‑250 (MQ=255)
                                                    cGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCAACTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCa                                                                                                                                                                                     >  1:527818/1‑250 (MQ=255)
                                                       tGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAAtct                                                                                                                                                                                 >  1:559159/1‑251 (MQ=255)
                                                           tCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTGGGCGACGGTAGCGGCAATctcct                                                                                                                                                                              <  2:336165/250‑1 (MQ=255)
                                                               aTCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCtcttc                                                                                                                                                                          <  2:173386/250‑1 (MQ=255)
                                                                       gATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAgtgt                                                                                                                                                                 >  1:487837/1‑251 (MQ=255)
                                                                           taaGTTACCAAAACCAGAGAGTGTCACCTGTTCGCCGGTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTGGGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCGGGAATCGCTCTTTTAATGCCTCTTCACATTTGGCGACGTTAGCGGCAATCTCCTCTGCCTCGAGTGTa                                                                                                                                                                <  2:545574/247‑1 (MQ=255)
                                                                           gAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCACTCTCCTCTTCTTCGCGTGTACg                                                                                                                                                              >  1:545574/1‑250 (MQ=255)
                                                                                      aaCCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACGCATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTg                                                                                                                                                  <  1:395043/251‑1 (MQ=255)
                                                                                          agagagTTTCAACTGTTCGCCGTTTTACAGAGCGCGACGGATCTCTTCGCAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATCGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAgg                                                                                                                                              >  2:391348/1‑251 (MQ=255)
                                                                                                      cTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATg                                                                                                                                  >  1:299588/1‑251 (MQ=255)
                                                                                                               gTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGctc                                                                                                                         >  1:345719/1‑251 (MQ=255)
                                                                                                                ttttCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATTCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGctct                                                                                                                        >  2:317676/1‑251 (MQ=255)
                                                                                                                 tttCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCtt                                                                                                                       >  2:214693/1‑251 (MQ=255)
                                                                                                                 tttCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTGCAATGCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCGCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTAGGGCGAGGctct                                                                                                                        <  2:527818/250‑1 (MQ=255)
                                                                                                                     cAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGAGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATAc                                                                                                                   >  1:337891/1‑251 (MQ=255)
                                                                                                                          cgcgACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTc                                                                                                              >  2:508752/1‑251 (MQ=255)
                                                                                                                           gcgACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGTATCAGGCGTATGTCGAGTCTCTTATACCCCTcc                                                                                                             <  2:345719/251‑1 (MQ=255)
                                                                                                                                  aTCTCTTCGAAAAGCAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAACCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCACTCTCCTCTTCTTCGCGTGTACGGCTGGTATCTTGCAGGATCAGCCTTATGGCGAGGCTCTTATCCCCCTCCGcaccac                                                                                                      >  2:171746/1‑251 (MQ=255)
                                                                                                                                         cGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTAcc                                                                                               <  1:322864/251‑1 (MQ=255)
                                                                                                                                         cGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTAcc                                                                                               <  1:214693/251‑1 (MQ=255)
                                                                                                                                          gAAAAACAGTTCAAACAGTTCTTTGGCATCCCGCTTGCTAAGCCAAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACCCCCTTAccc                                                                                              >  2:278351/1‑250 (MQ=255)
                                                                                                                                                     cAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAAGTCCTCAATGATGCCTGGAATCGCTCTGTTAATGCCTCTACACATGGGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTGGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGt                                                                                    <  2:337891/250‑1 (MQ=255)
                                                                                                                                                          aGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCGTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATa                                                                              <  1:251419/251‑1 (MQ=255)
                                                                                                                                                           gTTCTTTGGCATCCCGCTTGCTACGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCACATa                                                                              >  2:524037/1‑250 (MQ=255)
                                                                                                                                                              cTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTACTGCCTCTACACATTTGGCTACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACTGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGtt                                                                          >  1:307430/1‑251 (MQ=255)
                                                                                                                                                                ttGGCATCCCGCTTGCTTAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGTTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCACTTATCCCCCTCCGCAACGCCCTTACCGCAGTACACGTCAAATAAGttt                                                                         <  2:330245/250‑1 (MQ=255)
                                                                                                                                                                  ggCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGTTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATTCCCCTCCGCAACACCCTTACCGCGGTGCACGTCAAAGAAGTTTACg                                                                      <  2:532869/251‑1 (MQ=255)
                                                                                                                                                                  ggCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACg                                                                      <  1:520249/251‑1 (MQ=255)
                                                                                                                                                                     aTCCCGCTTGCTAAGCCCAAGCTGAGCAAACAGATATTCTGACATTTCAGCTTGTGTAAGCGCCATAGGTTCAATTCCTCTATGATGCCTGGAATCGCTCTTTGAATGCCTCTACACATGTGGCGTCGGTAGCGGCAATCTCCTCTTCGTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTTTGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCa                                                                   <  1:233757/251‑1 (MQ=255)
                                                                                                                                                                              gCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCACCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTg                                                           <  1:317676/250‑1 (MQ=255)
                                                                                                                                                                                    cccAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACg                                                    >  1:202088/1‑251 (MQ=255)
                                                                                                                                                                                      cAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGCTGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTCAACGTCAAATAAGTTTACGCCAACTACCTGATTTACGcc                                                  >  2:160638/1‑251 (MQ=255)
                                                                                                                                                                                          cTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACt                                              >  1:440038/1‑251 (MQ=255)
                                                                                                                                                                                            tATCAAACAGATATTCTGACATTTCTGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGAAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACttt                                            >  2:576001/1‑251 (MQ=255)
                                                                                                                                                                                                     gATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTc                                   >  1:311187/1‑251 (MQ=255)
                                                                                                                                                                                                        atTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGGGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTAAGCAAACCACCTGATTTACGCCAACTTTATTACATTCgg                                 >  2:136101/1‑250 (MQ=255)
                                                                                                                                                                                                          tCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGAt                               >  1:499897/1‑250 (MQ=255)
                                                                                                                                                                                                            tGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGCTCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATaa                             >  2:127991/1‑250 (MQ=255)
                                                                                                                                                                                                              aCATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATaaaa                           <  2:487837/250‑1 (MQ=255)
                                                                                                                                                                                                                     aGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGTCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGCGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCt                   >  1:133834/1‑251 (MQ=255)
                                                                                                                                                                                                                     aGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCt                   >  1:479186/1‑251 (MQ=255)
                                                                                                                                                                                                                     aGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAAAACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACATGATTTACGCCAACTTTATTCCATTCGGATAAAATATCCgat                   >  1:100191/1‑249 (MQ=255)
                                                                                                                                                                                                                             tAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTCCGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAAc           >  2:70789/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cgcCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGttt       <  1:199918/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCt    >  1:435727/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cATAGGATCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGACGACGGTAGCGCCAATCTCCTCTTCTTCGAGTGTACAGCTGGTATCTTGCAGGATCAGGATTATGGAGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACCCGTCAAATAAGTTTACGCCAACTACCTGAATTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGACCGTTTTCt    >  1:69651/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGTTAGAGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCCACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCTGc  >  2:33503/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
GCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCTGC  >  NC_000913/1795322‑1795802

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