Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,566,267 T→G 5.5% A206A (GCA→GCC eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,566,2670TG 5.5% 306.3 / 6.1 114A206A (GCA→GCCeutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base T (50/57);  new base G (2/5);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 114 total aligned reads displayed.

GCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGG  >  NC_000913/2566119‑2566415
                                                                                                                                                    |                                                                                                                                                     
gCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTg                                                                                                                                                      <  1:192238/150‑1 (MQ=255)
 cGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGc                                                                                                                                                     <  1:547339/150‑1 (MQ=255)
 cGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGc                                                                                                                                                     <  1:612082/150‑1 (MQ=255)
      cacCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTTAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGa                                                                                                                                                >  2:1028385/1‑150 (MQ=255)
      cacCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGa                                                                                                                                                >  2:1725752/1‑150 (MQ=255)
        ccTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAg                                                                                                                                              >  1:1306868/1‑150 (MQ=255)
           gCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGACGccc                                                                                                                                           <  1:288439/150‑3 (MQ=255)
           gCGGGCAGCAAACCGGTTTGCAGCAATGCCTGCGCGAGCGGTGAGAGCGGTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCAGGCACCGACGccc                                                                                                                                           <  1:46333/150‑3 (MQ=255)
           gCGGGCAGCAAACCGGTCTGCATCAACGCCTGCGCGAGCGGGGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGGCACCTGCACCAGCTGCGCGCCGGGCAGCGAACGGGCGTCGGGGCAGGCACCGACGccc                                                                                                                                           <  1:1517572/150‑3 (MQ=255)
           gCGGGCAGCAAACCGGGTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGGCACCTGCACCAGCTGCGCGCCGGGCAGCGAACGGGCTTCGGTGCCTGCACCGAAGCAc                                                                                                                                           <  1:120730/150‑1 (MQ=255)
               gcagcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGGGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCAt                                                                                                                                       >  1:1725290/1‑150 (MQ=255)
               gatgcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCTCCGAAGCACTCAt                                                                                                                                       >  2:1009435/4‑150 (MQ=255)
                 agcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATcc                                                                                                                                     <  2:1375821/150‑1 (MQ=255)
                  gcaAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCca                                                                                                                                    >  1:685991/1‑150 (MQ=255)
                     aaCCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACaa                                                                                                                                 >  2:756120/1‑150 (MQ=255)
                       ccGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGACCGGGCGTCGGTGCCTGCACCGAAGCACTCCTCCACAatc                                                                                                                               >  2:75883/1‑150 (MQ=255)
                         ggTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATAAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGACGCACGCATCCACaagcct                                                                                                                             >  1:63098/1‑146 (MQ=255)
                         ggTTTGCATCAATGCCTGCGCGAGCGGGGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGGCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAatcat                                                                                                                             <  2:75576/150‑1 (MQ=255)
                            ttGCATCAATGCCTGCGCGAGCGGTGATAGCGTTCCGTCAATCACTTCGACAATGAGTTCTGCCATCCGCCGGGTCACCTGCACGAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTg                                                                                                                          <  2:576780/150‑1 (MQ=255)
                               caacAATTCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCTGTTCTGCCATCCCCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACCGGCGTCGGTGCCTCCTCCGAAGCACTCATCCACAATCATCAGcgc                                                                                                                      >  2:452047/4‑148 (MQ=255)
                                     tGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGttt                                                                                                                <  2:196646/150‑1 (MQ=255)
                                     tGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCGGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGttt                                                                                                                <  2:295036/150‑1 (MQ=255)
                                      gCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTa                                                                                                               <  2:711942/150‑1 (MQ=255)
                                       ccTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTAt                                                                                                              >  1:997625/1‑150 (MQ=255)
                                        tttcgcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCAGCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCTGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATg                                                                                                             <  1:220355/147‑1 (MQ=255)
                                          gcgcGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAg                                                                                                           <  2:1229638/150‑1 (MQ=255)
                                          gcgcGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGGCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAg                                                                                                           <  2:104163/150‑1 (MQ=255)
                                            gcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCg                                                                                                         <  2:1491840/150‑1 (MQ=255)
                                            gcgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCg                                                                                                         <  2:1696536/150‑1 (MQ=255)
                                             cgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGt                                                                                                        >  1:138786/1‑150 (MQ=255)
                                             cgAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGt                                                                                                        >  1:566122/1‑150 (MQ=255)
                                               aGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTaa                                                                                                      >  1:1358499/1‑150 (MQ=255)
                                                 cGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAAc                                                                                                    >  1:428475/1‑150 (MQ=255)
                                                  ggTGAGAGCGTTCCTTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCCCCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAAcc                                                                                                   <  2:917909/150‑1 (MQ=255)
                                                  ggTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAAcc                                                                                                   <  1:975023/150‑1 (MQ=255)
                                                  ggTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAAcc                                                                                                   >  1:189795/1‑150 (MQ=255)
                                                   gTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGACATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCCCAATCATCTGCCCCGGTTTATGCGCGTAAAccc                                                                                                  >  1:1058307/1‑149 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCGCCGGGTTCTGAGCGTAATCCAcg                                                                                                >  2:334486/1‑150 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcg                                                                                                >  1:1160333/1‑150 (MQ=255)
                                                     gagaGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCAcg                                                                                                <  1:640367/150‑1 (MQ=255)
                                                      agagCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTACGCTCATAAAccgggg                                                                                               <  1:98799/150‑5 (MQ=255)
                                                           gTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCg                                                                                          >  2:1001601/1‑150 (MQ=255)
                                                              ccGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                       <  1:1328593/150‑1 (MQ=255)
                                                              ccGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTAAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATACACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTgg                                                                                       >  2:475831/1‑150 (MQ=255)
                                                               cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                      <  2:742396/150‑1 (MQ=255)
                                                               cGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGGCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGc                                                                                      <  1:418520/150‑1 (MQ=255)
                                                                 tCAATCACTTCGACACTCACTTCTGCCATCCGCCGGGTCACCTGCACCAGCTCCGCGCCTGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAATCCACGCGCCCGTGGCTg                                                                                    <  2:82699/150‑1 (MQ=255)
                                                                     tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCgg                                                                                <  1:78910/150‑1 (MQ=255)
                                                                     tCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGGCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCgg                                                                                <  1:1430492/150‑1 (MQ=255)
                                                                         ttCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTc                                                                            <  2:1290642/150‑1 (MQ=255)
                                                                           cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACACTCCTCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCca                                                                          >  2:1576251/1‑150 (MQ=255)
                                                                           cGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCca                                                                          >  2:579937/1‑150 (MQ=255)
                                                                            gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcag                                                                         <  2:1358499/150‑1 (MQ=255)
                                                                            gACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcag                                                                         >  1:400306/1‑150 (MQ=255)
                                                                               aaTCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCcagcag                                                                      <  1:1624593/150‑1 (MQ=255)
                                                                                 gcAGTGCTGCCATCCGCCGGGGCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGc                                                                    <  1:1312722/149‑1 (MQ=255)
                                                                                   aGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCCCCATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCgc                                                                  >  1:112887/1‑149 (MQ=255)
                                                                                   aGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCgg                                                                  >  1:193460/1‑149 (MQ=255)
                                                                                        tGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGACGCACTCATCCACAATCATCTGCCCCGGTTTATGTGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGaccacc                                                             >  1:512716/1‑150 (MQ=255)
                                                                                         gCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCg                                                            >  2:1691679/1‑150 (MQ=255)
                                                                                          ccATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGACCGGGCGTCGGTGCCTGCACCGCAGCACTCATCCACAATCCTCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGCCCACCgc                                                           >  1:1607374/1‑149 (MQ=255)
                                                                                              ccgccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGtt                                                       >  1:1597537/1‑150 (MQ=255)
                                                                                               cgccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGCCCCCCGCCGTTg                                                      >  1:786678/1‑150 (MQ=255)
                                                                                               cgccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTGTATGAGCGTAAACCACGCGCCCGTGGCGGTCGGTTTCCAGCAGGCGACCACCGACGTTg                                                      >  2:95950/1‑150 (MQ=255)
                                                                                                gccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGa                                                     >  1:725555/1‑150 (MQ=255)
                                                                                                 ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGag                                                    <  2:1160333/150‑1 (MQ=255)
                                                                                                 ccgGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGag                                                    <  1:10451/150‑1 (MQ=255)
                                                                                                 ccgGGGCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGGCGGTGCCGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGag                                                    <  1:1371928/150‑1 (MQ=255)
                                                                                                   ggggCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggc                                                  <  1:680488/150‑1 (MQ=255)
                                                                                                   gggTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCCCCGACGGTGcggc                                                  >  1:576764/1‑150 (MQ=255)
                                                                                                     gTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCCTCCACAATCATCTGCCCCGGTGTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGGTGaggcag                                                >  1:411329/1‑150 (MQ=255)
                                                                                                       cACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGCAGCACTCATCCACAATCATCTGCCCCGGTTGATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCCGGCGACCCCCGACGGTGcggcaggc                                              >  1:1103348/1‑150 (MQ=255)
                                                                                                        aCCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggcaggc                                              <  1:1194035/150‑1 (MQ=255)
                                                                                                         ccTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGgcaggcag                                            <  1:1725752/150‑1 (MQ=255)
                                                                                                              cccAGCTGCGCGCCGGTCAGCGAACGGGCGTGGGTGCCTGCACCGAAGCACTCATCCACAATCTTCTGCCCCGGTTTATGAGCGTTAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCg                                       <  2:814791/149‑1 (MQ=255)
                                                                                                                cAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCt                                     <  2:685991/150‑1 (MQ=255)
                                                                                                                cAGCTGCGCGCCGGGCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCt                                     <  2:823744/150‑1 (MQ=255)
                                                                                                                   cTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAt                                  >  2:577570/1‑150 (MQ=255)
                                                                                                                       gcgcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAttttt                              <  2:8613/150‑1 (MQ=255)
                                                                                                                       gcgcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGAttttt                              >  2:1220475/1‑150 (MQ=255)
                                                                                                                         gcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTcc                            <  2:400306/150‑1 (MQ=255)
                                                                                                                         gcCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTcc                            <  2:1597537/150‑1 (MQ=255)
                                                                                                                          ccGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCg                           <  2:1306868/150‑1 (MQ=255)
                                                                                                                          ccGGGCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTAGGAGCGTAAACCACGCGCCCGTGGCTGTCGGGTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCg                           <  1:42764/150‑1 (MQ=255)
                                                                                                                           cGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCgg                          <  1:791085/150‑1 (MQ=255)
                                                                                                                           cGGTCAGCGAACGGGCGTCGGAGCCTGCACCGAAGCACTCATCCACAATCCTCTGCCCCGGTTTAGGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCgg                          <  1:75883/150‑1 (MQ=255)
                                                                                                                                 gCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGaggcaggc                                              <  2:1648010/124‑1 (MQ=255)
                                                                                                                                  cGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGaa                   >  2:965750/1‑150 (MQ=255)
                                                                                                                                   gAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAAc                  <  2:428475/150‑1 (MQ=255)
                                                                                                                                      cGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAgg               >  1:895527/1‑150 (MQ=255)
                                                                                                                                       gtgcttCGGTGCAGGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAggg              >  2:722752/6‑150 (MQ=255)
                                                                                                                                       gggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAggg              <  2:2470/150‑1 (MQ=255)
                                                                                                                                        ggcggcGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGTACAGGGc             <  1:1063591/150‑1 (MQ=255)
                                                                                                                                        ggCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGc             >  1:16224/1‑150 (MQ=255)
                                                                                                                                         gcttCGGTGCAGGCACCGACGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCg            >  1:359485/4‑150 (MQ=255)
                                                                                                                                          cGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGt           >  2:1507281/1‑150 (MQ=255)
                                                                                                                                                   cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCgg  >  2:408289/1‑150 (MQ=255)
                                                                                                                                                   cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCgg  >  2:402036/1‑150 (MQ=255)
                                                                                                                                                   cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCgg  >  1:1574625/1‑150 (MQ=255)
                                                                                                                                                   cTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCGGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCgg  >  1:1114247/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTAACACCTGCGGGCAGCAAACCGGTTTGCATCAATGCCTGCGCGAGCGGTGAGAGCGTTCCGTCAATCACTTCGACAATCAGTTCTGCCATCCGCCGGGTCACCTGCACCAGCTGCGCGCCGGTCAGCGAACGGGCGTCGGTGCCTGCACCGAAGCACTCATCCACAATCATCTGCCCCGGTTTATGAGCGTAAACCACGCGCCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGG  >  NC_000913/2566119‑2566415

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: