Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 465,005 G→A V181V (GTC→GTT)  queC ← 7‑cyano‑7‑deazaguanine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913465,0050GA95.1% 125.8 / NA 41V181V (GTC→GTT) queC7‑cyano‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/16);  minor base C (2/0);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAGGGGAATTGCGCGAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGGACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCATATATTGCCGCCAGCGTCAGGAACAAA  >  NC_000913/464759‑465245
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
gAAGGGGAATTGCGCGAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACt                                                                                                                                                                                                                                                >  2:135489/1‑250 (MQ=255)
          tGCGCGAAAAGCTGGAGCAGATGGTAAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGtt                                                                                                                                                                                                                                     <  2:347394/250‑1 (MQ=255)
               gAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGt                                                                                                                                                                                                                               >  1:119996/1‑251 (MQ=255)
               gAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTAGAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGt                                                                                                                                                                                                                               <  2:80924/251‑1 (MQ=255)
                           cAGATGGGTAAGTAAACGTTTTGTGGTGCCGGATGTTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCa                                                                                                                                                                                                                   <  2:119996/251‑1 (MQ=255)
                               tGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCCACGTTTCGTTACGAACTAAATC‑AAGTTTGCCGTAATAATCTGCCAgcgc                                                                                                                                                                                                               >  1:299893/1‑251 (MQ=255)
                                     attaagctgttggttgtgCCGGATACTCAAGCC‑TATACCGGCGACACGCGGAAGAATTTCCTCAACCCGGTTTGCTGCTTCATCGCTGCCATCCCCGTCGGTTTATCGGCCAGATAATGATTCAAACCGGTGGCGCGTAAATTGCATGCCGCACAATGACCGCAACCGTCGTCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAGCTAAATC‑CAGTTTGCCGTAAGACTCTGCCAACGCCCAGGt                                                                                                                                                                                                         <  2:25844/236‑1 (MQ=255)
                                                    gtgCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAAt                                                                                                                                                                                          <  1:483419/251‑1 (MQ=255)
                                                               tctaGCCGCAT‑CCGGCGACACACGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTAGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGAAACCGTCGCCTTTCAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATACACATCAgcg                                                                                                                                                                               >  2:252294/4‑251 (MQ=255)
                                                                cAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTGGCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAgc                                                                                                                                                                                <  1:52164/249‑1 (MQ=255)
                                                                       cAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTc                                                                                                                                                                        >  1:480206/1‑250 (MQ=255)
                                                                        aT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGGTTCa                                                                                                                                                                       >  2:431791/1‑250 (MQ=255)
                                                                         t‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCCGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCaa                                                                                                                                                                      >  1:74752/1‑250 (MQ=255)
                                                                               cGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGa                                                                                                                                                                 >  1:360275/1‑250 (MQ=255)
                                                                                 acacCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAatt                                                                                                                                                              >  2:252270/1‑250 (MQ=255)
                                                                                                aCCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCACCCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGCACTAAATC‑CAGt                                                                                                                                                                                                                                      >  1:198516/1‑164 (MQ=255)
                                                                                                aCCTCAACCCGGGTTTCTGCGTCATCGCTGCCATCTCTGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGTCGCGTAAATTACATGCCGCACAATTTCCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGGTTCGGTACGAACTAAATC‑TAGt                                                                                                                                                                                                                                      <  2:198516/164‑1 (MQ=255)
                                                                                                           gTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGAc                                                                                                                                       <  2:139351/248‑1 (MQ=255)
                                                                                                           gTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGAc                                                                                                                                       >  1:139351/1‑248 (MQ=255)
                                                                                                            ttttCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGc                                                                                                                                    >  1:258353/1‑250 (MQ=255)
                                                                                                                 tGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGtt                                                                                                                              >  2:540352/1‑251 (MQ=255)
                                                                                                                                  cACCGTCGGTTTATCGGCCAGATAATGGTTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACt                                                                                                             >  2:9719/1‑251 (MQ=255)
                                                                                                                                   cccttCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCATGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCGTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTc                                                                                                            <  1:88035/247‑1 (MQ=255)
                                                                                                                                      gTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATc                                                                                                         >  2:180521/1‑251 (MQ=255)
                                                                                                                                           tttatcacccAGATAATGATTCCAACCGTTGGCGCGTAAATTCCATGCCGCCAAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGATACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCAACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGAAGGCATGGTTTAGTGCTTTAACAAACTCATCGCgg                                                                                                     >  2:62175/9‑250 (MQ=255)
                                                                                                                                                cGGCCAGATCATGATTCCAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATcc                                                                                               >  2:489658/1‑251 (MQ=255)
                                                                                                                                                        tAATGATGCAAACCGGTGGCGCGTAAATTACATGCCGCACAATGACCGCAAACGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCACA‑TTTGCCGTAATAATATGCCAACGCCCAGGTTTCCGCTTTATCAATCCCAATCAGCGGCGTTTCAAAAAGAATATCTTTCGCCATGCACAAACTGACGGCATGGGTTAATGCTTTCACAAACTCCTCGCGGCatacagggaagcc                                                                                       >  2:292121/1‑239 (MQ=255)
                                                                                                                                                            gATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCATGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATAGCTATCGCCATGCCCAAACGGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGgag                                                                                   <  2:230372/251‑1 (MQ=255)
                                                                                                                                                                 aaaCCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCAt                                                                                                          >  1:209137/1‑223 (MQ=255)
                                                                                                                                                                 aaaCCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGGTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCAt                                                                                                          <  2:209137/223‑1 (MQ=255)
                                                                                                                                                                   aCCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCCGGTCAACGTTTCGTTACGCACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCACTCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCg                                                                            >  1:536002/1‑251 (MQ=255)
                                                                                                                                                                        tGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCg                                                                       >  2:137475/1‑251 (MQ=255)
                                                                                                                                                                        tGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAAGCCGGGTAGCCGGAGAAATCCGTTTc                                                                        <  2:518382/250‑1 (MQ=255)
                                                                                                                                                                                 aaTTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTCTCGCAGACGCCa                                                              <  1:444633/251‑1 (MQ=255)
                                                                                                                                                                                          ccGCACAAAGACCGCAACCGTCGACTTTACAGCCGTTATAGCAGGTCAACGGTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCTCCCAGGTTTCCGCTATATCAATCAACATCAGCGGCGTTTCAAAACGCATATCTTTAGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACACACTCATCGCGGCAATCCAGGTAACCGGAGAAATCCTTTTCGCAGACGCCAATCATTACg                                                     >  2:373360/1‑251 (MQ=255)
                                                                                                                                                                                             cacaATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCt                                                  >  1:443566/1‑251 (MQ=255)
                                                                                                                                                                                                 atgtccGCAACCGGTGCCTTTAAAGGCGTTATAGCAGGTCAACTGTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAAGCTGCCAGCGCCCAGGTTTCCGCTTTAGGAATCCACTTCAGCGGCGTTAGAAAACGAAGAGCATTCGCCATGTCCAAACTAACGGCATCGTTGAGTGCTTTCACAACCTCATCGCGGCAATCCGTGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTg                                              <  1:550736/247‑1 (MQ=255)
                                                                                                                                                                                                  ggACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGGCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTAGCAATCCACATCAGCGGCGTTTCAAAACGAATAGCTTTCGCCATGCCCAAACTGACGGCAGGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGAGGCCAGTAAGTACGGCTTCTGc                                             <  2:536002/250‑1 (MQ=255)
                                                                                                                                                                                                                      aaaGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCatat                          >  1:287289/1‑250 (MQ=255)
                                                                                                                                                                                                                                 agccggTCAACGCTTCGTTACGAACTAAATC‑CAGTTTGACGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACATGATACGCCAATATTgccgcccac              >  2:67132/5‑248 (MQ=255)
                                                                                                                                                                                                                                              ttCGTTACGAACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCATATATTGCCGCCAGCGTCAGGAACaaa  <  1:540352/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
GAAGGGGAATTGCGCGAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCAT‑CCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGGACTAAATC‑CAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCATATATTGCCGCCAGCGTCAGGAACAAA  >  NC_000913/464759‑465245

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: ‑