Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,665 A→T 5.9% intergenic (+54/+28) yegS → / ← gatR phosphatidylglycerol kinase, metal‑dependent/pseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,169,6650AT 5.9% 361.3 / 8.3 100intergenic (+54/+28)yegS/gatRphosphatidylglycerol kinase, metal‑dependent/pseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted
Reads supporting (aligned to +/- strand):  ref base A (43/51);  new base T (4/2);  total (47/53)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 111 total aligned reads displayed.

TGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGA  >  NC_000913/2169518‑2169813
                                                                                                                                                   |                                                                                                                                                    
tGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     <  2:160620/149‑1 (MQ=255)
tGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     <  2:1342048/149‑1 (MQ=255)
tGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     <  3:1057939/149‑1 (MQ=255)
   cGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAacac                                                                                                                                                  <  3:852053/149‑1 (MQ=255)
    gAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAacac                                                                                                                                                  >  2:1705656/1‑148 (MQ=255)
    gAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAacac                                                                                                                                                  <  1:1705656/148‑1 (MQ=255)
     aaCCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     <  1:951518/144‑1 (MQ=255)
     aaCCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     >  2:951518/1‑144 (MQ=255)
     aaCCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCt                                                                                                                                                >  4:1693010/1‑149 (MQ=255)
      aCCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCtt                                                                                                                                               >  2:1471494/1‑149 (MQ=255)
       ccGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTa                                                                                                                                              <  3:1750428/149‑1 (MQ=255)
            gAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCgg                                                                                                                                         <  2:211274/149‑1 (MQ=255)
                gggCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATa                                                                                                                                     <  3:1693010/149‑1 (MQ=255)
                   cAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAg                                                                                                                                  >  4:1567487/1‑149 (MQ=255)
                       aTTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACa                                                                                                                              <  3:77820/149‑1 (MQ=255)
                        ttttCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACaa                                                                                                                             <  2:787588/149‑1 (MQ=255)
                          ttCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAAta                                                                                                                           <  1:272890/149‑1 (MQ=255)
                             atatTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAAtacta                                                                                                                        >  4:351631/1‑149 (MQ=255)
                              tatTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATAAGGTGTTTACACCGCATCCGGCAt                                                                                                                                      <  1:288198/134‑1 (MQ=39)
                              tatTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATAAGGTGTTTACACCGCATCCGGCAt                                                                                                                                      >  2:288198/1‑134 (MQ=39)
                               atTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTa                                                                                                                                              >  2:1190011/1‑125 (MQ=255)
                               atTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTa                                                                                                                                              <  1:1190011/125‑1 (MQ=255)
                               atTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTaca                                                                                                                      >  3:1619361/1‑149 (MQ=255)
                                ttGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                <  4:157151/138‑1 (MQ=255)
                                ttGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                <  2:1332981/138‑1 (MQ=255)
                                ttGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                >  3:157151/1‑138 (MQ=255)
                                ttGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGaa                                                                                                                                >  1:1332981/1‑138 (MQ=255)
                                ttGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTacac                                                                                                                     >  3:629058/1‑149 (MQ=255)
                                    aaTACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCa                                                                                                                                       <  3:881485/127‑1 (MQ=255)
                                     aTACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaa                                                                                                                <  4:1619361/149‑1 (MQ=255)
                                       aCTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaatt                                                                                                              <  3:1567487/149‑1 (MQ=255)
                                          tCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaat                                                                                                           <  4:629058/149‑1 (MQ=255)
                                           ccGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTaattaatt                                                                                                          <  3:231119/149‑1 (MQ=255)
                                                       cgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAc                                                                                                      <  3:742589/141‑1 (MQ=255)
                                                       cgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAc                                                                                                      >  4:742589/1‑141 (MQ=255)
                                                        gttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            <  1:643843/134‑1 (MQ=255)
                                                        gttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            >  2:643843/1‑134 (MQ=255)
                                                         ttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGt                                                                                                                  >  4:688601/1‑127 (MQ=255)
                                                                 tACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGc                                                                                    <  1:410752/149‑1 (MQ=255)
                                                                    caccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAc                                                                                                      >  4:756331/1‑128 (MQ=255)
                                                                    caccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAc                                                                                                      <  3:756331/128‑1 (MQ=255)
                                                                    caccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                 >  3:1880315/1‑149 (MQ=255)
                                                                    caccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                 >  2:1679350/1‑149 (MQ=255)
                                                                    caccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                 >  1:611204/1‑149 (MQ=255)
                                                                     accaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGt                                                                                <  2:1561383/149‑1 (MQ=255)
                                                                      ccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAAtac                                                                                                                          <  4:1851605/106‑1 (MQ=255)
                                                                       caGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAacac                                                                                                                                                  <  2:377137/81‑1 (MQ=255)
                                                                       caGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAacac                                                                                                                                                  >  1:377137/1‑81 (MQ=255)
                                                                          aTTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                           <  1:1471494/149‑1 (MQ=255)
                                                                           ttGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                          <  2:611204/149‑1 (MQ=255)
                                                                              tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            <  4:1611731/112‑1 (MQ=255)
                                                                              tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            >  3:1611731/1‑112 (MQ=255)
                                                                              tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                          <  4:518507/146‑1 (MQ=255)
                                                                              tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                          >  3:518507/1‑146 (MQ=255)
                                                                               ccACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACct                                                                      <  3:351631/149‑1 (MQ=255)
                                                                                cACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAAta                                                                                                                           >  1:939906/1‑95 (MQ=255)
                                                                                cACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAAta                                                                                                                           <  2:939906/95‑1 (MQ=255)
                                                                                    gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            <  2:1641039/106‑1 (MQ=255)
                                                                                    gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            >  1:1641039/1‑106 (MQ=255)
                                                                                    gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTg                                                                                          <  4:1203973/124‑1 (MQ=255)
                                                                                    gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTg                                                                                          >  3:1203973/1‑124 (MQ=255)
                                                                                     ttGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCa                                                                                  >  4:1401632/1‑131 (MQ=255)
                                                                                     ttGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCa                                                                                  <  3:1401632/131‑1 (MQ=255)
                                                                                     ttGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGaca                                                                >  1:428250/1‑149 (MQ=255)
                                                                                        cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAg                                                             <  2:428250/149‑1 (MQ=255)
                                                                                         gTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     <  3:1220362/60‑1 (MQ=255)
                                                                                         gTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                                     >  4:1220362/1‑60 (MQ=255)
                                                                                         gTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTaattaatt                                                                                                          <  2:929892/103‑1 (MQ=255)
                                                                                         gTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTaattaatt                                                                                                          >  1:929892/1‑103 (MQ=255)
                                                                                               cAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCatcatc                                                      <  4:1880315/149‑1 (MQ=255)
                                                                                               cAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCatcatc                                                      <  1:1679350/149‑1 (MQ=255)
                                                                                                aaaCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCaa                                                                                       >  4:1241367/1‑115 (MQ=255)
                                                                                                aaaCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCaa                                                                                       <  3:1241367/115‑1 (MQ=255)
                                                                                                aaaCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGc                                                            <  4:537498/142‑1 (MQ=255)
                                                                                                aaaCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGc                                                            >  3:537498/1‑142 (MQ=255)
                                                                                                 aaCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCgg                                                    >  3:573569/1‑149 (MQ=255)
                                                                                                   cTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                  <  3:504780/133‑1 (MQ=255)
                                                                                                        ctctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGcc                                                                                   >  2:72348/1‑111 (MQ=255)
                                                                                                        ctctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGcc                                                                                   <  1:72348/111‑1 (MQ=255)
                                                                                                         tctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                          <  1:1391696/119‑1 (MQ=255)
                                                                                                         tctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                          >  2:1391696/1‑119 (MQ=255)
                                                                                                                  ccAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCaa                                                                                       <  3:1469635/97‑1 (MQ=255)
                                                                                                                    aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            >  1:956777/1‑74 (MQ=255)
                                                                                                                    aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                            <  2:956777/74‑1 (MQ=255)
                                                                                                                    aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                           <  4:1823689/107‑1 (MQ=255)
                                                                                                                    aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                           >  3:1823689/1‑107 (MQ=255)
                                                                                                                        agaATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaat                                                                                                           <  2:119617/71‑1 (MQ=255)
                                                                                                                        agaATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaat                                                                                                           >  1:119617/1‑71 (MQ=255)
                                                                                                                           aTCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGAtt                          <  4:573569/149‑1 (MQ=255)
                                                                                                                            tCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGAATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaataataa                                           <  1:1261839/131‑1 (MQ=255)
                                                                                                                                    gCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGAt                           <  2:900106/139‑1 (MQ=255)
                                                                                                                                    gCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGAt                           >  1:900106/1‑139 (MQ=255)
                                                                                                                                       ggataaggTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGa              >  2:848287/7‑149 (MQ=255)
                                                                                                                                         ataaggTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGAtt            >  1:995028/5‑149 (MQ=255)
                                                                                                                                              gtgtAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGAtt                          >  1:1054119/1‑130 (MQ=255)
                                                                                                                                              gtgtAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGAtt                          <  2:1054119/130‑1 (MQ=255)
                                                                                                                                                 tAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACcgc    >  4:695277/1‑149 (MQ=255)
                                                                                                                                                   aaCACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                  <  4:557630/85‑1 (MQ=255)
                                                                                                                                                   aaCACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                  >  3:557630/1‑85 (MQ=255)
                                                                                                                                                   aaCACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGa  >  1:839516/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGA  >  NC_000913/2169518‑2169813

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