Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,433 C→A 6.3% intergenic (‑256/+141) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,4330CA 6.3% 225.4 / 2.3 74intergenic (‑256/+141)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base A (3/2);  total (40/34)
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.63e-01

CGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCA  >  NC_000913/699286‑699576
                                                                                                                                                   |                                                                                                                                               
cGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATgcg                                                                                                                                                <  8:201774/149‑1 (MQ=255)
 gTAATGGTGATTGTCGATTGGTTAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATgcgg                                                                                                                                               <  6:206094/149‑1 (MQ=255)
          aTTGTCGATTGGGGAAAAAAGGGCTGGTGCGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTGTCTTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                      <  6:309905/149‑1 (MQ=39)
           ttGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATGGAGAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                     >  4:156590/1‑149 (MQ=255)
             gTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                    >  3:103239/1‑148 (MQ=255)
              tCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTCTATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                   <  8:47868/148‑1 (MQ=255)
                gATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                >  3:248119/1‑149 (MQ=255)
                   tcgTGAAAAAAGGTCTGCTCAGACGATATATTTATTGATGACTCGATAATTTTTAGCGGGTTTTATTGAATGTTTATAGTTTACGGGGGCCAAATTGCTGACAAAGTGAGATTTGTTCATGCCGGATGCGGCGTGAACGCCTAATCCgg                                                                                                                             >  2:109823/3‑149 (MQ=37)
                     gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCTAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                           >  5:359767/1‑149 (MQ=255)
                        aaaaaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                         >  4:288923/1‑148 (MQ=255)
                        aaaaaaGGTCTGGTGTGAGGATATATTTATCGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                         >  2:152639/1‑148 (MQ=255)
                         aaaaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                        >  7:324679/1‑148 (MQ=255)
                           aaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                                                                                      >  2:167836/1‑148 (MQ=255)
                            aaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                    <  5:269308/149‑1 (MQ=255)
                                tCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCAt                                                                                                                >  6:82926/1‑149 (MQ=255)
                                    gtgtGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGaa                                                                                                             >  1:296449/1‑148 (MQ=255)
                                     tgtgAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGaaaa                                                                                                           >  7:306598/1‑149 (MQ=255)
                                          ggATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGCGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCaa                                                                                                      <  1:135409/149‑1 (MQ=255)
                                          ggATATATTTATTGAGGAATCGATAATTTTTAGCGGGTGTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCaa                                                                                                      <  8:374421/149‑1 (MQ=255)
                                           gATATATTTATCGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAat                                                                                                     <  1:152639/149‑1 (MQ=255)
                                            atatatTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCATTATCCGGCCTACAAAAGCATGAAAATTCAata                                                                                                    <  6:168205/149‑1 (MQ=255)
                                             tatatTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAata                                                                                                    >  2:70995/1‑148 (MQ=255)
                                             tatatTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTATATATTTTACGGGGGCAAAATTTCTGACAAAGTGCTATTTGTTCATTCCGGATGAGACGTGAACGCCTAATCCGGACTACAAAAGCATGAAAATTCAata                                                                                                    >  5:162071/1‑148 (MQ=255)
                                              atatTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAatata                                                                                                  >  8:82352/1‑149 (MQ=255)
                                                 tttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTg                                                                                               <  4:333349/149‑1 (MQ=255)
                                                        atgaATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGc                                                                                        >  5:253822/1‑149 (MQ=255)
                                                             tCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGt                                                                                   <  4:103239/149‑1 (MQ=255)
                                                              cGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCg                                                                                    >  7:176118/1‑147 (MQ=255)
                                                                   aTTTTTAGAGGGTTTTATTGAATTTTTATATTTTACGGGGGCCAAATTGCTGACAAATTGCGATTTGTTCATGCCGGATTCTGCGTGAACGCCTTATCCGGCCGACCAAAGCATGAAAATTCAATATATTGCAGTAGCTGCGTAGGCCt                                                                             <  4:32248/149‑1 (MQ=25)
                                                                      tttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAATATTCAATATATTGCAGGAGCTGCGTAGGCCTGa                                                                           <  8:176118/148‑1 (MQ=255)
                                                                      tttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGa                                                                           <  5:198274/148‑1 (MQ=255)
                                                                         aGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                        <  2:296449/148‑1 (MQ=255)
                                                                          gCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCg                                                                                    >  3:69188/1‑135 (MQ=255)
                                                                          gCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCg                                                                                    <  4:69188/135‑1 (MQ=255)
                                                                          gCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGc                                                                      <  6:253822/149‑1 (MQ=255)
                                                                          gCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTAATAAg                                                                       <  1:70995/148‑1 (MQ=255)
                                                                            gggTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCg                                                                     <  1:167836/148‑1 (MQ=255)
                                                                             ggTTTTATTGAATGTTTATATTTTACGTGGGCCAAATTGCTTACAAAGTGCGATTTGTTCATGGAGGATGAGGCGTGAACGCCTTATCCGGCCTAGAAAAGCATGAAAATTCAATATATTGCAGGAGATGCGTAGGCCTGATAAGCGTa                                                                   <  6:162071/149‑1 (MQ=25)
                                                                                     tGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCAt                                                             >  7:39425/1‑147 (MQ=255)
                                                                                     tGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTACAGGAGCTGCGTa                                                                                  <  4:44329/126‑1 (MQ=255)
                                                                                     tGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTACAGGAGCTGCGTa                                                                                  >  3:44329/1‑126 (MQ=255)
                                                                                       aaTGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAg                                                          <  8:39425/148‑1 (MQ=255)
                                                                                         tGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  2:38535/1‑147 (MQ=255)
                                                                                         tGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       <  8:306598/149‑1 (MQ=255)
                                                                                             tataTTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                                                   >  1:12706/1‑149 (MQ=255)
                                                                                             tataTTTTACGGGGGACAAATTGCTGACAAAGTGCGATTTGTTCCTGCCGGATGCGGAGTGAACGCCTTATACGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTTCGTAGTCGTGATAAGCGTAGCGCATCAGGCAGttt                                                   >  1:18874/1‑149 (MQ=21)
                                                                                                 ttttACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTa                                                                   >  5:314617/1‑129 (MQ=255)
                                                                                                 ttttACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTa                                                                   <  6:314617/129‑1 (MQ=255)
                                                                                                 ttttACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCg                                               >  8:152465/1‑149 (MQ=255)
                                                                                                  tttACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGt                                              >  3:305955/1‑149 (MQ=255)
                                                                                                   ttACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGt                                              <  7:11163/148‑1 (MQ=255)
                                                                                                          ggCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGTCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcatc                                      <  4:248119/149‑1 (MQ=255)
                                                                                                           gCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcat                                       >  4:272401/1‑147 (MQ=255)
                                                                                                             cAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcatca                                     >  4:319582/1‑147 (MQ=255)
                                                                                                                   gCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCaa                              <  3:288923/148‑1 (MQ=255)
                                                                                                                     tGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGcc                                                                              >  8:371768/1‑98 (MQ=255)
                                                                                                                      gACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGATAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCAc                          >  6:313082/1‑149 (MQ=255)
                                                                                                                        cAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAAcc                            >  2:229374/1‑145 (MQ=255)
                                                                                                                        cAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAAcc                            <  1:229374/145‑1 (MQ=255)
                                                                                                                        cAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACAACGt                        >  2:101025/1‑149 (MQ=255)
                                                                                                                          aaGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTcc                      <  8:324679/149‑1 (MQ=255)
                                                                                                                          aaGTGCGATTTGTACATGCCGGATGCGTCGTGAACGCCTTTTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCGGCGTAGGCGTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGTGCCAACCACGTcc                      <  7:360513/149‑1 (MQ=18)
                                                                                                                                  tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGTAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGt              >  4:32203/1‑149 (MQ=255)
                                                                                                                                      ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAgcg                                                                <  7:104918/95‑1 (MQ=255)
                                                                                                                                      ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAgcg                                                                >  8:104918/1‑95 (MQ=255)
                                                                                                                                      ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGtt          <  4:305955/149‑1 (MQ=255)
                                                                                                                                         aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGc        <  2:115832/148‑1 (MQ=255)
                                                                                                                                         aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGc        >  1:61569/1‑148 (MQ=255)
                                                                                                                                          tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTa      <  7:82352/149‑1 (MQ=255)
                                                                                                                                          tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTa      >  6:352485/1‑149 (MQ=255)
                                                                                                                                          tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTa      >  6:347211/1‑149 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTa      >  7:127373/2‑149 (MQ=25)
                                                                                                                                           gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTAt     >  3:261818/1‑149 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTg         <  5:109374/145‑1 (MQ=25)
                                                                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTg         >  6:109374/1‑145 (MQ=25)
                                                                                                                                              ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa  <  3:156590/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCA  >  NC_000913/699286‑699576

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