Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,142,505 G→C 5.3% R148R (CGC→CGG udk ← uridine/cytidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,142,5050GC 5.3% 422.3 / 12.4 131R148R (CGC→CGGudkuridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (61/63);  new base C (4/3);  total (65/66)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 138 total aligned reads displayed.

CGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTAT  >  NC_000913/2142356‑2142650
                                                                                                                                                     |                                                                                                                                                 
cGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcgg                                                                                                                                                   >  2:1368038/1‑150 (MQ=255)
 ggTTTTTCCGGCCGCGCGGCACGATAATGTCCTCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                                                                  <  2:730380/150‑1 (MQ=255)
 ggTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                                                                  >  2:393192/1‑150 (MQ=255)
 ggTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                                                                  >  2:10737/1‑150 (MQ=255)
 ggTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                                                                  <  1:1515190/150‑1 (MQ=255)
 ggTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                                                                  >  1:546893/1‑150 (MQ=255)
    ttttCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCa                                                                                                                                               <  1:700984/150‑1 (MQ=255)
         ccgccgCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAgg                                                                                                                                          <  1:1129905/150‑1 (MQ=255)
         ccgccgCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAgg                                                                                                                                          >  1:1587866/1‑150 (MQ=255)
         ccgccgCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAgg                                                                                                                                          >  2:135511/1‑150 (MQ=255)
          cgccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGc                                                                                                                                         <  2:453284/150‑1 (MQ=255)
           gccgcGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCa                                                                                                                                        <  2:18421/150‑1 (MQ=255)
             cgcgcgGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGa                                                                                                                                      <  1:68211/150‑1 (MQ=255)
             cgcgcgGCACGATAATGTCAGCATATTGTTTAGAAGGCTCAATGAATTGCATGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGa                                                                                                                                      <  2:671926/150‑1 (MQ=255)
              gcgcgGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGTAACATCGGGCTCACGGGTTTTTGATATTGCGCCAGCACTGAATCCATTGAACGCCCACGCGCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGat                                                                                                                                     <  1:1234052/150‑1 (MQ=255)
                 cgGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATc                                                                                                                                  >  2:551463/1‑150 (MQ=255)
                      cGATAATGCCCGCATATTGTTTAGAAGGCTCAATGAATAGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCg                                                                                                                             <  2:649010/150‑1 (MQ=255)
                       gATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                            <  2:1527838/150‑1 (MQ=255)
                       gATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                            <  1:979226/150‑1 (MQ=255)
                       gATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgg                                                                                                                            <  1:972896/150‑1 (MQ=255)
                        aTAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCggg                                                                                                                           <  1:430090/150‑1 (MQ=255)
                            tGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGATTTTTGATATTGCGACATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTAt                                                                                                                       >  1:1214961/1‑150 (MQ=255)
                                  cATATTGTTTAGCAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCCCGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGa                                                                                                                 >  1:314210/1‑150 (MQ=255)
                                  cATATTGTTTAGAAGGCTCAATGTATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGa                                                                                                                 <  1:1008710/150‑1 (MQ=255)
                                   atatTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCAATGAACGCCCACGCTCCTTAACTTCACGCTTGATGCGGAGCATCAGGCAGATATACAGCTGGGTATCAACGaa                                                                                                                >  2:799730/1‑150 (MQ=255)
                                      ttGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAAt                                                                                                             >  1:1194646/1‑150 (MQ=255)
                                        ggttATAAGGCGCAAAGAATTGCAGGAAAATAGGGCGCACCGGTTTTTGATATTTCGCCATCACAGAATACATTGGACGCCCCCGCGAGTTAACGTCACGCTTGATGCGGCGCCTCAGGCAGATATGCCGCGGGGGTTCAACGAAAATgg                                                                                                           >  2:1808357/3‑150 (MQ=255)
                                        gTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATgg                                                                                                           <  2:655374/150‑1 (MQ=255)
                                         tttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                                                                          >  1:1679821/1‑150 (MQ=255)
                                         tttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                                                                          >  1:176368/1‑150 (MQ=255)
                                         tttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                                                                          >  2:508816/1‑150 (MQ=255)
                                         tttAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCAAGCGGGGTATCAACGAAAATGga                                                                                                          >  2:1724390/1‑150 (MQ=255)
                                           tAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgaga                                                                                                        >  2:1270586/1‑150 (MQ=255)
                                           tAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgaga                                                                                                        >  2:396173/1‑150 (MQ=255)
                                             taaGGCACAATGAATTGCATGAACAGCGGGCGCACGGTTTTTTGATCTTGCTCCATCACTGAATCCATTGAACGCACACGCTCGTTAACGTCACGCTTGAGGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaag                                                                                                      <  2:72581/149‑1 (MQ=255)
                                             gAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaag                                                                                                      <  2:3697/150‑1 (MQ=255)
                                             gAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaag                                                                                                      >  1:457549/1‑150 (MQ=255)
                                                ggCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATTTTGCGCCATCACTGAATCCAATGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagtta                                                                                                   <  1:325844/150‑1 (MQ=255)
                                                 gCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaa                                                                                                  >  1:1432500/1‑150 (MQ=255)
                                                  cTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaag                                                                                                 <  2:1675528/150‑1 (MQ=255)
                                                    cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaagtt                                                                                               <  1:34564/150‑1 (MQ=255)
                                                     aaTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTc                                                                                              <  1:463650/150‑1 (MQ=255)
                                                      aTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCg                                                                                             <  1:551463/150‑1 (MQ=255)
                                                        taaTTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTc                                                                                           <  2:20972/149‑1 (MQ=255)
                                                         aaTTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCa                                                                                          <  1:1793278/150‑1 (MQ=255)
                                                            tGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGc                                                                                       >  2:250260/1‑150 (MQ=255)
                                                             gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCa                                                                                      >  2:1228289/1‑150 (MQ=255)
                                                             gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCa                                                                                      <  1:250655/150‑1 (MQ=255)
                                                              cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCaa                                                                                     >  2:897988/1‑150 (MQ=255)
                                                              cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCAGCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCaa                                                                                     >  1:113127/1‑150 (MQ=255)
                                                                ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                   >  2:1612878/1‑150 (MQ=255)
                                                                ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                   >  2:1529691/1‑150 (MQ=255)
                                                                ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                   >  2:796719/1‑150 (MQ=255)
                                                                ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                   >  1:693637/1‑150 (MQ=255)
                                                                ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                                                                   >  1:1301254/1‑150 (MQ=255)
                                                                    cATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGGATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACgcgc                                                                               >  2:144176/1‑150 (MQ=255)
                                                                         ggCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCg                                                                          >  1:1499826/1‑150 (MQ=255)
                                                                                ggTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGcaa                                                                   >  1:1797245/1‑150 (MQ=255)
                                                                                    ttttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACaaa                                                               >  1:933489/1‑150 (MQ=255)
                                                                                       tGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATg                                                            <  1:2180/150‑1 (MQ=255)
                                                                                           atTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCtt                                                        <  2:497898/150‑1 (MQ=255)
                                                                                             tGCGCCATCACTGAATCCCTTGACCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCg                                                      <  2:294652/150‑1 (MQ=255)
                                                                                                 ccATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGaa                                                  <  1:1228289/150‑1 (MQ=255)
                                                                                                  cATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAat                                                 <  1:1314950/150‑1 (MQ=255)
                                                                                                  cATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAat                                                 <  2:116381/150‑1 (MQ=255)
                                                                                                   aTCACTGAATCCCAGGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCAGCCGTCCGCAACAAAATGCCTTCGAGAatg                                                <  2:879062/150‑1 (MQ=255)
                                                                                                   aTCACTGAATCCATTGAACTGACAGGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTATGTTCGTCACGCAAACGCGCTTCCGTCAGCAACAAAATGCCTTCGAGAttg                                                <  1:88019/150‑1 (MQ=255)
                                                                                                      aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatgatg                                             >  2:599954/1‑150 (MQ=255)
                                                                                                      aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatgatg                                             <  1:1598801/150‑1 (MQ=255)
                                                                                                      aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCAAGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatgatg                                             <  1:1591831/150‑1 (MQ=255)
                                                                                                         taaTCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGAcc                                          <  1:135511/149‑1 (MQ=255)
                                                                                                          aaTCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACct                                         >  1:869518/1‑150 (MQ=255)
                                                                                                           aTCCATTGAACGCCCACGCTCGTTAACGTCACGTTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACctt                                        >  2:1190154/1‑150 (MQ=255)
                                                                                                            tCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCCCGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACcttc                                       >  1:1631139/1‑150 (MQ=255)
                                                                                                              cATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACcttctt                                     >  2:1131498/1‑150 (MQ=255)
                                                                                                                 tGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCgg                                  <  1:508816/150‑1 (MQ=255)
                                                                                                                 tGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCAAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCgg                                  <  1:1724390/150‑1 (MQ=255)
                                                                                                                  gAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                                 <  2:279705/150‑1 (MQ=255)
                                                                                                                  caaCGCCCACGCTCGTTAACGTCACGCTAGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGc                                 <  2:1610464/149‑1 (MQ=255)
                                                                                                                   aaCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCATATATCCAGCGGGGTATCAACGAAACTGGATAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCt                                <  2:1370305/150‑1 (MQ=255)
                                                                                                                    aCGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTc                               <  2:1169265/150‑1 (MQ=255)
                                                                                                                     cGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCACCAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCa                              >  2:38414/1‑150 (MQ=255)
                                                                                                                     cGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCa                              >  2:1482100/1‑150 (MQ=255)
                                                                                                                      gCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCaa                             >  1:306499/1‑150 (MQ=255)
                                                                                                                           cGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTc                        >  2:239865/1‑150 (MQ=255)
                                                                                                                                   aaCGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttc                <  1:42497/150‑1 (MQ=255)
                                                                                                                                     cGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctt              <  1:1251782/150‑1 (MQ=255)
                                                                                                                                       tCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGtttctttc            <  2:1679821/150‑1 (MQ=255)
                                                                                                                                        cACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa           >  2:573902/1‑150 (MQ=255)
                                                                                                                                        cACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa           <  1:1368038/150‑1 (MQ=255)
                                                                                                                                        cACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTACGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCa           <  2:1432500/150‑1 (MQ=255)
                                                                                                                                          cGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATa         <  1:1626751/150‑1 (MQ=255)
                                                                                                                                          cGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATa         >  1:802613/1‑150 (MQ=255)
                                                                                                                                           gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAc        >  2:1015548/4‑150 (MQ=255)
                                                                                                                                           gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAc        >  2:378038/4‑150 (MQ=255)
                                                                                                                                              tGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAcgcg     <  2:707745/150‑1 (MQ=255)
                                                                                                                                              tGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATAcgcg     >  1:1399775/1‑150 (MQ=255)
                                                                                                                                                aTGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTa   <  1:796719/150‑1 (MQ=255)
                                                                                                                                                 tGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTAt  >  2:1640300/1‑150 (MQ=255)
                                                                                                                                                 tGCGGCGCATCAGGAAGATATCCAGCGGGGTATCAACGAAAATGGAGCAGTTAAGTTCGTCACGCAAACGCGCAGCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGGCACCGTTTCTTTCATACGCGTAt  >  1:1644060/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGGTTTTTCCCGCCGCGCGGCACGATAATGTCCGCATATTGTTTAGAAGGCTCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATGATGACCTTCTTCGGCTCAACCGTCACCGTTTCTTTCATACGCGTAT  >  NC_000913/2142356‑2142650

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