Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,195,080 A→G 7.0% Y340Y (TAT→TAC)  thiC ← phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,195,0800AG 7.0% 169.9 / 3.3 58Y340Y (TAT→TAC) thiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base G (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGA  >  NC_000913/4194933‑4195228
                                                                                                                                                   |                                                                                                                                                    
cgaAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGTAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACAt                                                                                                                                                     <  1:18318/149‑1 (MQ=255)
 gaAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGcaca                                                                                                                                                      >  2:340007/1‑147 (MQ=255)
         ctacgccgCACAGATTTCCCAAATTTCGCTGAAGTGTTGAGAGCGGAAATATTCCTGCTGATGGGAGAGGCACCATTACGCCATAATCGAACCGCAGCTAGAGACGATAATGTTCAGGCGTTTCGCGGTCATCTGCACATAGCgcagca                                                                                                                                            <  2:53421/142‑1 (MQ=255)
             gCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTAc                                                                                                                                        <  7:284850/149‑1 (MQ=255)
                gCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTACTGATGATGGGAGAGGCACCATTTCGCCATACTCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGAACATAGCGCAGCAGTACGcc                                                                                                                                     >  5:142288/1‑149 (MQ=255)
                 cacaGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGcaca                                                                                                                                                      <  4:133439/131‑1 (MQ=255)
                 cacaGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGcaca                                                                                                                                                      >  3:133439/1‑131 (MQ=255)
                  acaGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCg                                                                                                                                   >  7:354133/1‑149 (MQ=255)
                  acaGATTTCACAAATTTCCCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGccag                                                                                                                                   >  3:204475/1‑147 (MQ=255)
                   caGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGc                                                                                                                                  <  1:187300/149‑1 (MQ=255)
                      aTTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCCGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATg                                                                                                                               <  3:281562/149‑1 (MQ=255)
                           acaAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCgcag                                                                                                                                              >  6:233300/1‑129 (MQ=255)
                           acaAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCgcag                                                                                                                                              <  5:233300/129‑1 (MQ=255)
                           acaAATTTCGCGGAAGTGGTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATa                                                                                                                           <  4:121097/148‑1 (MQ=255)
                               aTTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTa                                                                                                                                         >  3:53729/1‑130 (MQ=255)
                               aTTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTa                                                                                                                                         <  4:53729/130‑1 (MQ=255)
                               aTTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGaa                                                                                                                      <  3:153010/149‑1 (MQ=255)
                               aTTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGaa                                                                                                                      <  1:340007/149‑1 (MQ=255)
                               aTTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGaa                                                                                                                      <  7:325357/149‑1 (MQ=255)
                                tttCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAg                                                                                                                     <  8:241160/149‑1 (MQ=255)
                                    gcgGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTg                                                                                                                        <  5:292130/142‑1 (MQ=255)
                                    gcgGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTg                                                                                                                        >  6:292130/1‑142 (MQ=255)
                                    gcgGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAAt                                                                                                                 <  6:235087/149‑1 (MQ=255)
                                     cgGAAGTGTTGATAGAGGAAATTTTCCTGATGCGGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCTAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATc                                                                                                                <  6:142288/149‑1 (MQ=255)
                                     cgGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATc                                                                                                                >  8:87901/1‑149 (MQ=255)
                                       gAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCca                                                                                                              >  7:192700/1‑149 (MQ=255)
                                            gtTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCt                                                                                                         >  7:104827/1‑149 (MQ=255)
                                                 tAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTc                                                                                                    <  7:36131/149‑1 (MQ=255)
                                                  agagGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTc                                                                                                    <  8:192700/148‑1 (MQ=255)
                                                    agGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGc                                                                                                 >  5:133831/1‑149 (MQ=255)
                                                    agGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGc                                                                                                 >  7:9981/1‑149 (MQ=255)
                                                     ggAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGcc                                                                                                >  7:204102/1‑149 (MQ=255)
                                                               cTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACGTAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCcag                                                                                        >  2:177941/1‑147 (MQ=255)
                                                               cTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACGTAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCcag                                                                                        >  2:172746/1‑147 (MQ=255)
                                                                 gatgatGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACCTAGCGCAGCAGTAAGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTTTTCCagcagc                                                                                    >  2:345431/1‑149 (MQ=255)
                                                                  atgttgGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCTCACCTTGCTCGGCCTGTTCCagcagc                                                                                    <  8:101098/148‑1 (MQ=255)
                                                                    gatgGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCgtg                                                                                 >  3:60020/1‑149 (MQ=255)
                                                                            agGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGACCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa                                                                         <  4:204475/149‑1 (MQ=255)
                                                                             ggCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAAc                                                                        <  2:17014/149‑1 (MQ=255)
                                                                                aCCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGCACGCt                                                                     >  4:241924/1‑149 (MQ=255)
                                                                                    tttCGCCATAATCGAACCGCCGCTAGAGACGATACCTGTCAGGCGGTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTccc                                                                 <  5:240997/149‑1 (MQ=255)
                                                                                          cATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATTGTGAAGTAATCCACACCTTTCTCGGCCTGTTCCAGCAGAGTGTCGCGGCACGCTTCCCAGGt                                                             >  2:223866/1‑147 (MQ=255)
                                                                                           aTAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGt                                                             >  1:322814/1‑146 (MQ=255)
                                                                                              aTCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATc                                                       >  3:180487/1‑149 (MQ=255)
                                                                                                 gAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTc                                                    <  1:282996/149‑1 (MQ=255)
                                                                                                  aaCCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCg                                                   <  2:102159/149‑1 (MQ=255)
                                                                                                   aCCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTATGATCTTCgg                                                  <  8:354133/149‑1 (MQ=255)
                                                                                                             gagaCGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGt                                         >  6:2950/1‑148 (MQ=255)
                                                                                                                aCGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAAc                                     >  7:227753/1‑149 (MQ=255)
                                                                                                                  gATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCt                                   >  2:989/1‑149 (MQ=255)
                                                                                                                  gATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGAGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCt                                   >  3:4688/1‑149 (MQ=255)
                                                                                                                      ccGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTtctc                               <  1:989/149‑1 (MQ=255)
                                                                                                                       cGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTtctc                               >  8:47052/1‑148 (MQ=255)
                                                                                                                         gTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAg                            >  2:298414/1‑149 (MQ=255)
                                                                                                                             ggCGTTTCGCGGTCATCGGCACGTAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGcc                        <  1:177941/149‑1 (MQ=255)
                                                                                                                             ggCGTTTCGCGGTCATCGGCACGTAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGcc                        <  1:172746/149‑1 (MQ=255)
                                                                                                                                   tCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTgg                     >  2:67805/1‑146 (MQ=255)
                                                                                                                                           aTCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCAc          <  8:104827/149‑1 (MQ=255)
                                                                                                                                            tCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACt         >  3:188432/1‑149 (MQ=255)
                                                                                                                                               gCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTa      >  8:158358/1‑149 (MQ=255)
                                                                                                                                               gCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTa      <  2:322814/149‑1 (MQ=255)
                                                                                                                                                cacaTAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTAc     <  1:298414/149‑1 (MQ=255)
                                                                                                                                                   aTAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGa  >  4:78814/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCTGATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGA  >  NC_000913/4194933‑4195228

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