Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,586 G→C 5.0% intergenic (‑51/+18) guaA ← / ← guaB GMP synthetase (glutamine aminotransferase)/IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,5860GC 5.0% 345.4 / 8.1 99intergenic (‑51/+18)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (41/53);  new base C (2/3);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 102 total aligned reads displayed.

TAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAC  >  NC_000913/2632438‑2632728
                                                                                                                                                    |                                                                                                                                              
tAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATg                                                                                                                                                <  4:35519/149‑1 (MQ=255)
tAAAAACCCAGCTCACGCACTCGTAGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGATTATGAATGTTTTCCGTCATTGACGCTTATTACGAGGCAAGTGAAACAGATAATATAAATCTCCATACATg                                                                                                                                                <  4:40137/149‑1 (MQ=255)
  aaaacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGCACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATg                                                                                                                                                >  4:148102/4‑147 (MQ=255)
    cacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGt                                                                                                                                            <  5:346181/149‑1 (MQ=255)
     acCCAGCTCACGCACGCGGCGCACAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTc                                                                                                                                           <  6:317105/149‑1 (MQ=255)
       ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCgg                                                                                                                                         >  4:309657/1‑149 (MQ=255)
           cTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGa                                                                                                                                     <  6:349666/149‑1 (MQ=255)
           cTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                                                                      >  5:286834/1‑148 (MQ=255)
            tCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                                                                      >  2:201833/1‑147 (MQ=255)
            tCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                                                                      >  1:166950/1‑147 (MQ=255)
              acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                                                                    >  8:295174/1‑147 (MQ=255)
              acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                                                                    <  7:286970/147‑1 (MQ=255)
              acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                                                                    >  7:69092/1‑147 (MQ=255)
              acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGGCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                                                                    <  3:225868/147‑1 (MQ=255)
              acccacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGGCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAaga                                                                                                                                  <  6:239609/146‑1 (MQ=255)
                  acgcGGCGCGCAACCAGTTGCGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAAt                                                                                                                              <  2:258249/149‑1 (MQ=255)
                  acgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAAt                                                                                                                              >  2:286458/1‑149 (MQ=255)
                  acgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAAt                                                                                                                              <  1:80168/149‑1 (MQ=255)
                   cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATc                                                                                                                             <  6:168364/149‑1 (MQ=255)
                       gcgcgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGa                                                                                                                         <  2:322246/149‑1 (MQ=255)
                        cgcgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAg                                                                                                                        <  8:218603/149‑1 (MQ=255)
                            cAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCa                                                                                                                    <  7:90462/149‑1 (MQ=255)
                             aaCCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAaga                                                                                                                                  >  1:91449/1‑134 (MQ=255)
                             aaCCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAaga                                                                                                                                  <  2:91449/134‑1 (MQ=255)
                               ccAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGAc                                                                                                                 >  1:218054/1‑149 (MQ=255)
                                cAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGGAAGTGAATCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                                                                <  8:303469/149‑1 (MQ=255)
                                 aGTTGAGTGTACTGAGACCCGAAGTCGTGAATGAGGATGCTATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgg                                                                                                               <  2:165091/149‑1 (MQ=255)
                                  gTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAag                                                                                                                                   >  1:205428/1‑128 (MQ=255)
                                  gTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAag                                                                                                                                   <  2:205428/128‑1 (MQ=255)
                                   ttGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCgg                                                                                                                                         >  1:350100/1‑121 (MQ=255)
                                   ttGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCgg                                                                                                                                         <  2:350100/121‑1 (MQ=255)
                                   ttGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGc                                                                                                                       <  8:212523/139‑1 (MQ=255)
                                   ttGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGc                                                                                                                       >  7:212523/1‑139 (MQ=255)
                                   ttGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgg                                                                                                               >  7:88601/1‑147 (MQ=255)
                                   ttGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTa                                                                                                             <  6:320405/149‑1 (MQ=255)
                                       gtgtACTGAGAACCGAAGTCAAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTc                                                                                                         <  8:360678/149‑1 (MQ=255)
                                          tACTGAGAACCGAAGTCCATCATGAGGATGCGATGCTTATTAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACATATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       <  6:100782/148‑1 (MQ=255)
                                          tACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       <  1:286458/148‑1 (MQ=255)
                                          tACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       <  1:201833/148‑1 (MQ=255)
                                             tGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                                                                   >  4:50194/1‑149 (MQ=255)
                                                 aaCCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCtt                                                                                               >  5:198622/1‑149 (MQ=255)
                                                       aGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTaa                                                                                         <  8:364073/149‑1 (MQ=255)
                                                       aGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTaa                                                                                         <  2:166950/149‑1 (MQ=255)
                                                        gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAAt                                                                                        <  8:339957/149‑1 (MQ=255)
                                                              aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAc                                                                                  <  2:218054/149‑1 (MQ=255)
                                                                tGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgt                                                                                >  6:101969/1‑149 (MQ=255)
                                                                   ggATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGCATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcgt                                                                             <  3:278832/149‑1 (MQ=255)
                                                                    gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                                                                >  2:302096/1‑113 (MQ=255)
                                                                    gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                                                                <  1:302096/113‑1 (MQ=255)
                                                                    gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGAAATGAAGTCGGGCGAAGAGAATCATGAGAACAGACGGTAGTTCCGGGACTCTTTAGTAATGTTCAcgtcgt                                                                             >  8:183125/1‑148 (MQ=255)
                                                                         gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       >  4:279462/1‑117 (MQ=255)
                                                                         gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       <  3:279462/117‑1 (MQ=255)
                                                                         gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                                                                       >  7:78911/1‑149 (MQ=255)
                                                                         gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                                                                       <  3:309657/149‑1 (MQ=255)
                                                                           tgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTgg                                                                     <  8:69092/149‑1 (MQ=255)
                                                                            gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                       <  6:348142/130‑1 (MQ=255)
                                                                            gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                       >  5:348142/1‑130 (MQ=255)
                                                                              ttATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCTCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCtt                                                                  >  3:90506/1‑149 (MQ=255)
                                                                                 tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCTCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                                                                <  4:90506/148‑1 (MQ=255)
                                                                                   aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACTGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                                                              <  5:101969/148‑1 (MQ=255)
                                                                                   aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActc                                                                                                 <  2:330804/113‑1 (MQ=255)
                                                                                   aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                                                              >  6:320801/1‑148 (MQ=255)
                                                                                     tGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa                                                           <  2:42938/149‑1 (MQ=255)
                                                                                          tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                                                                    >  4:119750/1‑119 (MQ=255)
                                                                                          tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                                                                    <  3:119750/119‑1 (MQ=255)
                                                                                            cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg                                                     >  8:269205/1‑148 (MQ=255)
                                                                                             gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGATCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCa                                                   >  2:142075/1‑149 (MQ=255)
                                                                                             gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                       >  8:348843/1‑113 (MQ=255)
                                                                                             gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                       <  7:348843/113‑1 (MQ=255)
                                                                                             gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCa                                                   <  3:50194/149‑1 (MQ=255)
                                                                                               cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGccc                                                      <  6:215098/144‑1 (MQ=255)
                                                                                               cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGccc                                                      >  5:215098/1‑144 (MQ=255)
                                                                                                aTTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       >  1:322598/1‑94 (MQ=255)
                                                                                                aTTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                                                                       <  2:322598/94‑1 (MQ=255)
                                                                                                   gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt                                              <  8:114230/148‑1 (MQ=255)
                                                                                                   gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt                                              >  1:350008/1‑148 (MQ=255)
                                                                                                    aCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGa                                            <  7:269205/149‑1 (MQ=255)
                                                                                                          aTTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGt                                       <  6:198622/148‑1 (MQ=255)
                                                                                                             ccGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACa                                    >  5:43681/1‑148 (MQ=255)
                                                                                                              cGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACaag                                  <  8:243495/149‑2 (MQ=255)
                                                                                                                    aaGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGc                            >  1:336145/1‑149 (MQ=255)
                                                                                                                      gTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                          >  6:87596/1‑149 (MQ=255)
                                                                                                                         aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                          >  5:245247/1‑146 (MQ=255)
                                                                                                                         aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTa                        <  6:245247/148‑1 (MQ=255)
                                                                                                                           aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTa                     <  6:43681/149‑1 (MQ=255)
                                                                                                                               ataataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACg                   >  2:187179/1‑147 (MQ=255)
                                                                                                                                   tataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtc             >  7:265550/1‑149 (MQ=255)
                                                                                                                                      aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                          >  2:17334/1‑133 (MQ=255)
                                                                                                                                      aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                          <  1:17334/133‑1 (MQ=255)
                                                                                                                                      aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgtc          >  4:156058/1‑149 (MQ=255)
                                                                                                                                        ttCGCCCGACTTCATGTCGGGTGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGa        >  3:43440/2‑149 (MQ=255)
                                                                                                                                           gCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGccc                                                      >  3:25273/1‑100 (MQ=255)
                                                                                                                                           gCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGccc                                                      <  4:25273/100‑1 (MQ=255)
                                                                                                                                           gCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATa      <  1:142075/148‑1 (MQ=255)
                                                                                                                                           gCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATa      <  7:295174/148‑1 (MQ=255)
                                                                                                                                           gCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATa      <  1:187179/148‑1 (MQ=255)
                                                                                                                                           gCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATa      <  8:88601/148‑1 (MQ=255)
                                                                                                                                             ccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa   >  6:146941/1‑149 (MQ=255)
                                                                                                                                              cGACATGAAGTCGGTCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAc  <  7:183125/149‑1 (MQ=255)
                                                                                                                                               gACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAc  >  4:172062/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAC  >  NC_000913/2632438‑2632728

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