Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 364,170 2 bp→AA coding (2135‑2136/3075 nt) lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913364,1700GA94.4% 106.6 / NA 37A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/25);  minor base C (2/0);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 8.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913364,1710GA94.6% 104.7 / NA 37A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/25);  minor base C (1/1);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCG  >  NC_000913/363927‑364407
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                            
gTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTg                                                                                                                                                                                                                                        <  1:346376/251‑1 (MQ=255)
                          aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACt                                                                                                                                                                                                              <  1:268001/251‑1 (MQ=255)
                             gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCTCATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCGGCTCTATGCAAAAATCCATTTCGCTGTTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACGGCGGCCAGGCGCTGATGTGCCCTACTTCTGACCATGCGATCGCGTTCGGTGGCACTa                                                                                                                                                                                                             <  2:193884/249‑1 (MQ=255)
                                   gcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAAGCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACt                                                                                                                                                                                                     <  2:348998/251‑1 (MQ=255)
                                   gcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACAACTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCGCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACt                                                                                                                                                                                                     <  2:121636/251‑1 (MQ=255)
                                         gcgTCAGCAGTTGTTTTTTATCGCCAAGCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGc                                                                                                                                                                                               <  1:166965/251‑1 (MQ=255)
                                                   ttgttTTTTATCGCCAATCCACAGCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGc                                                                                                                                                                                      <  2:343511/250‑1 (MQ=255)
                                                            tgcgtcAATCCACATTTGTGAAAGAAAGCCTGACTGGCGGTTATAGTGCCAACGCTTATTGCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCGGACGCCGCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGctc                                                                                                                                                                             <  2:13216/245‑1 (MQ=255)
                                                                                            aCTGGCGGTTAAATTGCCAACGCTTATTCCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGTTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAAc                                                                                                                                            <  1:212541/251‑1 (MQ=255)
                                                                                               ggCGGTTAAATTGCCAACGCTTATTACCCAGCGCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGt                                                                                                                                         <  2:316030/251‑1 (MQ=255)
                                                                                                      aaaTTGCTAACGCTTATTACCCAGCTCGATGCAAAAATCCCTTTCGCTGGTGGCCAGCTGCGGGCTGGCGTGGGACGCGGCGGGGCGGGGGCCCCTGCGGTTTTCCGCCCGACGCCACTGCTGCCCGGCGCTGGTGGTCACACATTCTGATCCTGCGGGAGCGTTTCGTTGCACTTCCCGTTCTGTGACCAAGAGTTGCCCGGCCCCCCCCGGCCGAGGTGGTTATGGCAGTTCCATCAAccttttacatt                                                                                                                                  >  2:449090/1‑241 (MQ=255)
                                                                                                             cAACGCTTATTACCCAGCTCGATGCAAAAATCCATGTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGc                                                                                                                           <  1:374347/251‑1 (MQ=255)
                                                                                                                    tattaCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATc                                                                                                                     <  1:42474/250‑1 (MQ=255)
                                                                                                                       taCCCAGCTCGATGCAAAAATCCATTTCGCTGGTTGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCTGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTTGAGCGACATCCag                                                                                                                  <  2:181934/250‑1 (MQ=255)
                                                                                                                           cAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTTCGCCCGACGCCCCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGGTCGGGTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGGAGTTCCGGCAGGTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGcc                                                                                                             >  2:32387/1‑250 (MQ=255)
                                                                                                                               tCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGAAGCCACTGCTGCCAGGCGCTGATGTGCCCCACTTCTGACCATGCGGTAGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTc                                                                                                         >  1:417405/1‑251 (MQ=255)
                                                                                                                                         aaTCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGcc                                                                                               <  1:33198/251‑1 (MQ=255)
                                                                                                                                            ccATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATTCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCCCTCCGGCTGCGGTAGGTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGc                                                                                            >  1:305328/1‑251 (MQ=255)
                                                                                                                                            ccATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGc                                                                                            <  1:49624/251‑1 (MQ=255)
                                                                                                                                                         ggtCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGAGTCACACTGAGGTTTTCCGCCAGACGCCCCCGCTGCCAGGAGCTGATGTGCCCCACTTCTGCCCATGCGGTCGCGTTCGGTTGCACTACTCGTACTGTGGGCCAGAGTTGCCCGGCGCTCCCCGGCTGCGGGAGTTCAGGCAGGTCAATCAACTGTTGACCTTGTTGAGCGACATCCAGAGGCGCTTCAACGCTTGCCAGCGGCTTACCATcac                                                                               >  2:47648/1‑249 (MQ=255)
                                                                                                                                                            cAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCcgcg                                                                            >  1:347848/1‑251 (MQ=255)
                                                                                                                                                                             gggACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCGCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTGGCCAGCGTCTTACCATCCAGCGCCACCATCCAGTGCAgg                                                           <  1:44451/251‑1 (MQ=255)
                                                                                                                                                                                      cggGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGtt                                                   >  1:327190/1‑250 (MQ=255)
                                                                                                                                                                                       ggggAGCGTCACACGGAGGTTGTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTAt                                                 <  2:303207/251‑1 (MQ=255)
                                                                                                                                                                                             cGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGGGCAGGAGCTCGTTCTCGCTAt                                           >  2:58995/1‑251 (MQ=255)
                                                                                                                                                                                                 acacTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACg                                       >  1:182928/1‑251 (MQ=255)
                                                                                                                                                                                                 acacTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACGTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTGTACCTTGTGGAGCGACATCCAGAGTCACTTCACCGCGTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGGTATCGCTATGACg                                       <  2:214134/251‑1 (MQ=255)
                                                                                                                                                                                                      gAGGTTTTCCGCCAGACGCCTCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa                                    <  1:146483/249‑1 (MQ=255)
                                                                                                                                                                                                         gTTTTCCGCCAGACGCCACTGCTGCCAGGCGATGAAGTGACTAACTTCTGACCCTGCGATCGCGTTCGGTTGCACTACGCGTACTGTGACCCATAGTTGACCGGCTCTCTCTGGCTGCGGAAGTTCAGGCAGGTCAATCAACAGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCTCTAGCCAGCGGCTTACCATCCAGCGGCACCATCCAGTGCAGGAGCTCGTTATCGCTTTGACgttacaga                                >  1:286090/1‑243 (MQ=255)
                                                                                                                                                                                                            ttCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTGCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTAt                             <  2:379991/250‑1 (MQ=255)
                                                                                                                                                                                                             tCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCGCCGGCTGCGGTAGTGCAGGCAGTTCAATCAACTGTTTACCGTGTGGAGCGACATCCAGAGGCACTGCACCGCTTGCCAGCGGCTTGCCTTTCATCGGCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTAtt                            <  1:59685/250‑1 (MQ=255)
                                                                                                                                                                                                                 ccAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTg                       <  1:53411/251‑1 (MQ=255)
                                                                                                                                                                                                                     aCGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCa                   <  1:71378/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccACTGCTTCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTGGCCCGGCGCTCTCCTGCCGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACtt                <  2:69515/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcCAGTCGCTGATGTGCCCACCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACGGTGAGCCAGAGTGGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCATCTGTTTACCTTGTGGAGCGACATCCATAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGATCTCGTTATCGCTATGATGGAACAGGTAGTCGCTGGTCACTACGATgg          <  2:389490/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCGCTGCTGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGcc    >  1:171750/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cgcTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTGCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTTGTTATCGCTATGACGGAGCAGGTATTCGCTGGTCACTTCGATGGTTTGCCCg  <  1:433077/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                            
GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCG  >  NC_000913/363927‑364407

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